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As  physics  continues  to  expand  its  knowledge  base,  physicists  seek  new  frontiers 
to  investigate.  Quantum-mechanical,  two-dimensional  systems  have  proven  to  be  a 
subject  that  is  not  only  rich  in  new  discoveries  (e.g.,  Kousterlitz-Thoules  transitions  and 
new  phases  of  matter)  but  also  filled  with  exciting  predictions  (e.g.,  new  superfluids).  The 
physisorption  of  a  gas  onto  the  surface  of  a  homogenous,  spacious  (on  a  molecular  scale) 
substrate  with  a  low  adsorption  potential  is  one  physical  analogue  to  the  much-analyzed 
theoretical  two-dimensional  system.  Hydrogen  in  reduced  dimensions  has  been  found  to 
have  suppressed  melting  and  freezing  points,  thus  possibly  permitting  the  onset  of  Bose 
condensation  and  a  new  superfluid  phase.     The  use  of  boron  nitride  as  an  adsorption 
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substrate  allows  for  the  study  of  physisorbed  systems  with  a  lower  adsorption  potential 
than  previous  studies  using  similar  substrates  (i.e.,  graphite  and  magnesium  oxide). 

This  dissertation  has  two  parts.  The  first  concerns  the  translational  properties  and 
adsorption  energies  of  hydrogen  adsorbed  onto  boron  nitride.  These  properties  are 
investigated  through  the  use  of  volumetric  adsorption  isotherm  techniques.  The  data 
suggest  that  the  adsorption  of  hydrogen  occurs  in  a  step- wise  manner  at  temperatures 
below  20  Kelvin.  Changes  in  the  translational  phases  (vapor-hquid-solid)  occur  between 
10  to  20  Kelvin  for  the  first  four  adsorbed  monolayers.  Isotopic  effects  are  investigated 
through  the  use  of  the  three  common  forms  of  hydrogen:  molecular  hydrogen  (H2), 
deuterium  hydride  (HD),  and  deuterium  (D2).  The  critical  temperatures  of  the  second, 
third,  and  fourth  layers  are  determined,  presented,  and  compared  with  the  known  phase 
diagrams  of  hydrogen  isotopes  adsorbed  onto  graphite  and  MgO.  While  the  adsorption 
potential  of  the  hydrogen-boron  nitride  system  is  found  to  be  less  than  that  of  previously 
studied  substrates,  it  does  not  translate  into  lower  critical  temperatures. 

The  second  part  of  this  thesis  concerns  the  orientational  ordering  of  H2  in  two- 
dimensions  as  investigated  through  the  use  of  nuclear  magnetic  resonance  (NMR) 
techniques.  While  ordering  in  three-dimensional  (bulk)  systems  involves  a  crystal  phase 
change  and  a  strong  ortho-hydrogen  concentration  dependence,  ordering  of  hydrogen  on 
boron  nitride  is  shown  to  be  more  subtle.  Evidence  of  partial  ordering  below  1  Kelvin  is 
uncovered.  This  ordering  is  shown  to  depend  not  only  on  the  ortho-hydrogen 
concentration  but  also  on  the  number  of  adsorbed  layers. 
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CHAPTER  1 
INTRODUCTION 


Naturea  non  nisi  pareudo  vincitur. 
Nature,  to  be  commanded,  must  be  obeyed. 

Novum  Organum  1620,  Sir  Francis  Bacon  1561  -  1626 


1.1   Quantum  Mechanical  Nature  and  Other  Properties  of  Hydrogen 

The  low  temperature  properties  of  hydrogen  are  dominated  by  the  quantum 
mechanical  nature  of  the  molecule  as  evidenced  by  the  large  zero  point  motion,  quantum 
rotator  properties,  and  the  existence  of  two  distinct  molecular  species.  The  interaction 
between  adjacent  hydrogen  molecules  is  a  weak  Van  der  Waals  attraction  and  can  be 
modeled  by  the  Lennard-Jones  (6-12)  potential 


V{r)  =  -AE 


\r  ) 


6 


(of 
\rJ 


Equation  1-1 


where  e  and  a  are  the  well  depth  and  the  hard-core  radius,  respectively,  and  are  known  as 
the  Lennard-Jones  (L-J)  parameters.  The  quantum  nature  of  a  molecule  can  be  quantified 

by  rescaling  the  Schroedinger  equation  in  terms  of  the  L-J  parameters: 

r  V(r) 

r*  =  —    and    V*  =  —^  Equation  1  -2 

a  e 


to  derive  the  de  Boer  parameter 


yjlmsa 
which  describes  the  importance  of  the  quantum  mechanical  effects. 


Equation  1-3 


In  terms  of  this  fundamental  parameter  measuring  the  translational  degrees  of  freedom, 
hydrogen  is  second  only  to  helium  in  its  quantum  mechanical  nature. 

The  quantum-mechanical  nature  of  the  rotational  degrees  of  freedom  are  even 
more  pronounced.  The  homonuclear,  diatomic  forms  of  molecular  hydrogen,  H2,  D2,  and 
T2,  must  be  described  as  quantum  rotators  with  good  quantum  numbers,  J,  for  the  angular 
momentum.  In  accord  with  the  foundations  of  quantum  mechanics,  the  Pauli  exclusion 
principle  and  indistinguishability,  molecular  hydrogen  forms  two  distinct  molecular 
species,  ortho-hydrogen  and  para-hydrogen  as  outlined  in  Table  1-1.  The  ortho-para 
terminology  refers  to  the  spin  state  with  the  largest  number  of  occupants  at  low 
temperatures.  At  high  temperatures,  room  temperature  and  above,  gaseous  H2  converts 
into  three  parts  ortho  and  one  part  para  as  governed  by  spin  statistics.  This  nomenclature 
can  be  difficult  to  follow  for  the  uninitiated.  Some  authors  in  the  field,  P.  Soures 
[Sou  86],  for  example,  have  therefore  suggested  dropping  the  ortho-para  nomenclature  in 
favor  of  a  direct  referral  to  the  spin  status  of  the  species,  eg.,  "spin  1  hydrogen,"  and  "spin 
2  deuterium,"  etc.  Unfortunately,  the  new  names  have  not  been  universally  adapted. 

Table  1-1   Ground  States  of  the  Homonuclear,  Diatomic  Molecular  Hydrogen  Isotopes. 


Ortho  H2 
(and  T2) 

Para  D2 

Para  H2 
(and  T2) 

Ortho  D2 

Angular 
Momentum 

J=l 

(odd) 

J=  1 

(odd) 

J  =  0 
(even) 

J  =  0 
(even) 

Nuclear 
Spin 

1=  1 
(even) 

1=  1 
(odd) 

1  =  0 
(odd) 

1  =  0,2 
(even) 

In  contrast,  atomic  helium  (''He)  is  spinless  and  the  question  of  indistinguishability 
does  not  anise  for  heteronuclear  HD.  Selected  properties,  including  A,  o,  and  e  of  the 
hydrogen  isotopes  are  provided  at  the  end  of  this  chapter. 


Boron  Nitride 


Graphite 


Magnesium  Oxide 


^ 


m 


•nitrogen   ^  boron  ^  carbon 

Hexagonal 

Figure  1-1.  Crystal  structure  of  graphite,  BN,  and  MgO. 


Simple  Cubic 


1.2  Adsorption  Substrates 

Of  the  various  substrates  used  for  the  study  of  gases  adsorbed  onto  surfaces,  three 
are  of  special  concern  here:  graphite,  boron  nitride  (BN),  and  magnesium  oxide  (MgO). 
As  illustrated  in  Figure  1-1,  the  crystal  structure  of  both  graphite  and  boron  nitride  is 
based  on  a  planar  structure.  Furthermore,  the  lattice  dimensions  of  BN  are  only  about  5% 
greater  than  those  of  graphite.  These  dimensions  and  others  are  presented  in  Table  1-2. 
For  BN,  these  planes  are  stacked  into  platelets  of  the  order  of  0. 1-5.0  microns  in  diameter 
and  0.01-0.1  microns  thick.  Figure  1-2  is  a  scanning  electron  microscope  photograph'  of 
the  boron  nitride  used  m  this  study.  In  contrast  to  the  hexagonal  structure  of  graphite  and 
BN,  MgO  has  a  simple  cubic  crystal  structiue  with  a  base  distance  of  4.23  Angstroms. 
Thus,  with  the  use  of  just  these  three  substrates,  it  is  possible  to  investigate  the  effects  of 


Image  produced  at  the  University  of  Florida. 


Figure  1-2.     Photograph  of  BN  from  a  scanning  electron  microscope  image  in  the 
reflective  mode. 


Figure  1-3.    Photograph  of  BN  from  a  scanning  electron  microscope  in  the  transmission 
mode.  Note  the  planes  of  BN  which  form  the  platelet. 


binding  energy,   lattice  spacing,   corrugation,   and   substrate  lattice  symmetry  on  the 
translational  and  orientational  properties  of  adsorbed  films. 


Table  1-2.  Intermolecular  distances,  in  Angstroms,  for  adsorbed  monolayers  and  bulk 
hydrogen. 


Substrate 

Graphite 

BN 

MgO 

Bulk  HCP 

in-plane 

4.26  A 

4.35  A 

4.23  A 

H2:3.79A 

out-of-plane 

6.71  A 

6.66  A 

4.23  A 

D2:  3.61  A 

1.3  Molecular  Orientation 


1.3.1  Overview 


Nuclear  magnetic  resonance  (NMR)  studies  at  high  frequencies  are  invaluable  as 
an  unambiguous  probe  of  the  molecular  dynamics  and  molecular  order  parameters. 
Specifically,  the  orientational  ordering  of  a  system  of  spin  1/2  molecules  is  uniquely 
determined  fi'om  the  NMR  line  shape  of  that  system.  At  high  temperatures,  the  hydrogen 
molecule  rotates  rapidly  compared  to  the  interaction  time  causing  the  intramolecular 
dipole-dipole  interaction  to  average  to  zero.  At  lower  temperatures,  the  rotation  of  the 
molecules  slows  down  and  the  dipole-dipole  interaction  becomes  measureably  significant. 
At  some  low  temperature  the  rotation  of  the  molecule  is  severely  inhibited  and  the  solid 
undergoes  a  phase  transition  to  the  orientationally  ordered  phase.  This  orientational  phase 
transition  (near  1  K)  should  not  be  confijsed  with  the  translational  phase  transition  (the 
liquid-to-solid  transition)  which  occurs  at  a  temperature  approximately  one  order  higher  in 
magnitude.  The  dipole-dipole  interaction  is  central  to  the  determination  of  the  orientation 
of  the  molecular  axis.     This  interaction  lifts  the  rotational  degeneracy  of  the  nuclear 


Zeeman  ground  state  allowing  NMR  to  probe  the  orientation  of  the  molecule.    The  line 
shape  can  be  calculated  exactly  starting  from  the  full  Hamiltonian. 


1.3.2  Dipole-Dipole  Interaction 

The  intramolecular  interaction  between  two  hydrogen  atoms  within  the  hydrogen 
molecule  is  due  solely  to  the  nuclear  dipole-dipole  forces.  The  Hamiltonian  for  this  system 
of  two  spins  in  a  large  magnetic  field  can  be  written  as  the  sum  of  two  terms.  Ho  and  Hdd- 

,       y-fi-  ,  ,  Equation  1-4 

=-rriBoir,  +  i:)^^[l\-r-xr-n)irn)] 

where  /  is  the  proton  nuclear  spin  operator  and  the  angles  are  defined  in  Figure  1-4.  The 

first  term,  Ho,  represents  the  large  Zeeman  splitting  of  the  nuclear  energy  levels  caused  by 
the  static  magnetic  field  applied  in  the  Z  direction.  The  second  term,  Hdd,  represents  the 
dipole-dipole  interaction  and  results  in  a  perturbation  of  the  energy  levels  due  to  the  first 
term.  It  is  elucidative  to  express  the  dipolar  Hamiltonian  in  terms  of  spherical  harmonics, 
namely 

//^o=^^^{/  '■/  '-3{/  'cos0  +  sin0(/  '-0030  +  /  'sinO)]*  r,       ,•       ,   c 

DD       ^i     <.!      1         Liz  Wx  ly  'i     ^         Equation  1-5 

[/^^  cos0  +  sin0(//  -cosO+Z  "  -sinO)]} 

This  expression  can  be  further  transformed  by  combining  the  x  and  y  components 
into  raising  (  /^  )  and  lowering  (  /    )  operators, 

//^~{/'-/^-3[/;cos0  +  ls.n0(/;.e-+/_'.e-)l*  Equation  1-6 

[//  cos0  +  sin0(//  -e"^  + 1  ^  -e  "*)]} 
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Figure  1-4.  Axes  used  to  define  the  dipole-dipole  interaction. 


In  general,  the  dipolar  Hamiltonian  can  be  written  as  the  sum  of  six  terms, 


^DD   ~^A+^B+^C+^D+^E+^F 


Equation  1-7 


where, 


/-/!) 


^.=/z/zO-3cos^0) 

//,  =  l(l-3cos^0)(/^/;-/'/)=i(l-3cos^0)(/;/: 

-}  Equation  1-8 

H,  =  --smQcosQe-^{j[j:-j[jl)  =  H,* 

In  the  limit  that  the  Zeeman  field  splitting  is  much  larger  than  the  dipole-dipole  interaction 

(i.e..  Ho  »  Hdd),  then 

H^  ^  Hg,         H^,Hg  »  Hc,H^  ,  Equation  1-9 

and 

H^ , H^  »  H^,Hp      .  Equation  1  - 1 0 

Furthermore,  H;\^  represents  the  classical  interaction  between  two  static  local 
dipoles.   The  second  term,  Hg  represents  the  primary  resonance  at  the  Larmor  frequency 


and  allows  a  simultaneous  reversal  of  two  neighboring  spins  in  opposite  directions.  The 
last  four  small  terms  represent  off-diagonal  matrix  elements  and  admix  higher  states. 
Thus,  a  good  approximation  is  to  drop  Hq,  Hq,  Hg,  and  Hp  while  retaining  the  secular 
terms,  H^  and  H3. 

Two  eigenstates  exist  for  this  simplified  Hamiltonian.  Since  the  singlet  state 
represents  the  spin  zero  para-hydrogen  molecules,  we  can  define  this  state  to  have  zero 
Zeeman  energy.  Physically,  these  molecules  act  as  inert  dilutents  and  are  ignored  by  the 
NMR  but  do  provide  other  physical  effects.  The  triplet  term,  however,  provides  the  NMR 
signal  and  can  be  most  easily  understood  by  changing  the  spin  variables  into  spin  up  and 
spin  down  operators.  By  defining  the  total  operators. 


r.(i-3cos^^),  fr  =  fjl+^(i:i:+i[i:) 

it  can  be  shown  that. 


Equation  1-11 


The  eigenstates  for  this  Hamiltonian  include  the  singlet  state. 


|o) 


^ 


and  the  triplet  state, 


|i)  =  k+);   |o) 


k-)+l-+) 


■i>=l-->. 


Equation 


-17. 


Equation  1-13 


Equation  1-14 


The  energies  of  each  level  can  be  determined  as 
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£-,-HW-i>=  r^^o^^T 


,2t2 


^o      \  \    \  I  2r'  Equation  1-14 
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Figure  1-5.  Effect  of  dipolar  interaction  on  Zeeman  splitting  of  spin  1/2  energy  levels. 


Therefore,  two  resonant  frequencies  exist,  co'  and  co",  where. 


.2*2 


rio}'  =  E-rEo  =  rf^Bo+-^r 


hco'=E,-E,  =  Yt^Ro-^7:^'^ 


3rlr 
4r' 


Equation  1-16 


The  form  of  continuous  wave  (CW)  NMR  chosen  for  use  utilizes  a  fixed  frequency  and  a 

variable  magnetic  field,  B,  to  measure  the  resonant  absorption.   In  this  case  the  resonances 

occur  at 

B  =  Bo±ar  Equation  1-17 


11 

where  g  =  ^ ^  and   5  =^  .   Thus,  for  a  pair  of  dipoles,  or  a  single  crystal  of  such  pairs, 

Ar  y 

the  NMR  line  shape  becomes  two  lines  centered  at  Bq  with  a  peak-to-peak  separation  of 
2ar. 

For  a  powdered  sample,  many  pairs  of  dipoles  are  oriented  randomly,  and  an 
average  over  all  possible  angles,  a  "powder  average",  must  be  performed.  Such  an 
average  yields  the  familiar  Pake  Doublet  Powder  Pattern  NMR  line  shape. 


1.3.3  The  Pake  Doublet  Line  Shape 

The  Pake  doublet  is  a  very  unique  line  shape  that  can  be  calculated  exactly.  The 
determination  of  the  line  shape  depends  upon  finding  the  intensity  of  each  frequency  given 
a  distribution  of  spins  (frequencies).  For  a  polycrystalline  (powdered)  sample,  each  crystal 
in  the  powder  yields  two  lines  due  to  the  dipole-dipole  interaction.  The  frequency 
separation  of  these  lines  from  the  resonance  center  is  given  by 

v  =  ^^— =  ±a'(3cos-0-l)  .  Equation  1-18 

where   a'  =y allii. 

The  intensity  of  the  signal  is  a  fijnction  of  the  frequency  of  applied  radiation  and  is 

related  to  the  angle  between  the  intramolecular  axis  and  the  externally  applied  field  (e.g., 

the  Zeeman  axis) 

I[v)  dv  =  P{(d)dQ.  .  Equation  1-19 

Since  the  interaction  is  axially  symmetric,  the  differential  solid  angle  depends  only  upon 

the  azimuthial  angle,  i.e., 

dO.  =  Ik  ^/(cosG)  .  Equation  1-20 

For  a  powdered  sample,  all  angles  are  equally  probable,  i.e., 

P(0)  =  — .  Equation  1-21 
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Thus,  the  intensity  of  the  signal  can  be  found  as  a  function  of  the  frequency,  i.e., 

1  1 


/(.) 


idv, 


i/(cos0) 


Equation  1-22 


-2a 


-a 


a 


2a 


Figure  1-6.  Pake  doublet  NMR  line  shape  for  a  powder  sample. 


The  desired  differential  can  be  calculated  from  Equation  1-18,  i.e., 

dv 


d{cosQ) 


Sa'cosQ 


Equation  1-23 


Thus,  the  intensity  diverges  at  two  specific  frequencies,    y-±a',  i.e., 

1 


/(v). 


2^/v± 


a' 


where 


2a  >  v> -a 


n 


Equation  1-24 


Equation  1-25 


as   0  <  0  <  — .  The  "powder  average"  Pake  Doublet  line  shape  is  presented  in 
Figure  1-6. 
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1.4  Orientational  Ordering  in  Hydrogen  and  Deuterium. 

1.4.1  Overview. 

Orientational  ordering  may  occur  for  any  system  of  atoms  or  molecules  which  have 
anisotropic  intermolecular  interactions  (e.g.,  electro-  or  magneto-static  potentials  of 
dipolar,  quadrupolar,  octopolar,  or  higher  order).  In  such  a  system,  when  the  energy  of 
the  random  motion  associated  with  the  temperature  of  the  lattice  decreases  below  the 
interaction  energy  due  to  the  static  potentials,  the  molecules  may  then  align  themselves 
with  respect  to  their  neighbors.  For  the  case  of  quadrupoles,  as  in  this  study,  the  lowest 
energy  state  is  a  "tee"  formation  between  pairs.  Complete  orientational  ordering  occurs  if 
this  alignment  is  both  long-term  and  long-ranged.  Some  systems  may  show  only  partial 
ordering  for  only  a  short  period  of  time  and  have  a  wide  distribution  of  local  order 
parameters.  Such  systems  are  said  to  reside  in  an  orientational  glass  state. 

The  properties  of  bulk  solid  hydrogen  at  low  temperatures  have  been  studied  for 
some  time.  For  example,  the  static  susceptibility  of  solid  hydrogen  at  2  K  was  measured 
as  early  as  1937  [Laz  37].  The  low  temperature  orientational  properties  of  hydrogen  have 
been  identified  and  investigated  almost  continuously  since  then.  Specific  heat 
measurements  identified  a  singularity  in  high  ortho-hydrogen  samples  indicating  the 
existence  of  either  a  transitional  or  orientational  transition  [Hil  54,  Ahl  64].  X-ray 
diffraction  revealed  a  change  in  the  crystal  lattice  spacing  at  a  similar  low  temperature 
[Mil  66]  as  did  neutron  diffraction  [Muc  66]  and  infrared  (IR)  absorption  spectra  studies. 
Beginning  in  the  early  1970's,  NMR  studies  identified  the  various  orientational  phases  and 
produced  a  phase  diagram  [Sul  72].  Most  bulk  properties  of  hydrogen  may  be  found  in  a 
fairly  comprehensive  review  article  by  Silvera  [Sil  80]  as  well  as  in  a  book  by  P.  Souers 
[Sou  86]. 
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The  orientational  ordering  of  hydrogen  and  deuterium  in  part  determines  the 
crystal  structure  of  the  bulk  solid.  If  no  orientational  ordering  is  present,  as  is  the  case  for 
a  system  of  pure  J=0  molecules  (para-hydrogen  and  ortho-deuterium),  then  the 
constituents  form  a  hexagonal  close  packed  (h.c.p.)  crystal  lattice  structure.  However, 
this  triangular  lattice  structure  is  incommensurate  with  the  cubic  symmetry  of  quadrupolar 
interactions  in  J=l  molecules  (ortho-hydrogen  and  para-deuterium).  Thus  orientational 
ordering  of  samples  with  J=l  molecules  drives  the  lattice  into  a  cubic  structure.  For  pure 
J=l  systems  this  occurs  at  approximately  3  K  for  hydrogen  and  4  K  for  deuterium  and  the 
resulting  lattice  structure  is  three  interpenetrating  cubic  lattices  each  with  orientational 
order  along  the  three  unique  body  diagonals  [Sou  86,  Sil  80].  For  the  more  general  case 
of  mixed  J  systems  the  transition  temperature  drops  with  decreasing  J=l  concentration. 
Below  55%  ortho-hydrogen  concentration,  the  ordering  is  in  the  form  of  a  quadrupolar 
glass  [Sul  72]. 


1.4.2  Review  of  Three-Dimensional  Results 

1.4. 2. a  Order  Parameters 

Bulk  ortho-hydrogen,  a  quantum-mechanical  molecule  with  J=l,  has  a  relatively 
large  electric  quadrupole  moment  and  exhibits  orientational  ordering  below  4  K.  The  low 
temperature  behavior  is  dominated  by  the  anisotropic  quadrupolar  interaction,  given  by 

Vqo  =  [^    rgF(Q,,n,)  Equation  1-26 

where  A  is  the  lattice  parameter,  R  is  the  intermolecular  distance,  F  is  an  angular  function 

describing  the  anisotropy  of  the  interaction,  and  Fq  is  the  quadrupolar  coupling  constant 

r„  = K-  ^  \K.  Equation  1-27 

^       25a' 
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Here  eQ  is  the  quadrupole  moment  of  the  molecule.  The  quadaipole-quadaipole 
interaction  can  be  expressed  either  in  terms  of  the  spherical  harmonics  or  in  terms  of 
spherical  tensor  operators.  The  components  of  the  second  rank  quadrupolar  tensor, 

Q     =i-J-5     --[j  J    +J  j)j  Equation  1-28 

define  the  parameters  which  describe  the  degrees  of  freedom  for  the  system.  Of  the  five 
unique  components,  three  are  Euler  angles  that  specify  the  orientations  of  the  principal 
axes  of  the  tensor  Qab,  and  two  describe  the  degree  of  ordering.  The  latter  two  are  of 
special  interest.  They  are  the  eccentricity, 

77  =  (j]  -  J]  )  Equation  1  -29 

which  is  zero  for  an  axially  symmetric  body  (a  good  approximation  for  hydrogen)  and  the 

alignment, 

o-  =  2„  =  -  -  {2  J]  -  2) .  Equation  1  -30 

The  alignment  varies  from  one,  for  Jz  =  0,  to  minus  one-half,  for  Jz  =  +1  or  -1. 
Since  the  nuclear  spin  is  uniquely  related  to  the  angular  momentum,  the  orientational 
degrees  of  freedom  can  be  studied  using  NMR  techniques  that  probe  the  nuclear  spin  of 
the  ortho-hydrogen  molecule. 

The  effect  on  any  one  hydrogen  spin  in  a  pure  molecular  hydrogen  system  subject 
to  a  large  magnetic  field  can  be  readily  understood  in  terms  of  the  order  parameters.  The 
energy  levels  of  each  atomic  spin  are  affected  by  three  forces:  the  external  Zeeman 
magnetic  field,  the  intramolecular  forces  due  to  the  other  proton  within  the  hydrogen 
molecule,  and  the  sum  of  all  of  the  intermolecular  forces.  The  intermolecular  forces  on 
any  one  spin  are  much  smaller  than  the  intramolecular  forces  since  the  dipole-dipole 
interaction  decays  as  the  cube  of  the  inverse  distance  between  particles  and  the  spacing 
between  molecules  is  about  four  times  the  intramolecular  distance,  i.e.. 
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^aY    fO.lSAY 


\.RJ 


J.ISAJ 


=  0.8%. 


Equation  1-31 


Of  the  two  remaining  interactions,  the  Zeeman  field  can  be  made  arbitrarily  large 

relative  to  the  intramolecular  interaction.     At  high  resonance  frequencies  (tens  and 

hundreds  of  Megahertz)  the  Zeeman  interaction  dominates.  Thus,  as  described  above,  the 

resulting  Hamiltonian  can  be  broken  into  two  parts: 

Equation  1-32 


H  =  H+H 


DD 


where  the  dipolar  term  acts  as  a  small  correction  to  the  Zeeman  term  and  splits  the  single 
resonance,  v^  -y  hB^,  into  two  resonances  as  mentioned  in  Section  1.3.2, 


y    =Vo±Vrf 

Equation  1-33 

where  [Sul  78], 

'      2 

a  (3 cos'  0  - 1)  +  -  77  sin^  0  cos^  0 

Equation  1-34 

and 

d 

1  rrv' 
-c.    /3V=54.2/r//z. 

Equation  1-35 

Note  that  not  only  the  alignment,  s  ,  but  also  the  eccentricity,  r\,  has  entered  into 
the  line  shape  and  that  the  exact  value  of  the  dipolar  constant,  d,  depends  on  the 
translational  and  rotational  properties  of  the  hydrogen  system  under  examination. 

Ortho-hydrogen  (J=l)  decays  into  its  para-hydrogen  (J=0)  ground  state  at  low 
temperatures  [Sil  80].  In  contradistinction  to  ortho-hydrogen,  para-hydrogen  has  no  spin 
and  no  higher  moments,  including  quadrupolar.  Thus,  this  J=0  hydrogen  does  not  yield  an 
NMR  signal  and  behaves  as  an  inert  diluent  to  the  system  of  interacting  J-1  hydrogens.  A 
pure  ortho-hydrogen  system  eventually  decays  into  a  disordered  system  and  the  effects  of 
disorder  on  a  fixed  sample  can  be  studied  by  simply  aging  the  sample.  A  partial  list  of 
some  of  the  other  properties  of  hydrogen  is  provided  at  the  end  of  this  chapter. 
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1.4.2.b.  Phase  EHagram 

Continuous  wave  nuclear  magnetic  resonance  (CW  NMR)  experiments  on  three- 
dimensional  (3D  or  bulk)  systems  have  been  performed  for  some 
time[Sil  80,  Sul  78,  Har  85].  It  was  found  that  the  phase  diagram  of  bulk,  solid  hydrogen 
includes  three  distinct  orientational  phases  below  4  K,  as  illustrated  in  Figure  1-7.  The 
appearance  of  the  nuclear  magnetic  frequency  response,  the  NMR  "hne  shape,"  depends 
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Figure  1-7.  Orientational  ordering  phase  diagram  for  bulk  hydrogen. 
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upon  the  crystalline  state  of  the  sample.  Both  polycrystalline  and  single  crystal  samples 
have  been  studied  and  much  literature  exists  describing  the  various  subtle  properties. 
However,  since  the  two-dimensional  samples  more  closely  parallel  the  three-dimensional 
polycrystalline  samples,  polycrystalline  line  shapes  will  primarily  be  mentioned. 

1.4.2.C  NMR  Line  Shapes 

Solid  hydrogen  above  3  K  shows  no  signs  of  ordering.  The  NMR  line  shape  is 
fairly  narrow,  approximately  50  kHz,  and  relatively  independent  of  ortho-hydrogen 
concentration  as  illustrated  in  Figure  1-8  (a).  Part  (b)  displays  the  distribution  of  the 
alignment  where  unity  represents  complete  ordering.  For  such  a  narrow,  disordered  line 
shape,  the  distribution  resembles  a  delta  function  at  zero  order  parameter.  The 
eccentricity  is  nearly  zero  for  all  systems  of  hydrogen  molecules  discussed  here. 

A  second  phase  exhibits  long  range  ordering  and  occurs  for  high  ortho-hydrogen 
concentrations  (X  >  55%)  and  low  temperatures  (T  <  3  K).  The  resulting  line  shape  is  the 
familiar  Pake  Doublet,  shown  for  a  polycrystalline  sample  in  Figure  l-8(c),  with  a 
characteristic  width  of  3D/2,  or  approximately  150  kHz.  The  line  shape  for  a  single 
crystal  sample  (composed  of  four  sublattices)  displays  two  lines  per  sublattice  whose 
frequency  separation  depends  upon  the  angle  between  the  local  crystal  axis  and  the  applied 
magnetic  field.  The  angular  dependent  splitting  has  been  studied  in  detail  [Har  85]  and 
magic  angle  spinning  used  to  contract  the  four  pairs  of  lines  into  one  set.  The  distribution 
of  order  parameters  resembles  a  delta  function  near  unity,  as  presented  in  Figure  l-8(d). 

The  third  and  most  interesting  phase,  the  glass  phase,  has  short-ranged  ordering 
for  intermediate  ortho-hydrogen  concentrations  (20  <  X  <  55%)  and  low  temperatures 
(T  <  IK).  The  CW  NMR  line  shape  depends  strongly  upon  both  the  temperature  and  the 
ortho-hydrogen  concentration.   Accordingly,  there  exists  a  broad  distribution  of  the  local 


NMR  LINE  SHAPE 


(a) 


-3d       -3d/2      0        3d/2       3d 


(c) 


-3d         -3d/2      0        3d/2       3d 


(e) 


-3d         -3d/2      0        3d/2       3d 


(g) 


19 
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Figure  1-8.  NMR  line  shapes  and  associated  distributions  of  the  order  parameter,  a,  for 
selected  hydrogen  systems. 
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order  parameters  and  the  shape  of  this  distribution  is  highly  temperature  and  ortho- 
concentration  dependent.  An  example  of  such  a  distribution  and  its  associated  line  shape 
can  be  seen  in  Figure  1-8  (f)  and  (e),  respectively.  Parts  (g)  and  (h)  of  Figure  1-8  are 
included  here  for  convenience  but  will  be  discussed  in  Chapter  4. 

Figure  1-9  (top)  illustrates  the  variation  of  the  line  shape  corresponding  to  the 
transition  from  the  ordered  to  disordered  phase.  Part  (a)  displays  the  ordered  Pake 
Doublet  while  part  (b)  shows  the  orientationally  disordered  solid  NMR  line  shape.  The 
mean  order  parameter,  shown  in  part  (c),  exhibits  a  first  order  phase  transition  as  it  rapidly 
changes  from  near  unity  to  zero. 

Figure  1-9  (bottom)  illustrates  the  spectral  change  for  the  transition  from  the 
quadrupolar  glass  to  the  unordered  phase  for  solid  hydrogen.  Part  (a)  displays  the  glass 
CW  NMR  line  shape  while  part  (b)  shows  the  unordered  line  shape.  In  this  case,  the 
phase  change  is  weakly  second  order  as  the  mean  order  parameter  undergoes  a  more 
subtle  transition,  as  displayed  in  part  (c). 

This  new  collective  orientational  phase,  the  quadrupolar  glass  phase,  was  an 
important,  exciting,  and  unexpected  discovery.  The  intermolecular  interactions  of  the 
hydrogen  system  are  relatively  simple  and  well  understood.  Further,  NMR  provides  for  a 
very  fundamental  probe  of  the  glass  state.  Thus,  this  physically  simple,  yet  theoretically 
rich,  system  has  been  thoroughly  examined  and  modeled.  Since  so  much  is  known  about 
the  three-dimensional  system,  the  effects  of  the  reduction  in  the  dimensionality  of  the 
system  can  be  readily  studied.  Hydrogen  adsorbed  onto  a  surface  provides  a  model  two- 
dimensional  system  for  this  purpose.  Knowledge  gained  from  this  relatively  simple  two- 
dimensional  system  may  aid  in  the  understanding  of  more  complicated  systems  in  reduced 
dimensions. 

A  second  phase  exhibits  long  range  ordering  and  occurs  at  high  ortho-hydrogen 
concentrations  (X  >  55%)  and  low  temperatures  (T  <  3  K).  The  resulting  line  shape  is  the 
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Figure  1-9.  Temperature  dependence  of  the  CW  NMR  line  shape  and  mean  order 
parameter  for  the  long-ranged  ordered  phase  (top)  and  the  glass  phase  (bottom)  for  bulk 
molecular  hydrogen. 
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familiar  Pake  Doublet  presented  for  a  poiycrystalline  sample  in  Figure  l-8(c),  with  a 
characteristic  width  of  3d/2  (approximately  150  kHz). 

1.4.3  Two-Dimensional  Orientational  Ordering  Model 

Figure  1-10  illustrates  one  model  of  orientational  ordering  for  a  simple  two- 
dimensional  system  of  quadrupoles.  Part  (a)  shows  four  quadrupoles  confined  to  a  square 
lattice.  In  this  case,  the  symmetry  of  the  underlying  lattice  matches  the  symmetry  of  the 
intermolecular  quadrupolar  interaction.  At  zero  temperature  the  system  relaxes  into  its 
lowest  energy  state  and  each  molecule  forms  a  "Tee"  configuration  with  each  and  every 
neighboring  quadrupole.  Part  (b)  shows  a  system  of  quadrupoles  confined  to  a  hexagonal 
(triangular)  lattice.  Such  a  lattice  symmetry  is  incompatible  with  the  intermolecular 
interaction.  At  zero  temperature,  the  system  is  frustrated  because  each  molecule  in  the 
lattice  cannot  attain  its  lowest  energy  state  with  respect  to  all  of  its  neighbors.  The  lowest 
energy  for  the  entire  system  may  be  a  herringbone  or  a  pinwheel  pattern.  Parts  (c)  and  (d) 
show  the  effect  of  the  addition  of  a  vacancy  or  a  neutral  molecule(e.g.,  a  noble  gas  or 
para-hydrogen).  Since  a  neutral  molecule  does  not  orientationally  interact  with  any  of  its 
neighbors,  these  molecules  behave  as  inert  diluents;  that  is,  they  act  to  dilute  the 
orientational  interactions  among  the  molecules.  Since  all  of  the  lattice  sites  are  no  longer 
equivalent,  the  system  is  said  to  be  disordered.  The  effect  of  the  inert  diluents  to  a  system 
on  the  square  lattice  of  part  (c)  is  negligible.  However,  the  effect  of  the  substitution  of 
inert  diluents  on  the  triangular  lattice  of  quadrupoles  is  quite  noticeable.  In  this  case,  the 
remaining  quadrupoles  may  align  more  completely  with  respect  to  each  other.  In  other 
words,  the  degree  of  fiustration  decreases  as  the  disorder  increases  for  a  non-square 
lattice. 
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Two  definitions  require  extra  clarity:  zero  point  motion  (ZPM)  and  disorder. 
Zero  point  motion  is  the  intrinsic  quantum  energy  associated  with  all  matter.  While  it  is 
quite  small  for  massive  objects,  it  is  a  very  important  property  for  hydrogen.  ZPM  exists 
for  each  degree  of  freedom  including  rotational,  translational,  and  vibrational  motion. 
Likewise,  the  term  disorder  has  multiple  uses  in  the  literature.  Disorder  describes  a  system 
which  has  substitutional  impurity  as  well  as  a  lack  of  a  type  of  ordering  (e.g.,  rotational  or 
orientational). 


1.5  Studies  of  Adsorbed  Hydrogen 

A  distinction  may  be  made  in  the  manner  in  which  a  gas  is  adsorbed  onto  a  surface. 
In  general,  two  modes  exist  for  the  adsorption  of  hydrogen  onto  a  surface:  physisorbtion 
where  the  adsorption  potential  is  rather  weak  and  the  adsorbed  hydrogen  retains  its 
molecular  form,  and  chemisorption  where  the  adsorption  potential  is  rather  strong  and  the 
adsorbed  hydrogen  dissociates  into  its  atomic  form.  As  will  be  discussed  later  in  Chapter 
2,  hydrogen  is  physisorbed  onto  boron  nitride. 

While  no  studies  have  been  previously  reported  concerning  molecular  hydrogen  adsorbed 
onto  boron  nitride,  many  studies  have  been  reported  for  hydrogen  physisorbed  onto  other 
substrates.   A  history  of  the  studies  of  the  hydrogen  isotopes  on  graphite  and 

magnesium  oxide  has  been  compiled  recently  elsewhere  [Liu  93].  The  following 
overview  is  intended  to  give  more  of  a  flavor  of  the  extent  and  types  of  the  studies  that  are 
particularly  related  to  this  work. 

Graphite  has  been  used  as  a  substrate  for  decades  and  many  types  of  experiments 
have  been  conducted  on  adsorbed  gases,  including  low  energy  electron  diffraction  (LEED) 
and  Auger  electron  spectroscopy  (LEED-Auger)  [Suz  77],  x-ray  scattering  [Ros  83],  and 
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Figure   1-10.      Two-dimensional  model  of  quadrupolar  orientational   ordering. 
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ellipsometry  [Vol  89].  Adsoqjtion  isotherms  with  graphite  as  the  substrate  have  been 
carried  out  for  most  noble  gases  [Tho  70,  Goo  93],  the  hydrogen  isotopes  [Vil  92],  and 
some  hydrocarbons  [Alk  92,  Ham  83].  NMR  experiments  have  been  few  but  at  least 
three  adsorbed  gases  have  been  studied  at  some  level:  "^He  [Rol  72,  Cow  84],  H2  and  D2 
[Kub  78],  and  '^2  [Sul  83].  Of  particular  note  is  the  fact  that  the  translational  modes  of 
hydrogen  adsorbed  onto  graphite  have  been  well  studied.  Methods  of  study  include  heat 
capacity  [Mot  85,  Fre  85],  inelastic  neutron  diffraction  [Fre  90],  adsorption  isotherms 
[Liu  93],  and  computer  models  [Got  90].  Figure  1-1 1  presents  the  phase  diagram  of  HD 
on  graphite  as  recently  reported  [Liu  93]. 

Magnesium  oxide  is  another  substrate  of  interest.  Crystals  of  MgO  form  a  cubic 
structure  which  has  the  same  underlying  symmetry  as  the  quadrupole-quadrupole 
interaction.  NMR  studies  of  HD  and  D2  adsorbed  onto  MgO  have  been  performed 
[Jeo  93]  and  orientational  ordering  identified  at  elevated  temperatures  [Jeo  93].  It  is 
believed  that  the  strong  corrugation  of  the  surface  overwhelms  the  adatom-adatom 
interaction  and  is  responsible  for  the  orientational  ordering.  While  an  interesting  finding  in 
itself,  it  indicates  that  MgO  is  not  a  good  surface  for  the  study  of  frustration  and  disorder 
in  orientationally  ordered  systems  in  two-dimensions.  Adsorption  isotherms  have  also 
been  carried  out  to  study  the  two-dimensional  translational  modes  [Ma  89]  identifying  the 
liquid-solid  transition. 

1.6  Orientational  Ordering  of  Adsorbed  Quadrupoles 

A  system  of  adsorbed  nitrogen  is  the  classical  analogue  of  our  adsorbed  hydrogen 
system.  Both  gases  are  homonuclear  and  diatomic  and,  therefore,  molecular  quadrupoles. 
Both  gases  have  been  shown  to  form  a  registered  solid  phase  on  graphite,  in  particular, 
one  molecule  occupying  every  third  carbon  ring,  i.e.,  a  (  ^  x  V3  )  R  30°  triangular  net. 

While  the  extremely  light  mass  of  the  hydrogens  makes  quantum  mechanical  behavior  very 
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important,  the  large  mass  of  molecular  nitrogen  yields  mainly  classical  behavior.  Two 
common  forms  of  nitrogen  exist:  ''*N2  which  has  no  nuclear  dipole  moment  but  a  rather 
large  nuclear  quadrupole  moment  and  '^N2  for  which  the  nuclear  quadrupole  moment  is 
zero  but  the  nuclear  dipole  moment  is  non-zero.  Like  hydrogen,  the  nuclei  of  '^N2  are 
identical  fermions  and  thus  there  exist  two  molecular  species:  ortho-nitrogen  (total  nuclear 
spin  I  =  1  and  odd  orbital  symmetry)  and  para-  nitrogen  (1  =  0  and  even  orbital 
symmetry).  Above  4  K  the  two  species  are  randomly  mixed  together  in  a  constant  3:1 
ortho-  to  para-  nitrogen  ratio.  This  is  in  contrast  to  hydrogen  for  which  the  ortho  to  para 
ratio  strongly  depends  upon  the  temperature  with  about  equal  mixtures  at  77  K. 

Nitrogen  adsorbed  onto  graphite  has  been  studied  for  some  time  both 
experimentally  and  theoretically.  Low  energy  electron  diffraction  (LEED)  [Die  82],  heat 
of  adsorption  [Pip  83],  adsorption  isotherms  [Liu  93],  heat  capacity  [Cha  84,  Chu  77],  and 
NMR  [Sul  83]  are  some  of  the  various  experimental  methods  used  to  study  this  rich 
system.  The  results  from  these  studies  confirm  that  the  orientational  ordering  transition 
for  the  two-dimensional  system  is  reduced  to  28  K  from  the  bulk  temperature  of  35.6  K 
[Chu  77].  The  NMR  study  directly  probed  the  local  order  parameter  and  determined  that 
the  order-disorder  transition  is  of  first  order.  A  number  of  theoretical  methods  have  been 
used  to  study  the  orientational  properties  of  N2  and  H2  adsorbed  onto  graphite  including 
classical  Monte  Carlo  simulations,  molecular  dynamics  simulations,  renormalization  group 
analysis,  mean-field-like  approximations,  universality  arguments,  ground-state  energy 
investigations,  quasiharmonic  lattice  dynamics,  time-dependent  Hartree  calculations,  and 
desity-fijnctional  theory  of  freezing  [Mar  93]. 

An  analogous  reduction  in  the  temperature  for  the  orientational  ordering  of  the 
hydrogen  isotopes  (2.8  and  3.8  K  for  ortho-H2  and  para-D2  in  bulk  form,  respectively)  is 
expected  for  a  two-dimensional  system  of  hydrogen.  One  study  [Har  79]  modeled  the 
adsorption  of  H2  on  graphite  in  the  mean  field  approximation  using  a  two-dimensional 
system  of  interacting  quadrupoles  in  the  presence  of  the  substrate  field.    Herringbone  and 
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other  orientational  staictures  were  predicted  depending  upon  the  strength  of  the  substrate 
field  relative  to  the  quadrupole-quadrupole  interaction.  Another  study  [Kub  85]  used 
NMR  to  measure  the  orientational  ordering  of  both  hydrogen  and  deuterium  adsorbed 
onto  graphite.  For  90%  ortho-H2,  the  authors  observed  a  splitting  of  the  NMR  absorption 
line  shape  near  0.6  K  and  the  addition  of  line  structure.  However,  for  90%  para-D2  (the 
I  =  1  and  odd  orbital  symmetry  species  for  deuterium),  the  authors  found  no  evidence  of 
orientational  ordering  down  to  0.3  K.  Given  that  the  two  isotopes  orientationally  order  in 
the  bulk  at  similar  temperatures,  the  very  different  behavior  was  quite  unexpected. 

1.7  Liquid  Hydrogen:  Supercooling  and  Superfluidity 

Bose  condensation  of  liquid  helium  into  the  superfluid  state  has  been  observed  and 
studied  for  some  time.  That  helium  alone  should  posses  superfluidity  is  doubted  by  many 
and  a  current  search  exists  for  other  superfluids. 

The  transition  temperature  for  an  ideal  Bose  gas  can  be  calculated  [Lan  69]  and  is 
given  by: 


h' 


Mk, 


n 


T\   -  ^^^  ~  Equation  1-36 


\gJ 


where  n  is  the  concentration,  M  the  molecular  mass,  and  g  is  the  degeneracy  for  a  single 
particle  state  (1  for  para-  and  9  for  ortho-hydrogen).  In  bulk  form,  the  superfliud 
transition  of  helium  occurs  at  2.17  K,  approximately  2  degrees  below  condensation,  and 
near  the  predicted  temperature  of  3  K. 

Hydrogen  is  the  lightest  of  all  elements.  While  it  has  large  zero  point  motion,  the 
intermolecular  interactions  are  relatively  weak,  as  can  be  seen  from  Table  1-6.  Thus, 
molecular  para-hydrogen  is  a  potential  candidate  for  superfluidity  if  the  right  conditions 


28 


can  be  achieved.  Calculations  have  been  made  by  Giozburg  and  Sobyanin  for  para- 
hydrogen,  and  a  superfluid  transition  temperature  of  6.5  K  reported  [Gin  72].  It  is 
reasonable  to  expect  this  calculation  to  be  an  upper  limit  on  the  actual  temperature  since 
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Figure    1-11.       Translational    phase    diagram    of  HD    on    graphite,    from    [Liu  93]. 
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the  same  formula  overestimates  the  helium  transition  temperature.  Maris  et  al.  have 
argued  that  the  onset  temperature  for  para-hydrogen  should  be  lower  than  4.5  K  [Mar  83]. 

However,  the  mehing  temperature  for  hydrogen  is  approximately  14  K  and  this 
precludes  the  possibility  of  a  superfluid  transition  in  normal  bulk  hydrogen.  Thus  studies 
have  been  performed  to  try  to  supercool  hydrogen  under  various  conditions.  For  example, 
heat  capacity  studies  of  hydrogen  in  vycor  have  been  performed  by  Tell  et  al.  and 
supercooling  to  9.9  K  reported  [Tel  83]. 

In  two  dimensions,  a  Kosterlitz-Thoules  transition  (a  second  order  transition;  no 
latent  heat  of  transformation)  for  helium  films  has  been  found  and  an  onset  temperature 
given  by: 

^'=f  a7^"'  Equation  1-37 

where  n'  is  the  superfliud  film  density  at  the  onset  temperature.  If  hydrogen  films  were  to 
follow  a  similar  behavior,  then  an  estimate  can  be  made  from  this  equation  for  hydrogen 
films.  One  such  estimate  by  O.  Vilches  for  a  monolayer  of  hydrogen  with  a  film  density  of 
0.08  A'^  places  the  transition  temperature  around  2  K  [Vil  92].  The  results  of  current 
investigations  of  helium  and  hydrogen  on  surfaces  will  be  usefiil  to  determine  the 
conditions  for  which  to  study  hydrogen  in  reduced  dimensions. 

NMR  is  an  ideal  instrument  for  the  study  of  translational  phase  transitions.  The 
width  of  the  resonance  is  directly  related  to  the  motion  of  the  nuclei  under  study  and  can 
be  measured  directly  by  CW  NMR.  Pulsed  NMR  can  also  be  used  to  measure  rates  of 
motion.  For  example,  NMR  has  been  used  to  study  the  supercooling  of  hydrogen  in  the 
pores  of  zeolite  [Ral  91].  Rail  measured  both  the  transverse  relaxation  time  (T2)  and  the 
spin  echo  amplitude  as  fianctions  of  temperature.  Both  measures  were  found  to  decrease 
suddenly  at  about  10  K.  The  temperature  depended  upon  the  para-hydrogen 
concentration    of   the    sample    with    higher    concentrations    having    lower    transition 
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temperatures.  From  an  extrapolation  of  the  data.  Rail  estimated  that  100%  para-hydrogen 
could  be  supercooled  to  approximately  6  K. 

Likewise,  the  supercooling  of  hydrogen  on  the  surface  of  boron  nitride  can  be 
studied  by  NMR.  However,  both  the  low  number  of  atoms  under  study  in  two  dimensions 
and  the  background  atomic  hydrogen  signal  precluded  using  the  apparatus  as  set  up  for 
this  thesis  for  temperatures  above  approximately  7  K.  The  richness  of  this  field  and  the 
importance  of  the  possible  results  suggests  that  the  apparatus  should  be  improved  and  the 
high  temperature  studies  performed.  Such  improvements  are  underway. 

1.8  Selected  Properties  of  the  Hydrogen  Isotopes  at  Low  Temperatures 

This  section  is  provided  to  familiarize  the  interested  reader  with  some  of  the 
relevant  properties  of  molecular  hydrogen.  The  properties  of  hydrogen  at  low 
temperatures  has  been  studied  extensively  since  the  use  of  bulk  liquid  hydrogen  in  bubble 
chambers  of  high  energy  physics  experiments  began  some  years  ago.  The  three-dimesional 
properties  of  the  hydrogen  family  of  isotopes  listed  here  may  be  compared  with  the  two- 
dimensional  properties  discussed  elsewhere  in  the  text. 


Table  1-3.  Heats  of  transformation  of  selected  hydrogen  isotopes  at  the  triple  points 
[Sou  86]. 


Isotope 

Transformation 

H2                     HD             D2 

Sublimation 

1.028  kJ/mole 

1.270 

1.470 

Vaporization 

0.911 

1.110 

1.270 

Fusion 

0.117 

0.159 

0.199 
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Table  1-4.  Specific  volumes  of  selected  hydrogen  isotopes  at  the  triple  points  [Sou  86]. 


Hz 

HD 

D2 

Solid 

23.14*10"^  mVmole 

21.80*10-* 

20.40*10"* 

Liquid 

26.11*10^ 

24.61*10^ 

23.16*10-* 

Vapor 

0.0158 

0.0194 

0.00885 

Table  1-5.  Triple  and  critical  points  of  selected  hydrogen  isotopes  [Sou  86]. 


Triple  Temp 

Triple  Press 

Critical  Temp 

Critical  Press 

e-H2 

13.81   (K) 

7030(Pa) 
52.74 
(Torr) 

32.98  (K) 

1.29  (MPa) 
9677  (Torr) 

n-H2 

13.96 

7200 
54.01 

33.19 

1.32 
9903 

HD 

16.604 

12,370 
92.80 

35.91 

1.484 
11,133 

e-Dz 

18.69 

17,130 
128.51 

38.26 

1.650 
12,378 

n-D2 

18.73 

17,150 
128.66 

38.25 

1.665 
12,491 

Table  1-6.  Quantum  parameters  (A,  e,  a,)  and  bulk  triple  and  critical  point  temperatures 
of  selected  gases  [Liu  93]^ 


Gas 

A 

e/kb(K) 

a(A) 

T,(3D)  (K) 

Tc(3D)  (K) 

^He 

0.555 

10.22 

2.556 

- 

3.3102 

^He 

0.481 

10.22 

2.556 

- 

5.1899 

H2 

0.334 

37.0 

2.928 

13.804(e) 
13.956  (n) 

33.976(e) 
33.19  (n) 

HD 

0.273 

37.0 

2.928 

16.60 

35.91 

D2 

0.236 

37.0 

2.928 

18.69(e) 
18.73  (n) 

38.262  (e) 
38.34  (n) 

Ne 

0.108 

36.7 

2.788 

24.553 

44.40 

Ar 

0.0386 

119.80 

3.405 

83.806 

150.70 

Kr 

0.0221 

164.0 

3.624 

115.763 

209.5 

Xe 

0.0143 

230.4 

3.921 

161391 

289.72 

Table  data  is  from  [Liu  93]. 


CHAPTER  2 

ADSORPTION  ISOTHERMS 

Nam  et  ipsa  scientia  potestas  est. 
For  knowledge  itself  is  power. 
Religious  Meditations,  Of  Heresies,  Sir  Francis  Bacon  J  56  J  -  J  626 

2.1  Volumetric  Adsorption  Isotherms 

A  number  of  experimental  methods  have  been  used  in  recent  years  to  study  the 
adsorption  of  gases  onto  surfaces,  including  quartz  microbalance  [Kri  84,  Mis  92, 
Tab  92],  dielectric  measurements  [Lah  82],  NMR  [Cra  94],  third  sound  [Zim  92,  Spr  94], 
neutron  diffraction  [Wei  91b],  LEED  [Wliu  93],  and  ellipsometry  [You  90].  Computer 
models  have  also  been  built  and  simulations  run  [Phi  93].  However,  the  most 
straightforward  and  common  method  of  observation  is  the  volumetric  adsorption 
isotherm.  The  most  simple  version  of  this  method  utilizes  only  one  pressure  gauge  and  a 
number  of  valves  at  room  temperature.  An  amount  of  gas  is  loaded  into  a  known 
"dosing"  or  "loading"  volume,  also  at  room  temperature,  and  the  pressure  of  that  volume 
is  measured.  The  valve  between  the  loading  volume  and  the  sample  volume  is  then 
opened  and  the  gas  is  allowed  to  expand  into  the  sample  chamber.  The  final  pressure  of 
the  gas  above  the  substrate  is  then  measured.  From  these  two  measurements  and  the  ideal 
gas  law,  the  number  of  atoms  that  are  adsorbed  onto  the  surface  can  be  calculated.  The 
number  of  atoms  adsorbed  onto  the  substrate  at  a  fixed  temperature  is  then  plotted  against 
the  measured  final  pressure  in  a  graph  routinely  called  an  adsorption  isotherm  curve. 
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One  advantage  of  volumetric  isotherms  is  their  simpUcity.  The  method  requires 
only  the  repeated  measurement  of  the  gas  pressure.  For  a  fixed  temperature,  the  chemical 
potential  is  a  monotonic  function  of  the  equilibrium  vapor  pressure  of  the  gas  above  the 
surface.  Thus,  additional  properties  of  the  two-dimensional  system  can  be  calculated  fi"om 
one  or  more  of  the  adsorption  isotherm  curves.  For  example,  the  two-dimensional 
corapressibiUty  is  directly  related  to  the  slope  of  the  adsorption  isotherm  and  the  isosteric 
heat  of  adsorption  can  be  calculated  fi"om  two  isotherms.  It  is  usually  assumed  that  the 
chemical  potential  of  the  film  is  the  two-dimensional  chemical  potential  (interactions 
within  the  film)  plus  the  substrate  binding  energy  [Das  75].  The  quantity  adsorbed  onto 
the  surface  is  a  measure  of  the  average  film  density. 

More  than  one  isotherm  can  be  measured  at  once  thereby  maximizing  the  amount 
of  information  gained  per  experimental  run.  As  shown  in  Figure  2-1,  after  measuring  the 
first  final  pressure  (Pa)  and  calculating  the  associated  number  of  molecules  adsorbed  (Na), 
the  temperature  (Ta)  can  be  changed  to  a  new  (lower)  value  (Tb)  and  the  final  pressure 
measured  (Pb).  The  pressure  at  the  first  temperature  can  then  be  remeasured  (a  vs.  a')  to 
check  for  hysteresis  before  adding  more  gas  to  the  system  (a'  ->  c).  Note  that  point  (a) 
measures  adsorption  while  point  (a')  measures  desorption. 
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Figure  2-1.   Procedure  for  the  simultaneous  measurement  of  two  isotherms. 
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2.2  Sample  Preparation 

All  samples  of  boron  nitride  powder  used  in  this  thesis  come  from  the  single 
purchase  of  50  grams  of  boron  nitride  (99.5%  metals  basis)  from  JM  Electronics'  with  lot 
number  BI3B03.  The  boron  nitride  sample  was  pressed  through  a  50  micron  mesh  by  the 
manufacturer  and  was  received  in  the  form  of  a  fine  white  powder.  Boron  nitride  is  also 
sold  m  the  form  of  a  hot  pressed  rod,  but  this  form  was  not  obtained  for  these 
experiments. 

As  discussed  above,  since  the  hydrogen  sample  adsorbed  onto  the  surface  is 
relatively  small,  it  was  essential  that  additional  forms  of  hydrogen  be  minimized  on  the 
boron  nitride  surface.  Just  as  important,  the  molecular  hydrogen  itself  is  sensitive  to  any 
contaminants  with  large  magnetic  moments.  Adsorbates  with  large  magnetic  molecular 
moments  can  artificially  increase  the  rate  of  ortho-para  hydrogen  conversion.  Thus,  care 
was  taken  to  keep  the  sample  free  from  additional  water,  oils,  and  even  airborne 
contaminants.  All  equipment  was  cleaned  with  rubbing  alcohol  and  then  heated  in  a 
vacuum  to  remove  the  alcohol.  Tlie  boron  nitride  was  then  loaded  into  an  oven  and 
heated  to  300°  C  under  difilision  pump  vacuum  (approximately  10"  Torr)  overnight.  This 
procedure  is  similar  to  that  prescribed  by  Regnier,  Thomy,  and  Duval  [Reg  79]  in  their 
comparative  adsorption  isotherm  study  of  Xenon  and  Krypton  on  boron  nitride  and 
graphite. 

Samples  with  various  bake-out  temperatures  and  times  were  investigated  using 
methane  isotherms  to  gauge  the  sample  heterogeneity  and  cleanliness.  A  dirty  sample  was 
expected  to  have  either  additional  small  steps  due  to  the  adsorption  of  certain  surface 
contaminates  (e.g.,  water),  more  rounded  steps  due  to  inhomogeneity,  very  few  distinct 
steps,  or  a  shorter  first  step  due  to  occupation  of  adsorption  sites  by  contaminants. 
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However,  the  differences  in  the  isotherms  between  differently  prepared  samples  was 
negligible.  In  fact,  after  bake-out,  the  samples  were  surprisingly  stable.  Even  short 
exposiu-e  (hours)  to  air  did  not  seem  to  aUer  the  isotherms.  Such  stabihty  is  in  sharp 
contrast  to  some  other  substrates  (e.g.,  magnesium  oxide)  which  readily  adsorb  water 
from  the  air  and  must  be  both  manufactured  and  retained  in  vacuum  in  order  to  insure 
sample  cleanliness. 

An  attempt  was  made  to  increase  the  NMR  signal  by  maximizing  the  surface  area 
per  unit  volume.  Some  boron  nitride  samples  were  pressed  inside  a  cyHndrical  former 
during  bake-out.  The  result  was  a  white  rod  resembling  a  stick  of  chaUc  with  a  density 
approximately  5  times  that  of  the  powder  form.  Unfortunately,  upon  attempting  to 
measure  an  adsorption  isotherm  at  77  K  on  the  pressed  sample,  it  was  discovered  that  the 
equUibrium  tkne  increased  from  tens  of  minutes  to  tens  of  hours.  Tlie  sample  was  then 
broken  into  millimeter  size  pieces  and  the  isotherm  repeated.  However,  the  fragmentary 
sample  still  had  an  equilibrium  time  of  a  few  hours  at  77  K  which  was  deemed  to  be  too 
long  to  be  usefiil.  All  fiirther  experiments  were  performed  on  non-pressed  samples. 

A  recent  study  at  Southern  Illinois  University  by  Migone  et  al.  [Mig  94a]  has 
investigated  diffefent  sample  preparations  of  boron  nitride.  Samples  were  prepared  by 
baking  under  vacuum  at  temperatures  up  to  900  C  for  several  hours.  Some  samples  were 
also  washed  in  methanol.  Adsorption  isotherms  were  then  performed  utilizing  Argon  as 
the  adsorbate.  These  isotherms  were  then  analyzed  for  any  additional  features  and  the 
clarity  of  the  known  steps  and  substeps.  Additionally,  the  curvature  of  the  beginning  of 
the  first  step  was  monitored.  Added  curvature  along  this  section  of  the  isotherm,  shown 
in  the  Figure  2-2  (top),  would  indicate  the  presence  of  impurities  on  the  surface  providing 
adsorption  sites  with  an  adsorption  potential  stronger  than  that  of  rest  of  the  BN  surface 
[Mig  94b]. 

The  researchers  found  that  the  77.3  K  argon  isotherms  measured  on  the  dirty  BN 
samples  differed  from  those  measured  on  the  clean  samples  in  two  ways.  First,  an 
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Figure  2-2.  Adsorption  isotherm  of  argon  on  boron  nitride  taken  in  December  1993.  The 
first  and  part  of  the  second  layers  (top).  Detail  of  hquid-sohd  transition,  approximately 
1.5%  of  a  monolayer,  m  the  first  layer  (bottom).  The  hnes  serve  only  as  guides  to  the  eye. 
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Figure  2-3.  Detail  of  Argon  on  boron  nitride  isotherm  taken  in  June  1994.  Note  that 
neither  the  isotherm  (top)  nor  the  calculated  compressibihty  (bottom)  show  any  sign  of  a 
substep  at  six  Torr  as  identified  by  other  authors  for  a  dirty  sample. 
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additional  compressibility  peak  of  approximately  10  M/N  was  observed  near  6  Torr  with  a 
half- width  of  about  2  Torr.  This  low  pressure  peak  was  beUeved  to  be  due  to  the 
inadvertent  adsorption  of  oxygen  onto  the  surface  of  the  boron  nitride  yielding  boric 
oxide.  The  contaminant,  as  measured  by  the  compressibility  peak,  could  be  removed  by 
heat  pumping  at  900  C.  Second,  the  hquid-soUd  substep  centered  at  approximately  36 
Torr,  which  is  approximately  1.5  -  2%  of  a  monolayer  for  a  clean  sample,  diminished  m 
size  and  sharpness  for  a  dirty  sample  to  the  point  where  it  was  barely  measurable  from  the 
backgroimd  noise. 

Similar  isotherms  have  been  measured  at  77  K  on  the  boron  nitride  sample  used  m 
this  laboratory  for  hydrogen  studies  at  approximately  six  month  intervals  for  two  years. 
The  top  of  Figure  2-2  illustrates  an  isotherm  measured  in  December  1993  for  the  first  and 
part  of  the  second  layers  of  the  adsorbed  film  while  the  bottom  details  the  hquid-solid 
transition.  The  substep  is  approximately  1.5%  of  a  monolayer.  Figure  2-3  shows  a  low 
pressure  detail  from  a  June  1994  argon  isotherm  as  well  as  the  calculated  compressibihty 
of  the  adsorbed  film.  No  substep  is  evident  in  this  pressure  range  for  this  data  or  any 
earher  data  at  the  0.75%  of  a  monolayer  level.  Argon  measurements  were  not  performed 
on  untreated  or  dirty  samples.  It  is  possible  that  repeated  exposure  to,  or  flushing  with, 
hydrogen  removes  many  contaminants,  including  oxygen,  from  the  surface. 

Figure  2-4  presents  a  detail  of  an  argon  isotherm  measured  in  June,  1994.  A 
substep  near  38  Torr  is  quite  prevalent  in  this  data  as  well  as  all  earlier  and  later 
measurements.  Since  the  size  of  the  substep  remains  unchanged  in  time,  it  is  very  similar 
to  the  same  substep  measured  six  months  earlier,  the  sample  quality  does  not  seem  to 
deteriorate  under  experimental  conditions.  Tlie  compressibility  has  a  fiiU  width  at  half 
maximum  of  approximately  5  Torr  and  a  height  of  approximately  35  M/N.  Tliis  height  is 
very  similar  to  that  reported  elsewhere  [Mig  93]  though  somewhat  wider. 

Finally,  the  Southern  Illinois  group  reported  a  rather  large  (approximately  5%) 
hysteresis  between  adsorption  and  desorption  isotherms  for  argon  on  boron  nitride.   In 
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Figure  2-4.  DetaQ  of  Argon  on  boron  nitride  isotherm  measured  in  June,  1994.  The 
substep  at  38  Torr  (top)  signaling  the  Uquid-solid  transition  and  the  associated 
compressibihty  peak  (bottom)  are  clearly  evident  in  the  data. 
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other  words,  the  argon  adsorbed  at  a  higher  pressure  than  it  desorbed  in  that  region  of  the 
isotherm.  Similar  hysteresis  has  been  documented  for  other  systems.  Lysek  et.  al. 
[Lys91]  have  used  heat  capacity  techniques  to  verify  their  own  adsorption  isotherm 
studies.  They  found  the  surprismg  resuh  that  capillary  condensation  m  graphite  foam  can 
begin  at  1. 1  layers,  much  sooner  than  previously  thought.  It  should  be  noted  that  the  BN 
powder  used  m  this  study  is  much  less  dense  than  any  of  the  typical  forms  of  graphite 
currently  used  as  substrates. 

A  similar  adsorption-desorption  cycle  was  measured  for  approximately  ten  layers 
of  argon  on  boron  nitride  and  is  presented  in  Figure  2-5  (top).  A  pressure  dependent 
hysteresis  was  found  to  be  centered  around  layers  three  and  four.  A  measure  of  the 
systematic  errors  associated  with  the  adsorption  isotherm  experiments  may  be  made  from 
this  cycle.  As  shown  m  Figure  2-5  (bottom),  the  hquid-solid  substep  mentioned  above 
occurs  at  the  same  pressure  along  both  the  adsorption  and  desorption  curves.  Tlius  the 
displacement  shown  in  the  figure  is  a  measure  of  the  error  in  the  number  of  atoms 
adsorbed  onto  the  surface.  The  two  curves  are  displaced  approximately  0.05  cc-STP  after 
a  total  adsorption  of  18  cc-STP  m  100  steps.  Thus,  an  order  of  magnitude  estimate  for 
the  systematic  error  in  the  number  of  atoms  adsorbed  is  0.003  %  per  step. 

The  measurement  error  for  the  argon  isotherms  was  greater  than  that  for  hydrogen 
for  two  reasons.  First,  the  tune  to  vapor  pressure  equilibrium  was  found  to  be  greater  for 
argon  than  for  hydrogen.  It  is  believed  that  this  time  increase  was  a  direct  result  of  the 
more  massive  argon  traveling  at  a  slower  speed  than  the  extremely  light  hydrogen  which 
required  additional  time  to  pressure  equiUbrate  across  the  capillary.  Second,  the 
apparatus,  including  the  temperature  controller,  was  designed  for  low  temperature  work. 
Thus  a  liquid  nitrogen  bath  was  used  to  achieve  and  maintain  the  higher  temperature. 
However,  the  vapor  pressure  of  the  argon  above  the  substrate  was  found  to  depend  upon 
the  heiglit  of  the  bath  (and  hence  the  temperature)  and  this  height  was  not  stable  over  the 
week  of  data  taking.  Tims  small  fluctuations  occurred  in  the  data  as  the  bath  was  refilled. 
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Figure  2-5  Hysteresis  loop  in  the  adsorption-desorption  cycle  (top).  The  data  was 
recorded  at  two  different  times  along  the  same  curve  as  a  measure  of  residual  pressure 
drift.  The  hquid-soUd  substep  may  be  used  as  a  measure  of  systematic  errors  (bottom). 
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As  a  practical  matter,  it  is  difficult  to  distinguish  between  a  hysteresis  loop  with 
this  particular  pressure  dependence  and  a  systematic  measurement  of  the  final  vapor 
pressure  before  equUibrium  has  been  fully  estabhshed.  The  computer  algorithm  used  to 
gather  the  data  for  this  thesis  contained  a  built-in  redundancy  (shown  above  in  Figure  2-1). 
Specifically,  even  for  the  measurement  of  a  single  isotherm  (a  one  temperature  isotherm), 
the  final  pressure  was  recorded  twice:  once  after  all  criterion  were  met  and  again  at  some 
time  later  (approximately  5  minutes).  The  isotherm  could  then  be  analyzed  with  either  or 
both  final  pressure  values  and  any  differences  easily  compared.  This  method  was  found  to 
be  quite  usefiil  in  determining  the  equilibrium  criterion  for  the  different  pressure  ranges 
during  the  development  of  the  computer  program. 

In  Figure  2-5,  both  sets  of  data  are  presented  with  the  first,  or  early,  data 
presented  as  a  small  circle  and  the  second,  or  later,  data  shown  as  an  open  triangle.  Note 
that  the  two  symbols  are  displaced  more  for  the  desorption  curve  than  for  the  adsorption 
curve.  While  this  may  have  physical  significance,  the  nature  of  the  experimental  apparatus 
also  contributes  to  the  effect.  Namely,  the  pressure  gauge  is  located  at  room  temperature 
in  the  caUbrated  volume  which  is  separated  fi'om  the  sample  cell  by  a  long  capillary.  Thus, 
the  gauge  is  more  sensitive  to  the  adsorption  than  the  desorption  process.  Some  of  the 
hysteresis  shown  in  Figure  2-5  (top)  may  be  due  to  a  systematic  measurement  error,  but 
this  effect  is  not  large  enougli  to  account  for  the  entire  loop. 

2.3  Adsorption  Isotherm  Apparatus  and  Procedures 

Tlie  adsorption  isotherm  experiment  was  developed  to  characterize  the  boron 
nitride  for  use  in  the  low  temperature  hydrogen  NMR  ex-periment.  Tlie  isotherm 
apparatus  started  out  as  a  very  simple  set  of  valves  connected  with  common  copper  pipe 
mounted  on  a  plywood  sheet  and  a  brass  sample  cell  inserted  into  a  liquid  nitrogen  bath. 
This  setup  was  usefiil  for  methane  and  argon  isotherms  and  yielded  some  published  data. 


However,  the  isosteric  heat  of  adsorption  could  not  be  measured  very  precisely  as  only 
single  temperature  isotherms  could  be  performed.  More  importantly,  taking  data  by  hand 
is  inherently  slow  relative  to  a  computer  operated  design  that  can  take  data  day  and  night. 

The  first,  second,  and  third  generations  of  hquid  heUum  temperature  isotherm 
sticks  were  shown  to  suffer  fi-om  small  temperature  gradients  across  the  sample  cell.  The 
first  stick  had  a  100  mil  brass  sample  cell  wall.  The  second  generation  increased  this  to 
1/4  inch.  The  third  generation  replaced  the  brass  (chosen  for  ease  of  machining)  with 
copper  which  has  a  larger  thermal  conductivity.  A  temperature  difference  of  100  mK  at 
14  K  was  found  to  be  too  large  for  isothermal  measurements.  Such  a  temperature 
gradient  causes  one  part  of  the  sample  to  be  undergoing  a  transition  before  another 
portion  therefore  washing  out  the  finer  details  of  the  data. 

The  final  stick  utilized  a  two-wall  system  analogous  to  heat  capacity  measurement 
systems.  The  1/4  inch  copper  outer  wall  serves  as  an  adiabatic  shield  for  the  thin- walled, 
50  mil,  copper  sample  cell.  The  space  between  the  adiabatic  wall  and  the  sample  cell  was 
then  loaded  with  Uquid  or  solid  hydrogen  to  be  used  in  the  experiment.  The  low 
temperature  hydrogen  storage  had  three  purposes.  First,  the  hydrogen  vapor  minimized 
the  temperature  gradients  across  tlie  sample  cell  even  if  the  adiabatic  wall  had  a  small 
temperature  gradient.  Second,  the  low  temperature  storage  acted  as  an  elemental  filter. 
The  only  gases  with  substantial  vapor  pressures  at  hydrogen  temperatures  are  the 
hydrogen  isotopes  themselves  (i.e.,  H2,  HD,  D2),  helium-3,  hehum-4,  and  neon.  Thus,  any 
other  gas  commg  fi-om  either  the  hydrogen  supply  bottle  or  the  supply  plumbing  would 
freeze  out  in  the  low  temperature  hydrogen  supply  cell.  Third,  while  hydrogen  above 
100  K  consists  of  25  percent  para-hydrogen  (spin  0)  and  75  percent  ortho-hydrogen, 
(spin  1)  at  lower  temperatures  the  ratio  of  these  two  species  depends  of  the  temperature. 
Noting  that  the  vapor  pressure  difference  between  ortho-  and  para-hydrogen  is 
measurable,  it  is  obvious  that  the  use  of  room  temperature  hydrogen  creates  an  additional 
time  dependent  pressure  change  on  a  time  scale  similar  to  that  of  the  experiment. 
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Figure  2-6.    Schematic  of  gas-handling  system  showing  the  automated  valves.    Additional 
hand-operated  valves  (not  shown)  were  used  ordy  for  safety  and  pumping  out  the  sample. 


2.3.1  Computer  Automated  Valve  System  and  Adsorption  Algorithm 

The  adsorption  isotherm  apparatus  consisted  of  five  computer  operated  valves 
with  manual  override,  as  well  as  various  hand-operated  valves.  The  valves  were  standard 
Nupro^  "H"  series  bellows  valves.  Two  pressure  gauges  were  used  at  room  temperature. 
The  first,  by  MKS^,  had  a  fiill- scale  reading  of  10  Torr  while  the  second,  by  Druck^,  had  a 
fiill  scale  readmg  of  1000  Torr. 

For  simpHcity,  as  well  as  economy,  both  pressure  gauges  were  located  directly  on 
the  loadmg  volume  and  served  the  double  purpose  of  measuring  the  loading  and  final 
pressures.  This  was  made  possible  by  reducing  the  loading  volume  which  includes  the 
pressure  gauge  volume  to  a  minimum,  about  3.3  cubic  centimeters.  Tliese  pressures  were 
recorded  by  the  computer  and  presented  on  a  strip  chart  on  the  computer  screen. 


^  Nupro  Company,  4800  East  345th  Street,  Willoughby,  Ohio  44094. 
^  MRS  Instruments,  Inc.,  24  Walpole  Park  South,  Walpole,  MA,  02081. 


Druck  Limited,  Groby,  Leicester,  England. 
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After  a  waiting  period  of  approximately  one  minute,  the  computer  began  to  sample 
the  pressure  above  the  substrate  about  once  every  5  seconds.  The  program  then 
calculated  both  a  first,  as  well  as  a  second  order  polynomial  fit  to  the  previous  20 
measured  pressures.  For  the  acceptance  of  a  pressure  reading,  all  orders  of  the  polynomial 
fits  had  to  meet  certain  criterion.  EquUibrium  was  determined  when  all  of  the  criteria  were 
met  at  least  one  dozen  times  (more  for  lower  pressures). 

The  criteria  utilized  to  measure  equihbrium  were  a  fiinction  of  the  pressure. 
Lower  pressure  data  required  longer  equihbrium  tunes  and  thus  had  tighter  criteria. 
Equilibrium  times  varied  between  2  hours  for  equUibrium  pressures  below  0. 1  Torr  to 
about  5  minutes  for  pressures  above  100  Torr.  Argon  required  longer  equilibrium  tunes 
than  the  hydrogen  isotopes  for  a  given  final  pressure.  For  simuhaneous  isotherm 
experiments,  reproducibility  was  measured  for  every  dose  by  returning  to  the  first 
temperature  and  remeasuring  the  final  pressure.  Hysteresis  was  never  observed  in  the 
data. 

2.3.2  Temperature  Control 

The  ability  to  measure,  under  the  right  conditions,  the  layer  by  layer  growth  of  the 
adatoms  on  the  surface  through  a  pressure  measurement  is,  due  to  the  thermal  properties 
of  the  gas,  hmited  by  the  control  of  the  temperature.  For  precision  isotherms,  the  strong 
temperature  dependence  of  the  vapor  pressure  of  a  gas  must  be  taken  into  account. 

Tlie  temperature  controller  used  for  the  adsorption  isotherm  experiments  was 
modified  from  an  in-house  designed  and  built  seven  channel  temperature  control  unit. 
Figure  2-7  presents  the  schematic  diagram  of  the  main  elements  of  the  unit.  A  silicon 
diode  served  as  the  thermometer  for  the  controller.  The  band  gap  of  the  diode  was 
approximately  0.87  Volts  at  room  temperature  and  4.5  Volts  at  4  K.  The  band  gap  in  a 
semiconductor  is  a  strong  (approximately  exponential)  ftinction  of  temperature  and  is  due 
to  the  reduction  m  the  number  of  carrier  electrons  and  holes  in  the  conduction  bands.  Tlie 
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Figure  2-7.  Electrical  schematic  of  tenqperature  controller. 
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particular  diode  used  for  the  thermometer  was  chosen  from  a  set  of  diodes  after  extensive 
tests  and  measurements  of  their  thermal  properties. 

The  sensitivity  of  the  sihcon  thermometer  was  in  large  part  determined  by  the  noise 
of  the  constant  current  source  and  the  operatmg  temperatures.  The  stabihty  achieved  at 
low  temperatures  was  approximately  1  part  in  10^.  At  higher  temperatures,  the  stabihty 
was  decreased  by  approximately  an  order  of  magnitude.  The  system  was  not  useftil  for 
temperatures  approximately  greater  than  30  K. 

The  temperature  control  problem  for  the  adsorption  isotherms  were  comphcated 
by  the  length  of  the  experiments,  namely  1-3  weeks.  Over  this  long  time  period  the 
cooling  power  from  the  cold  bath  to  the  sample  cell  decreased  greatly  as  the  liquid  helium 
reservoir  evaporated.  Since  temperature  control  is  by  its  very  nature  a  balancing  act 
between  all  of  the  warming  and  cooling  power  inputs,  such  a  strong  drift  in  the  power  of 
one  of  the  components  had  a  dramatic  effect  upon  the  long-time  stability  of  the 
temperature.  Active  electronic  temperature  control  with  the  usual  feedback  systems  was 
found  to  be  inadequate  for  the  precision  required  for  the  adsorption  isotherm  ex^jeriments 
and  other  means  of  control  were  searched  out. 

A  solution  was  found  in  the  use  of  a  two  wall  dipstick  to  form  a  passive 
temperature  controller.  Tlie  outer  wall  of  the  dipstick  was  put  in  contact  with  the  liquid 
helium,  the  "cold  bath,"  in  the  transfer  dewar.  Tlie  inner  wall  was  themially  isolated  from 
the  outer  wall  with  one  important  exception.  A  residual  amount  of  helium  gas  was  placed 
between  the  two  walls  of  the  dipstick  to  act  as  a  thermal  exchange  gas.  Tliis  is  a  common 
procedure  in  low  temperature  physics.  In  fact,  only  hehum  can  exist  as  a  gas  in  this  space 
since  all  other  gases  would  "freeze  out"  when  put  in  contact  with  the  4  K  bath. 

Tlie  system  was  made  into  a  passive  temperature  controller  by  the  choice  of 
hydrogen  as  the  exchange  gas  to  provide  the  thermal  conduction  between  the  inner  wall  of 
the  dipstick  and  the  outer  wall  of  the  supply  cell  in  which  the  sample  cell  resided.  As  the 
level  of  the  hquid  helium  in  the  transfer  dewar  lowered  with  time,  the  thermal  conductance 
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through  the  heUum  exchange  gas  diminished.  This  eflfect  caused  the  pressure  of  the 
hydrogen  exchange  gas  to  increase  which,  in  turn,  led  to  an  increase  in  the  thermal 
conductance  of  the  hydrogen!  In  short,  the  decreasing  cooling  power  of  the  cold  bath  as  a 
flmction  of  hquid  hehum  height  was  offset  by  the  increasing  cooling  power  of  the 
exchange  gas. 

The  use  of  a  two  wall  dipstick  m  a  dewar  and  the  specific  choice  of  exchange  gases 
acting  to  form  a  passive  temperature  controller  has  uses  beyond  the  hehum/hydrogen 
system  described  above.  If  one  were  studying  an  argon  system  in  a  liquid  nitrogen  dewar, 
for  example,  then  one  may  choose  nitrogen  and  argon  as  the  thermal  exchange  gases  to 
provide  a  similar  passive  temperature  controller  at  a  very  different  temperature  range. 

Hydrogen  exchange  gas  was  chosen  over  helium  for  this  inner  space  for  two  more 
reasons.  Tlie  first  reason  concerned  limiting  the  cooling  power  of  the  cold  bath  while  the 
second  reason  involved  maintaining  the  sample  cell  at  a  temperature  less  than  any  portion 
of  the  fill  capillary.  Tlie  vapor  pressure  of  bulk,  sohd  hydrogen  is  a  strong  fiinction  of 
temperature  which  is  in  the  millitorr  range  at  approximately  10  K.  Thus,  the  use  of 
hydrogen  as  an  exchange  gas  provides  a  natural  thermal  conductance  limit  between  tlie 
experiment  and  the  cold  bath.  Tlie  same  thennal  properties  of  hydrogen  conspire  to  make 
hydrogen,  when  used  as  an  exchange  gas,  provide  a  natural  low  temperature  limit  on  tlie 
sample  cell,  approximately  10  K.  Tlie  thermal  conductance  between  the  two  walls  and 
between  the  inner  wall  and  the  ex-periment  was  set  and  controlled  by  keeping  the  pressure 
of  both  exchange  gases  in  the  millitorr  range.  In  practice,  it  was  found  that  a  helium 
exchange  gas  pressure  of  approximately  60  millitorr  and  a  hydrogen  exchange  gas 
pressure  of  1-200  millitorr  worked  well. 

Tliemial  uniformity  across  the  sample  was  achieved  in  two  ways.  First,  the  sample 
cell  was  placed  inside  the  supply  can.  Second,  the  walls  of  the  supply  can  were  made 
relatively  thick,  0.1  inches,  and  of  oxygen  fi^ee  high  conductivity  (OFHC)  copper.  As  side 
benefits,  the  hydrogen  supply  used  in  the  adsorption  experiments  both  achieved  ortho-para 
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equilibrium  and  most  impurities  in  the  sample  gas  were  filtered  out.  Finally,  the  supply 
can,  the  hydrogen  supply  gas,  and  the  sample  cell  provided  enough  thermal  mass  to  act  as 
a  filter  blocking  out  the  high  frequency  thermal  changes. 


2.4  Hydrogen  Isotope  Isotherms 

The  adsorption  of  hydrogen  isotopes  onto  graphite  has  been  extensively  studied  by 
H.  Freknuth  and  H.  Wiechert  et  al.  [Fre  90,  Wie  91],  J.M.  GottUeb  and  L.W.  Bruch 
[Got  90],  and  others  [Cha  90],  and  two-dimensional  phase  diagrams  have  been 
constructed  (see  Figure  2-9).  All  three  common  hydrogen  isotopes,  H2,  HD,  and  D2,  were 
used  in  this  study  to  allow  for  isotopic  comparison.  Two  goals  existed  for  the  isotherms: 
One,  to  map  out  the  translational  phase  diagrams  (e.g.,  the  critical  and  triple  points  and 
translational  phase  transitions);  two,  to  measure  the  isosteric  heat  of  adsorption. 

2.4.1   Interpretation 

The  adsorbate  gas  is  at  first  adsorbed  onto  the  surface  directly  into  a  solid  state 
with  a  very  low  residual  pressure.  Build-up  of  the  first  layer  continues  with  little  change  in 
the  vapor  pressure  and  is  marked  by  the  nearly  vertical  climb  in  an  isothenn  plot.  Tlie 
monolayer  is  complete  near  a  point  historically  marked  as  "B"  or  "Bl"  which  occurs  soon 
after  the  first  vertical  rise.  Tlie  following,  more  horizontal,  nearly  linear  region  e.xlends  to 
the  beginning  of  the  second  layer  where  the  adatoms  usually  adsorb  into  a  two- 
dimensional  gaseous  state.  Adsorption  into  the  gaseous  state  is  manifested  by  a  decrease 
in  the  isosteric  heat  of  adsorption  (q^t)  to  a  value  close  to  the  three-dimensional  heat  of 
condensation  (0.91  IkJ/mole  or  109.6  K  for  H2).  Tlie  second  and  higher  layers  are  marked 
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Figure  2-8.  Temperature  control  as  a  function  of  time  for  two  experimental  runs, 
BNHD.3T  (top)  and  BND2.28  (bottom). 
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Figure  2-9.  Translational  phase  diagrams  of  the  second  layer  of  H2  on  MgO  [Ma  89]  (top) 
and  graphite  [Sch  89]  (bottom). 
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by  similar  changes  both  in  the  slope  of  the  isotherm  and  the  magnitude  of  the  isosteric  heat 
of  adsorption.  After  a  number  of  layers,  these  slope  changes  fade  out  and  the  isosteric 
heat  of  adsorption  approaches  the  three-dimensional  heat  of  condensation. 

2.4.2  HP  Isotherms 

The  adsorption  of  deuterium  hydride  (HD)  was  chosen  over  H2  or  D2  for  detailed 
investigation  due  to  the  possibihty  of  comparing  these  results  with  the  interesting  results 
of  two  earUer  studies,  HD  on  graphite  [Liu  93]  and  HD  on  MgO  [Ma  89].  Both  systems 
were  not  only  investigated  via  adsorption  isotherms  but  also  with  heat  capacity  techniques. 
Both  systems  were  also  further  explored  via  quasielastic  neutron  scattering  (QENS),  and 
phase  diagrams  were  reported  [Wei  91,  Vil  94].  Both  studies  found  substeps  near  the 
completion  of  the  second  layer  indicative  of  the  two-dimensional  hquid-solid  transition. 
However,  the  size  of  the  substep  varied  greatly.  For  HD  on  MgO  the  substep  was 
approximately  11%  of  a  monolayer  wlrile  the  same  step  was  found  to  be  only  2%  of  a 
layer  for  HD  adsorbed  onto  graphite.  The  change  in  entropy  associated  with  the  making 
of  HD  on  graphite  corresponds  to  approximately  8.5%  of  the  bulk  melting  entropy  at  its 
triple  point  [Vil  94]. 

Figure  2-10  presents  the  resuUs  of  a  muUiple-temperature  (simultaneous) 
adsorption  isotherm  run  of  HD  on  BN  starting  with  the  end  of  the  first  layer.  Tiie  data 
fi-om  eleven  temperatures  are  plotted  and  five  steps  are  resolvable  in  the  lowest 
temperature  data.  Note  that  the  pressures  at  the  higliest  temperature  were  remeasured  as 
described  above  and  no  hysteresis  encountered.  Tlie  vapor  pressure  above  the  surface 
during  the  formation  of  the  first  layer  is  too  low  to  be  measured  reproducibly  with  our 
apparatus. 

To  check  for  possible  substeps,  another  set  of  isotherms  was  collected  using  finer 
steps  and  more  temperatures  ( 1 7).   Over  2500  data  points  were  measured  over  a  period  of 
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Figiire2-10.  Volumetric  vapor  pressure  isotherms  of  HD  on  BN.  Five  steps  are  visible  in 
the  low  temperature  data.  The  data  from  all  temperatures  was  recorded  simultaneously; 
lines  are  guides  to  the  eye. 
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Figure  2-11.  Details  of  two  layers  in  the  HD  isotherms.  Note  the  lack  of  a  substep  m  the 
middle  of  the  second  layer  (bottom)  and  a  possible  substep  near  the  completion  of  the 
third  layer  (top).  Lines  are  guides  to  the  eye. 
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22  days.  As  shown  in  Figure  2- 11  (bottom),  no  substep  was  found  at  the  1%  of  a 
monolayer  level  for  the  middle  of  the  second  layer  of  HD  adsorbed  onto  boron  nitride. 
Some  data  points  (near  2.72  cc-STP  and  elsewhere)  were  omitted  from  this  figure  due  to 
temperature  controller  error. 

The  top  portion  of  Figure  2-10  focuses  on  a  possible  substep  near  the  end  of  the 
third  layer,  between  4.85  and  4.95  cc-STP.  The  slope  change  in  this  region  can  not  be 
attributed  to  temperature  controller  error  since  the  pressure  changes  associated  with  the 
temperature  fluctuations  measured  by  the  thermometer  for  this  region  are  an  order  of 
magnitude  smaller  than  the  observed  pressure  deviations.  More  importantly,  the  change  in 
slope  is  observed  over  three  cycles  of  data  points.  Thus,  the  slope  change  can  not  be 
attributed  to  one  random  event. 

Two  more  items  of  interest  are  manifest  in  the  bottom  portion  of  Figure  2-11. 
First,  an  example  of  a  single  event  error  in  the  data  can  be  seen  by  the  deviation  of  the  line 
comprised  of  filled  circles  and  measured  at  an  average  temperature  of  11.61  K.  Second, 
the  slope  of  the  rise  portion  of  the  isotherm  is  sUglitly  negative  for  some  isotherms.  Wliile 
an  isotherm  with  a  vertical  rise  (an  infinite  slope)  has  a  straiglitforward  physical 
interpretation  (it  corresponds  to  a  path  through  a  two  phase  region)  a  negatively  sloped 
isotherm  is,  in  fact,  unphysical  Note  that  such  negative  slopes  are  found  in  the  data  only 
for  the  lowest  temperatures  and  pressures  measured  and  that  the  noise  in  the  data  is  due  to 
the  resolution  of  the  low  pressure  (10  Torr)  pressure  gauge,  0.001  Torr.  It  is  believed 
that  the  thermomolecular  effect  (a  difference  in  the  equilibrium  pressure  of  a  gas  across  a 
capillary  with  the  ends  at  different  pressures)  plays  some  role  here.  Thus,  the  means  of 
achieving  such  low  temperatures,  a  decrease  in  the  fill  line  heater,  is  a  possible  culprit  as 
the  lack  of  direct  heat  may  allow  some  portion  of  the  fill  line  to  become  cooler  than  the 
sample  cell.  Tlie  final  pressure  (measured  at  room  temperature)  may  be  the  pressure 
corresponding  to  the  thermomolecular  temperature  (for  the  capillary  diameter  and  the 
particular  gas)  with  the  pressure  of  the  sample  cell  as  a  lower  bound.    As  the  number  of 
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atoms  admitted  to  the  sample  cell  increases  (i.e.,  traveling  along  the  isotherm  in  the 
adsorption  direction),  the  pressure  in  both  the  sample  ceU  and  the  capillary  increases 
causing  the  thermal  conductance  of  the  system  along  the  capillary  to  increase  and  the 
temperature  gradients  to  return  to  normal  (linear  from  hot  to  cod  end). 

Much  eflFort  was  expended  in  attempting  to  eliminate  cold  spots  in  the  fill  line.  For 
example,  the  line  heater  was  made  more  uniform  and  a  shield  added  to  the  fill  line.  In 
addition,  the  direction  of  the  cycle  of  pressure  measurements  was  reversed.  The  earhest 
attempts  at  measuring  muhi-temperature  isotherms  first  measured  the  final  pressure 
corresponding  to  the  lowest  temperature,  jumped  to  measure  the  pressure  corresponding 
to  the  highest  temperature,  then  stepped  back  to  the  lowest  temperature.  After 
experimenting  with  all  four  possible  temperature  cycles  it  was  found  that  a  cycle  which 
measured  the  pressure  at  the  highest  temperature,  jumped  to  the  lowest,  then  stepped  back 
to  the  highest  produced  the  most  satisfactory  data.  This  change  in  the  temperature  cycle 
reduced  the  magnitude  of  the  negative  slope  approximately  5-7  times  compared  to  the 
original  procedure.  However,  this  effect  marks  the  lower  limit  of  the  experiment  at 
approximately  0.01  Torr. 

However,  the  best  solution  involved  redesigning  the  fill  line  to  be  inside  an 
evacuated  tube.  Theoretically,  such  an  arrangement  allows  for  the  most  linear 
temperature  gradient  along  a  tube  at  two  temperatures  of  any  passive  method.  Active 
methods,  for  example,  a  series  of  thermometers  and  heaters,  were  not  attempted.  A  new 
stick  was  built  but  not  yet  used  at  the  time  of  this  writing  because  it  was  determined  that 
the  errors  involved  in  calibrations  between  two  sticks  were  greater  than  the  loss  of  this 
non-essential  information.  Thus  the  plan  was  laid  to  finish  the  current  experiments  before 
the  introduction  of  the  new  stick.  Tlie  improved  stick  would  also  allow  for  concurrent 
beat  capacity  ex-periments  which  was  not  possible  with  the  old  stick. 

The  top  of  the  second  layer  was  not  investigated  as  thorouglily  as  the  same  area  in 
the  third  layer.   Further  investigations  are  currently  underway. 
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2.4.3  H?  and  D;  Isotherms 

Figures  2-12,  2-13,  and  2-14  exhibit  the  resuUs  of  isotherms  from  hydrogen  and 
deuterium  on  boron  nitride.  For  both  systems  at  least  four  steps  are  present  and  the 
critical  temperatures  have  been  calculated  (see  Table  2-1).  Neither  system  showed  any 
sign  of  Uquid-soUd  subsets  at  the  2%  level  nor  any  hysteric  effects. 


2.5  Compressibility,  Critical  Temperatures,  and  Larher's  Method 

At  temperatures  near  the  critical  point,  the  formation  of  each  layer  on  a  substrate 
surface  involves  the  two-dimensional  condensation  of  each  layer  from  a  dilute  to  a  dense 
adlayer.  Such  a  transition  is  of  first  order  and  must  end  at  a  critical  point.  Below  the 
critical  temperature  (e.g.,  in  the  hquid-vapor  coexistence  region),  the  chemical  potential 
and,  therefore,  the  vapor  pressure  is  constant  for  a  frxed  temperature.  Thus,  the  isotherm 
is  parallel  to  the  number  axis  in  this  region.  Since  the  compressibility  of  a  film  is 
proportional  to  the  slope  of  the  isotherm,  the  inverse  compressibility  vanishes.  However, 
for  real  systems,  surface  inhomogeneity  causes  isotherms  in  this  region  to  have  a  finite 
slope  [Eck  86].  Tlie  inverse  compressibihty  for  a  film  at  fixed  temperature  has  been  given 
by  [Das  75]  as 

/(_    =  n   kf,  r— —  I  Equation  2- 1 

where  n  is  the  areal  film  density.  Tlie  decrease  in  the  inverse  compressibility  was  first  used 
by  Y.  Larher  [Lar  79]  in  his  study  of  argon  on  CdCU  to  determine  the  critical  point  of  the 
film.  He  found  a  nearly  linear  rise  in  the  maximum  inverse  compressibility  as  a  fiuiction  of 
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Figure  2- 12.  Volumetric  vapor  pressure  isotherms  of  H2  on  BN.  Five  steps  are  visible  in 
the  low  temperature  data  with  the  last  four  shown.  The  data  from  all  temperatures  was 
recorded  simultaneously;  lines  are  guides  to  the  eye. 
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PRESSURE  (Torr) 

Figure  2-13.    Low  precision  volumetric  vapor  pressure  isotherms  of  D2  on  BN  detailing 
the  third  step.  Data  is  from  two  runs;  lines  are  guides  to  the  eye. 
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Figure  2-14.  High  precision,  volumetric  adsorption  isotherms  of  Dj  on  BN  detaihng  the 
third  step.  Data  is  from  one  multi-temp eratiu-e  isotherm  run.  The  gap  in  the  data  at 
approximately  3.5  cc-STP  is  due  to  experimental  difficulties  and  has  no  physical 
significance.  The  Usted  errors  represent  the  standard  deviation  of  the  temperature 
associated  with  each  volume/pressure  point. 
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temperature  above  the  critical  point  and  a  near  zero  inverse  compressibility  below  it.  The 
intersection  of  these  two  lines  was  determined  to  be  the  critical  temperature.  Later  work 
has  shown  that  Larher's  method  of  critical  point  determination  overestimates  the  critical 
temperature  by  approximately  4%  [Eck  86].  Comparisons  between  adsorption  isotherm 
and  heat  capacity  resuhs  bear  this  difference  out.  The  authors  flirther  pomted  out  that  the 
actual  overestimation  for  a  real  system  depends  on  the  temperature  range  of  the  data  and 
the  equation  of  state  of  the  adsorbed  system. 

Figure  2- 1 5  presents  a  Larher  plot  for  the  second  layer  of  HD  on  BN.  Two  lines 
are  drawn  and  for  clarity  only  some  of  the  error  bars  are  e?diibited.  The  horizontal  Une  is 
drawn  to  represent  the  positive  finite  slope  expected  of  real  systems.  In  this  region,  the 
data  is  liiglily  scattered  due  to  apparatus  diiBculties  in  measuring  low  vapor  pressures. 
Tlie  positive  sloped  hne  corresponds  to  the  linear  region  beyond  the  critical  temperature. 
At  higher  temperatures,  the  linear  approximation  breaks  down.  From  this  plot  the  critical 
temperature  is  estimated  to  be  at  12.3  ±  0.5  K.  This  temperature  is  similar  to  but  slightly 
higher  than  the  value  for  HD  on  graphite  of  1 1.8  +  0.3  K  [Liu  93]. 

Similar  plots  for  the  third  and  fourth  layer  of  HD  on  BN,  top  and  bottom  of 
Figure  2-16,  respectively,  yield  critical  temperatures  of  12.2  ±  0.5  K  and  11.8±0.5K, 
respectively.  Tlie  critical  temperatures  of  H:  and  D2  on  boron  nitride  have  been  similarly 
calculated  and  representative  plots  are  presented  in  Figure  2-17.  Tlie  critical  temperatures 
are  also  reported  in  the  tables  at  the  end  of  this  section. 

It  is  important  to  remember  that  graphite  has  been  the  subject  of  many 
investigations  using  a  number  of  methods  over  the  last  few  decades.  Further,  some  studies 
have  used  multiple  methods  to  determine  the  critical  points.  For  example,  the  HD  on 
graphite  studies  of  [Liu  93]  used  both  heat  capacity  measurements  and  adsorption 
isotherm  techniques  to  detennine  the  critical  temperatures.  Tlie  data  presented  herein 
represents  the  first  study  of  hydrogen  on  BN  and  only  adsorption  isothenii  data  is 
available.  Therefore,  the  critical  temperatures  hsted  here  should  be  viewed  as  preliminary 
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Figure  2-15.  Larher  plot  of  the  second  layer  of  HD  on  BN.  The  symbols  m  the  legend 
denote  different  experimental  runs.  The  numeric  code  represents  the  month  and  day  of 
1993  on  which  the  experimental  run  was  begun. 
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Figure  2-16.  Larher  plot  of  the  third  (top)  and  fourth  (bottom)  layers  of  HD  on  BN. 


14- 


64 


Osi- 


ers - 


C3 


1]    '    I    I    '    f~i    I    I"!    '    I    I    I    I    I    I    n — I — I — I — I — 1 — I — I — I — r — I — I — I — I — I — I — I — I — 1 — I — I — I — I — I — I — 1 — I — I — 1 — r 


/\^^^^  BND2.28 
<itt«3»ia  BND2.2ia 


/ 


y 


i^  _L 


• 


Layer    3 
D2     /    BN 


Tc     =     13.9     H —     0.5     K 


'  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

13  14-  15  16  17 

Tem.pereLtu.re     (K) 


CD  _ 


CD 


v-^       CD 


P^ 


OJ_ 
CD  _ 


CD 


CD 

CD- 


I    I    I    1 — I — I — I — I    I    I — I — I    I    I — I    I    I    I    I — I    I    I    I    I — I    I    I — n — I    I    I    I    I — 1    I    i — r— I — I — r—i — n — i    i    i — i — r 


nr****  bn.li.2.116 
/s^^A/s  bnl:i2.21 
ooooo  bxib.2.722 


.>y.r 


_         Layer     2 
H2     /     BN 

Tc     =      11.5     H-—      l.O 


CD — I   I    I   I   I    I    I    I    I   I   I   I   I   I  1   I   r   I   I   I   1   ''   I   I    'I   I    I   I    I   '"'''   I   '   'I    I   I   '    I    I    I   '   I    '   I 
MO  11  12  13  14-  15 


Temper  a  tiJire     (K) 


Figure  2-17.  Larher  plot  of  the  third  layer  of  D2  (top)  and  the  second  layer  of  H2  (bottom) 
onBN. 
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until  other  methods  verify  these  results.     Finally,  note  that  Larher's  method  tends  to 
overestimate  the  critical  temperatures,  as  mentioned  above. 

Adsorption  isotherms  have  been  used  to  identify  the  two-dimensional  melting  line 
and  triple  point.  A  substep  in  the  adsorption  isotherm  occurs  upon  crossing  the  melting 
line.  At  the  triple  point,  the  melting  line  intersects  the  condensation  line  in  two- 
dimensions  just  as  it  does  for  three-dimensional  systems.  Thus,  this  determination  of  the 
triple  point  makes  the  isotherm  search  for  the  melting  line  of  HD  on  BN  all  the  more 
exciting  and  important. 


Table  2-1.   Critical  temperatures  of  monolayers  of  hydrogen  isotopes  adsorbed  onto  BN. 


Critical  Temperature  (Kelvin) 

Isotope 

Layer  2 

Layer  3 

Layer  4 

H2 

11, 5±  1.0 

10.9  ±0.5 

HD 

12.3  10.5 

12.2  ±0.5 

11.8  ±0.5 

D2 

14.1  ±0.6 

13.9±0.5 

Table  2-2.  Critical  temperatures  of  HD  adsorbed  onto  three  substrates. 


HD  Critica 

Temperature  (Kelvin) 

Substrate 

Layer  2 

Layer  3 

Layer  4 

Layer  2 
Triple 
Temp 

BN 

12.3  ±0.5 

12.2  ±0.5 

11  8  ±0.5 

Graphite^ 

11.45'-11.8'' 

12.1  ±0.3'' 

11.8  ±0.4'' 

8.44' -8.58" 

MgO' 

10.5  ±0.3 

12.1  ±0.3 

12.4  ±0.4 

9.95  ±0,1 

'[Liu  92a]  on  graphite  foil,  'heat  capacity,  "isotherms,  ''[Ebe  94c],  "[Vil  91] 
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Table  2-3.  Critical  temperatures  of  H2  adsorbed  onto  three  substrates. 


H2  Critical  Temperature  (Kelvin) 

Substrate 

Layer  2 

Layer  3 

Layer  4 

Triple  Temp 
Layer  2 

BN 

11. 5±  1.0 

10.9  +  0.5 

Graphite^ 

9.31 

10.08 

5.96 

Graphite*" 

lO.OiO.l 

11.0  +  0.5 

6.5  +  O.i 

MgO'^ 

10.05  +  0.05 

10.3  ±0.3 

10.4  +  0.2 

7.20  +  0.05 

D2  on  Graphite 

9.44 

"[Wie  91b]  from  heat  capacity,  "[Per  92]  from  isotherms,  '[Vil  91] 
Table  2-4.  Critical  temperatures  of  D2  adsorbed  onto  three  substrates.     ! 


D2  Critical  Temperature  (Kelvin) 

Substrate 

Layer  2 

Layer  3 

Layer  4 

Triple  Temp 
Layer  2 

BN 

14.1+0.5 

13.9  +  0.5 

Graphite" 

11.8  +  0.2 

12.1+0.3 

11.8  +  0.4 

Graphite'' 

13.10 

13.86 

11.04 

MgO' 

12.4  +  0.4'^ 

13+  l-* 

12.9  +  0.4 

12.1  +  0.2 

'[Liu  92a],  ""[Wie  91],  \W\\  91],  '[Deg  88] 

Table  2-5.  Triple  point  temperatures  of  second  layer  hydrogen  isotopes  adsorbed  onto 
selected  substrates. 


Trip  e  Point  Temperature  (Kelvin) 

Isotope 

Substrate 

H2 

HD 

D2 

Graphite" 

5.96 

'8.44- "8.58'' 

11.04 

MgO' 

7.20  +  0.05 

9.9  +  0.1 

12.1  +  0.2 

D2  on  Graphite" 

5.74 

6.25-6.4 

HD  on  Graphite** 

6.25 

'[Wie  91b],  ''[Vil  91],  '[Liu  92a],  ''[Liu  93]  from  'heat  capacity,  "isotherms 
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2.6  Isosteric  Heat  of  Adsorption 

The  isosteric  heat  of  adsorption  can  be  derived  directly  from  a  pair  of  isotherms 
and  is  defined  as 

^-"^"^^7^  'v....  Equation  2-2 

where  the  fihn  density,  Nf,  and  the  substrate  area.  A,  are  frxed  [Das  75].  In  general,  this 
quantity  is  both  a  function  of  temperature  and  film  density.  Usually,  one  is  interested  in 
the  heat  of  adsorption  for  the  first  adatom  into  the  gaseous  state  of  an  adlayer.  Smce  the 
vapor  pressure  above  the  first  layer  for  hydrogen  at  low  temperatures  is  too  low  to  be 
practically  measured,  the  second  layer  value  is  often  reported  as  in  Table  2-6.  The 
isosteric  heat  is  a  measure  of  the  binding  energy  between  the  adatom  and  the  surface.  If 
the  adatom-adatom  interactions  can  be  neglected  (i.e.,  at  low  coverages),  then  the 
isosteric  heat  is  simply  the  sum  of  the  binding  energy,  s^  ,  and  the  kinetic  energy 

q^^  =  £„  +  a(kf^T)  Equation  2-3 

where  OC  =  1/2  for  a  gas  molecule  hannonically  bound  to  an  adsorption  site  and  CC  =  3/2 
for  a  two-dimensional  ideal  gas  on  a  flat  surface  [Liu  93]. 

Figure  2-18  shows  the  isosteric  heat  of  adsorption  for  the  second,  third,  and  fourth 
layers  of  HD  on  boron  nitride.  Tliese  values  are  extracted  from  two  isotherms  for  which 
the  temperature  difference  may  be  too  large  (  ^  T=  0.33,  0.49  K)  to  be  a  true  differential  as 
defined  in  Equation  2-2.  However,  many  features  are  still  apparent  in  the  data.  Using 
Equation  2-3,  the  binding  energy  for  the  second  layer  of  HD  on  BN  is  found  to 
he  Sg=  120  ±  10  K.  Tliis  result  is  somewhat  lower  than  the  value  of  160  K  reported  for 
HD  on  graphite  [Liu  93]  and  much  larger  than  the  HD-HD  binding  energy  of  37  K. 
Figure  2-18  fiarther  manifests  that  the  change  in  isosteric  heat  as  a  fiinction  of  film  density 
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Figure  2-18.  Isosteric  heat,  q^,  for  the  second,  third,  and  fourth  layers  of  HD  on  BN  for 
two  runs. 
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Figure  2-19.  Isosteric  heat,  q^,  for  the  second,  third,  and  fourth  layers  of  H2  (top)  and  the 
second  layer  of  D2  (bottom)  on  BN.  Lines  are  guides  to  the  eye. 
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is  smaller  by  approximately  20%  than  HD  on  graphite.  These  two  features  are  consistent 
with  the  lack  of  a  large  melting  substep.  Finally,  note  that  the  change  in  the  heat  of 
adsorption  for  the  third  layer  is  similar  to  the  second  layer  (approximately  20%)  and  that 
the  fourth  layer  change  (approximately  10%)  is  much  smaller  than  the  previous  two  layers. 
After  the  fifth  layer  (not  shown),  the  isosteric  heat  of  adsorption  tends  to  a  value  of  144  K 
which  can  be  compared  to  the  three-dimensional  heat  of  vaporization  (134  K)  and  the  heat 
of  sublimation  (153  K). 

Figure  2-19  exhibits  the  isosteric  heat  of  adsorption  for  H2  and  D2  on  BN.  Note 
that  relative  to  HD,  the  heat  is  less  for  H2  and  more  for  D2  as  expected. 

Vidali  has  compiled  a  table  of  potentials  of  first  layer  physical  adsorption  [Vid  91]. 
This  table  lists  the  adsorption  potentials  and  equilibrium  distances  for  a  variety  of  gases  on 
a  number  of  substrates,  including  BN,  graphite,  and  MgO.  Sources  of  information  include 
the  results  of  thermodynamic  (e.g.,  adsorption  isotherms  and  calorimetric  measurements) 
and  LEED  experiments  as  well  as  theoretical  estimates.  While  no  entries  exist  for  the 
hydrogen  isotopes  adsorbed  onto  BN,  the  tabulations  for  the  rare  gases  is  interesting. 
Bare  graphite  adsorbs  Ar,  Ne,  Xe,  Kr,  N2,  and  CH4  more  strongly  than  either  bare  BN  or 
MgO  by  10-20%. 

Table  2-6.  Binding  energy  of  hydrogen  into  the  second  adlayer. 


BN 

Graphite 

MgO 

H2 

100±20K 

120  ±  lOK' 

HD 

130 ± lOK 

177 ±  lOK 

D2 

170 ± lOK 

'  Taken  from  a  graph  in  [Per  92]. 


A  very  simple,  phenomenological  method  of  comparing  the  adsorption  potentials 
of  various  substrates  involves  directly  relating  the  temperatures  required  for  a  given 
adsorbate  vapor  pressure.  For  an  HD  vapor  pressure  of  0.1  Torr  at  the  beginning  of  the 
second  layer,  a  temperature  of  12.6  K  is  required  for  graphite,   12.1  K  for  MgO,  and 
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11.7  K  for  BN.    This  simple  method  emphasizes  the  relatively  weak  adsorption  potential 
of  boron  nitride. 


2.7  Summary 

The  use  of  boron  nitride  as  a  substrate  for  the  translational  phase  study  of  the 
hydrogen  isotopes  has  accomphshed  two  goals:  more  clues  have  been  added  to  the  search 
for  superfluid  hydrogen  in  two  dimensions  and  the  binding  energies  have  been  determined. 
Although  the  precision  of  the  determination  of  these  values  from  the  adsorption  isotherms 
of  this  study  is  less  than  ideal,  the  isotopic  trends  are  consistent  with  expectations  and 
reinforce  the  accuracy  of  each  measurement.  It  is  important  to  note  that  Larher's  method 
has  been  shown  to  overestimate  the  critical  temperatures. 

The  binding  energy  of  hydrogen  isotopes  on  boron  nitride  is  less  than  that  on 
graphite  and  MgO.  However,  this  decreased  attraction  does  not  translate  into  lower 
critical  pomt  temperatures  as  was  hypothesized.  Further,  the  magnitude  of  the  Hquid-soHd 
substep  appears  to  have  been  suppressed  below  those  of  hydrogen  on  either  grapliite  or 
MgO.  It  may  be  surmised  that  the  surface  potential  corrugation,  lattice  spacing  and 
substrate  geometry  play  a  significant  role  here.  In  the  search  for  superfluid  hydrogen, 
much  could  be  learned  about  the  role  of  binding  energies  on  critical  temperatures  through 
the  investigation  of  a  substrate  which  shares  the  geometry  of  graphite  and  boron  nitride 
but  possesses  a  larger  adsorption  potential  than  graphite. 

If  the  translational  phase  diagrams  of  hydrogens  adsorbed  onto  boron  nitride  are 
similar  to  those  of  hydrogens  adsorbed  onto  other  substrates,  then  the  triple  point 
temperatures  are  probably  also  higher  than  on  other  substrates.  However,  until  a  substep 
in  the  isotherms  is  identified,  other  methods  are  necessary  to  determine  the  triple  points 
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and  to  confirm  this  hypothesis.  The  investigation  of  the  triple  points  by  heat  capacity 
techniques  will  probably  have  to  wait  untU  difiBculties  associated  with  the  low  thermal 
conductivity  of  boron  nitride  (approximately  100-1000  tunes  less  than  graphite)  are 
overcome.  One  possible  solution  is  to  add,  in  addition  to  copper  wires,  an  inert  thermal 
link  (i.e.,  ^He)  to  the  system  such  as  was  done  for  the  NMR  apparatus. 

The  use  of  NMR  to  map  out  the  translational  phase  diagram  for  this  system  has 
two  advantages.  First,  NMR  is  not  only  ideally  suited  to  determine  translational  phase 
transitions  but  also  to  identiiy  the  transitional  phases.  Second,  as  opposed  to  the  high 
electrical  conductivity  of  graphite,  the  low  conductivity  of  the  boron  nitride  allows  for 
unhindered  NMR  measurements.  However,  such  a  search  would  require  even  more 
sensitivity  than  that  required  in  the  efforts  to  map  out  the  orientational  phase  diagram,  as 
discussed  in  the  following  chapters  of  this  thesis,  due  to  the  inverse  temperature 
dependence  of  the  magnetic  susceptibihty  (Curie's  Law).  However,  with  the  possibility  of 
superfluid  hydrogen  in  the  balance,  the  search  for  the  triple  point  of  the  isotopes  may  have 
large  rewards  indeed. 


CHAPTER  3 


NMR  APPARATUS  AND  PROCEDURES 

Reading  maketh  a  full  man;  conference  a  ready  man;  and  writing  an  exact  man. 
Of  Studies  Essay  50,  Sir  Francis  Bacon  1561  -  1626 


3.1  Overview 

Historically,  NMR  studies  have  been  performed  on  bulk  materials  [Pak48, 
Sul  67].  Recently,  however,  new  opportunities  for  the  study  of  two-dimensional  systems 
have  emerged  [Har  88,  Cow  93].  Even  one  dimensional  systems  have  begun  to  be 
studied  to  test  the  effects  of  reduced  dimensionality  [Ave  92].  However,  the  signal  from  a 
single  nucleus  is  very  small  and  provides  a  natural  limit  as  to  the  number  of  atoms 
observable  by  standard  NMR  techniques. 

An  NMR  signal  is  directly  proportional  to  the  number  of  atoms  being  investigated. 
Since  the  number  of  atoms  in  a  two-dimensional  system  is  proportional  to  the  number  of 
atoms  in  a  similar  three-dimensional  volume  to  the  two-thirds  power,  or 

N^d'^NTd  Equation  3-1 


a  similar  reduction  in  the  NMR  absorption  amplitude  is  expected.  In  other  words,  the 
signal-to-noise  ratio  of  two-dimensional  systems  is  inherently  low.  Fortunately,  hydrogen 
is  the  second  most  sensitive  nucleus,  after  tritium,  in  the  entire  periodic  table.  '^F  and  He'' 
are  the  third  and  fourth  most  sensitive  elements,  respectively. 

The  intensity  of  an  NMR  signal  also  depends  on  the  polarization  of  the  nuclei 
under  study,  which  in  turn  depends  upon  the  temperature  of  the  system  and  the  local 
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magnetic  field.     This  dependence  can  be  seen  by  starting  with  the  expression  for  the 
polarization  of  a  spin  system;  i.e., 

H+  -  ^L         1  -  g  *''^ 

P  =  -—^  = ZK^  Equation  3-2 

•"       -       1  +  g  *6r 

and  noting  that  the  denominator  is  essentially  unity.  One  can  then  expand  the  exponential 
in  the  numerator. 


. }  Equation  3-3 


na}„      1   na)„ 

P=l-{1- — -+ -.. 

k^T     2\kjj 

hco„ 
and  obtain  an  expression  for  the  spm  polarization  vahd  for  small  "r~:r ,  namely, 

P «  .  Equation  3-4 

Since  the  NMR  fi:equency,  cOo,  of  a  nucleus  is  directly  related  to  the  magnetic  field, 

B,  by  the  gyromagnetic  ratio,  ;' ,  through  the  expression, 

(Oo  =  yB,  Equation  3-5 

we  then  have  an  expression  for  the  polarization  in  terms  of  5  and  T: 

hy    B 


Ik,   T 


Equation  3-6 


For  small  B/T,  that  is,  at  high  temperatures  or  low  fields,  the  polarization  is 
approximately  hnear  in  B/T.  For  large  B/T  ,  that  is,  at  low  temperatures  and  high  fields, 
this  hnear  approximation  breaks  down  and  the  exponential  behavior  becomes  evident  as 
the  polarization  approaches  unity.  For  a  field  strength  of  6.29  Tesla,  corresponding  to  a 
resonance  frequency  of  268  Megahertz  for  hydrogen,  the  exponential  behavior  becomes 
apparent  for  temperatures  below  approximately  0.5  K  as  the  premise  for  the  linear 
approximation  of  Equation  33  breaks  down. 

It  is  important  to  estimate  the  signal-to-noise  ratio  before  considering  an 
experiment   lest   one   attempts   a   theoretically  impossible   ex'periment.      According   to 
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Abragam  [Abr  61],  the  signal-to-noise  ratio  for  an  NMR  measurement  can  be  shown  to 
be: 


S         TT   \  T]  f    V 


N     /[2VAW 


\ 


V. 


\ 


1/2 


QZolZo^  T^\  Equations-? 


where  f  is  the  noise  from  the  NMR  apparatus  excluding  the  coil;  v  is  the  frequency;  Av  is 
the  bandwidth;  Q  is  the  quality  factor  of  the  coil;  Xo  i^  the  nuclear  magnetic  susceptibihty 
of  the  nucleus  being  probed;  B  is  the  magnetic  field;  and  Vs  is  the  volume  of  the  sample. 
From  this  equation,  it  can  be  seen  that  the  NMR  signal  is  directly  proportional  to  the 
fiUing  factor  of  the  coil,  TJ  ,  which  is  usually  on  the  order  of  unity  for  a  typical  three- 
dimensional  experiment. 

However,  for  two-dimensional  experiments  such  as  this  one,  the  fill  factor  is 
approximately  10"*  lower.  The  small  filling  factor  arises  from  the  fact  that  the  apparent 
density  of  the  boron  nitride  powder  is  very  low  and  the  hydrogen  molecules  occupy  only 
the  surface  of  the  powdered  platelets.  The  actual  number  of  hydrogen  molecules  in  the 
sample  cell  is  determined  from  the  adsorption  isotherms  discussed  earher.  The  values  of 
the  variables  to  be  used  in  Equation  3-7  have  been  estimated  from  the  experimental 
conditions  to  be: 

7  «  10^ ,    V  «  268A///r,     A  v  =  03Hz,    Q^50, 

Equation  3-8 
B^  «  62,900gauss,    V^  «  3cm' ,     T«AK. 

The  nuclear  magnetic  susceptibility  of  the  system  may  be  calculated  from: 

N{y  nf  !{!  +  {) 
X  a " ^Tr ~  1-7*10'"   (M.K.S.  units)  Equation  3-9 

where  we  have  used  the  fact  that  N  ~  10*'  spins,  y ^  h  =  4.47*10'^'  Joules/Tesia,  and  k^T 
=  4.14*10''  Joules  at  4  K.  Wliile  the  signal-to-noise  ratio  may  now  be  calculated  from 
Equation  3-7,  it  is  more  instmctional  to  utilize  an  equation  which  ex-plicitly  shows  the 
dependence  of  the  signal-to-noise  ratio  relative  to  the  ex-perimental  NMR  standard  of  the 
protons  in  one  cubic  centimeter  of  water  at  room  temperature  in  a  field  of  one  Tesla: 
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where  Np  is  the  number  of  protons  in  a  cubic  centimeter  of  water  at  room  temperature, 
7*10-"',  and  Vc  is  the  volume  of  the  coU.  Using  the  values  for  our  two-dimensional 
system,  we  find 

—  xino^/f.  Equation  3-11 

While  at  first  this  may  seem  more  than  sufficient,  in  most  cases  the  theoretical 
values  are  rarely  achieved.  For  comparison,  the  calculated  NMR  signal-to-noise  ratio  for 
protons  in  water  at  room  temperature  is  approximately  10^//  [Abr  6 1]    but  ratios  of  lO'- 

10^  are  usually  achieved.  We  also  note  that  the  calculated  signal  for  our  two-dimensional 
experiment  is  approximately  1000  times  smaller  than  that  of  bulk  hydrogen  experiments, 
where  the  fill  factor  would  be  between  0.1  and  1  depending  upon  the  construction  of  the 
sample  cell. 

Another  calculation  one  can  make  concerning  the  experimental  feasibility  of  the 
experiment  is  the  determination  of  the  voltage  across  the  NMR  coil: 

£•      =-7j4mo^  ^  B^AQ*\0~^  xl.S'^lO''  volts  Equation  3-12 

where  A  is  the  area  of  the  coil  and  the  other  parameters  are  as  above.  This  low  voltage 
demands  the  use  of  higli  ampHfication  and  sensitive  measurement  techniques.  According 
to  Souers  [Sou  86],  the  limit  of  sensitivity  (i.e.,  the  state  of  the  art)  of  an  NMR  signal  for 
hydrogen  at  4  K  corresponds  to  a  signal-to-noise  ratio  of  4:1  for  a  spin  density  of  lO" 
spins  per  cubic  centimeter  (bulk  hydrogen  has  a  spin  density  of  10'  per  cc).  Data 
acquired  for  this  thesis  yielded  a  similar  signal-to-noise  ratio  which  improved  slightly  as 
the  temperature  decreased.    However,  the  expected  1/T  gain  in  the  signal-to-noise  ratio 
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was  realized  only  for  temperatures  down  to  1  K.  Below  1  K,  the  magnetization  was  found 
to  be  independent  of  temperature. 

As  the  technology  for  high  magnetic  fields  improves,  we  will  be  able  to  study  ever 
smaller  systems  of  spins.  Organizations  such  as  the  National  High  Magnetic  Field 
Laboratory  in  Tallahassee,  Florida,  and  the  European  collaborations  in  Grenoble,  France, 
continue  to  produce  ever  higher  fields.  Furthermore,  recent  advances  in  semiconductor 
technology  allows  for  the  use  of  spectrometers  immersed  in  the  low  temperature 
enviroiunent.  Since  thermal  noise,  also  known  as  Johnson  or  "shot,"  noise,  is  often  a 
limiting  factor  in  experiments  with  low  signal,  operating  the  core  electronics  at  low 
temperatures  is  a  great  advantage.  Recently,  an  NQR  experiment  utilized  new  GaAs 
transistors  in  a  marginal  oscillator  circuit  operating  at  high  milliKelvin  temperatures  and 
hundreds  of  Megahertz.  Tlie  oscillator  was  located  inches  away  from  the  MnSb  sample 
inside  the  resonance  coil  [And  92b].  For  this  quantum  thermometer  experiment,  the  low 
temperature,  marginal  oscillator  spectrometer  worked  extremely  well. 


3.2  Experimental  Apparatus 

The  heart  of  the  nuclear  magnetic  resonance  system  used  for  these  experiments 
consisted  of  a  continuous  wave  (CW)  NMR  spectrometer.  Tliis  arrangement  was  first 
used  by  Purcell,  Torrey  and  Pound  [Pur  46]  to  measure  the  first  magnetic  resonance  of 
protons  in  paraffin  in  1946.  Early  NMR  experiments  using  this  technique  were  soon 
replaced  in  most  research  by  pulsed  NMR  techniques.  These  pulsed  systems  have  many 
advantages  under  a  wide  variety  of  circumstances.  However,  pulsed  systems  yield  data  in 
the  time  domain,  which  is  the  Fourier  transform  of  the  frequency  domain  in  whicli  the  CW 
NMR  line  shape  is  measured.  In  a  world  of  infinite  signal-to-noise  ratios,  such  as 
mathematics,  these  two  domains  are  equivalent.    In  the  real  world  of  small  signals  and 
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large  noise  it  can  be  advantageous  to  record  the  data  in  the  domain  of  interest.  In  our 
case,  the  NMR  hne  shape  (the  resonance  response  as  a  function  of  magnetic  field)  was 
examined  to  obtain  information  concerning  the  orientational  ordering  in  the  fi"equency 
domam. 


3.2.1   Spectrometer  Schematic  Diagram 

The  CW  NMR  spectrometer  utilized  for  this  thesis  consisted  of  two  parts:  the 
magnetic  field  sweep  and  the  RF  absorption  detector,  connected  in  theory  only  through 
the  spins  of  the  hydrogen  nuclei  imder  investigation.  The  sample  cell,  detection  coil,  and 
sample  can  were  boUed  onto  the  cold  finger  of  a  dilution  refiigerator  resting  inside  a 
6.7  Tesla  superconducting  magnet. 

The  magnetic  field  sweep  consisted  of  four  main  components:  a  lock- in  amplifier;  a 
programmable  power  supply;  a  summing  ampUfier;  and  a  modulation  magnet.  An  audio 
fi"equency  modulation  signal  was  generated  by  the  lock-in  amplifier  and  sent  to  a  summing 
amplifier  where  it  was  combined  with  a  linear  ramp  (generated  by  an  IBM  compatible  PC) 
and  sent  to  the  current  feed  of  a  programmable  power  supply.  Tlie  power  supply 
generated  the  current  for  the  vertical,  superconducting  modulation  magnet  located 
between  the  sample  cell  and  the  main  magnet. 

Tlie  central  tenet  of  the  experimental  endeavor  concerned  the  width  of  the  NMR 
line.  Tims,  the  calibration  of  the  modulation  magnet  was  critical.  Calibration  of  the 
modulation  magnet  was  achieved  by  varying  the  fi'equency  of  the  radio  frequency  signal 
(known  to  1  part  in  lO'')  and  measuring  the  change  in  current  required  to  center  the  NMR 
signal  in  the  field  sweep.  Tliis  procedure  led  to  a  frequency-to-current  ratio  of  50  kHz  per 
amp.   Tlie  well  known  gyromagnetic  ratio  of  hydrogen,  42.5759  Megahertz  per  Tesla,  can 
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then  be  used  to  convert  current  to  magnetic  field.  The  current-to-field  ratio  was  found  to 
be  0.085  Amp  per  Gauss. 

In  practice,  the  much  larger  '^F  NMR  signal  fi"om  the  sample  holder  was  used  to 
obtam  a  more  precise  figure.  Fortuitously,  the  gyromagnetic  ratio  of  fluorine, 
40.0541  Megahertz  per  Tesla  [CRC  91],  is  only  sUghtly  smaller  than  hydrogen.  Thus,  the 
fluorine  NMR  signal  required  only  a  sUghtly  larger  static  magnetic  field  fi'om  the  main 
superconducting  magnet  (approximately  59  versus  55  amps)  for  the  same  resonant 
fi'equency,  268  Megahertz.  Furthermore,  the  NMR  hne  width  of  fluorine  in  Teflon 
(tetrafluoroethylene)  is  similar  to  the  line  width  of  molecular  hydrogen.  Thus,  another 
sample  holder  was  manufactured  out  of  Teflon  to  replace  the  original  KelF  sample  holder. 
While  KelF  is  a  commonly  used  material  in  low  temperature  physics,  a  fluorine  NMR 
signal  was  not  found  under  our  experimental  conditions  This  lack  of  signal  was  at  first 
surprising  since  KelF  is  chemically  very  similar  to  Teflon  with  the  exception  that  one  out 
of  every  four  fluorine  atoms  is  replaced  with  a  chlorine  atom. 

Tlie  radio  fi'equency  spectrometer  consisted  of  five  main  components:  an  RF  signal 
generator;  a  hybrid  tee;  an  RF  amplifier;  audio  amplifiers;  and  a  lock-in  amplifier.  Tlie 
inherent  low  absolute  signal  of  a  two-dimensional  experiment  necessitates  the  use  of  a 
high  sensitivity  spectrometer.  The  hybrid,  or  "magic,"  tee  configuration  of  a  CW  NMR 
spectrometer  is  the  most  sensitive  of  the  various  CW  NMR  configurations  and  is 
illustrated  in  Figure  3-1.  As  in  most  CW  NMR  arrangements,  the  NMR  signal  is  a  small 
perturbation  on  a  larger  RF  carrier  signal.  Wliile  measuring  a  small  change  in  a  relatively 
large  signal  is  inherently  difficult,  measuring  the  same  small  change  fi'om  a  null  is  much 
easier.  Thus,  the  RF  energy  coming  into  the  hybrid  tee  from  the  RF  generator  is  split  into 
two  parts. 
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Figure  3-1    Schematic  diagram  of  the  Continuous  Wave  Nuclear  Magnetic  Resonance 
(CW  NMR)  spectrometer  utilizing  a  hybrid  tee. 


One  half  of  the  RF  energy  resonates  the  NMR  coil  at  port  (b)  while  the  phase  of 
the  other  half  is  shifted  by  180  degrees  and  open  to  a  real  resistive  load  at  port  (c).  The 
high  sensitivity  of  the  hybrid  tee  is  achieved  by  opposing  the  approximately  50  ohm 
impedance  tank  circuit  with  a  resistance  (real  impedance)  of  50  ohms.  The  recombination 
of  these  two  signals  with  opposite  phases  generally  resuUs  in  a  small  signal  and,  in  fact,  no 
signal  if  the  two  loads  are  perfectly  matched.  Upon  sweepmg  the  magnetic  field  tlirough 
resonance,  both  the  inductance  and  resistance  of  the  NMK  coil  are  altered  very  shghtly  by 
the  nuclear  spins  inside  the  coil.  This  change  in  unpedance  alters  the  dehcate  balance  (null 
status)  of  the  hybrid  tee,  is  picked  up  by  the  amphfiers,  and  sent  on  to  the  recording 
equipment. 

In  practice,  the  hybrid  tee  is  tuned  just  off  of  the  null  to  allow  a  constant  carrier 
signal  (voltage)  which  establishes  a  very  important  background  reference.  Usually,  the 
hybrid  tee  is  further  timed  to  measure  the  change  in  circuit  resistance  rather  than  the 
change  in  inductance.  This  has  the  effect  of  measuring  the  absorption  signal  rather  than 
the  dispersion  signal.  However,  tliis  is  merely  a  convention  as  the  two  signals  are  directly 
related  through  the  Kramers-Kronig  relations.  Some  experimenters  choose  to  look  at  the 
signal  "m  quadrature";  that  is,  to  measure  both  the  in-phase  and  out-of-phase  signal 
simultaneously.  While  this  method  has  the  advantage  that  any  error  in  phase  can  then  be 
corrected  througli  the  use  of  the  Kramers-Kronig  relations  after  the  data  has  been 
collected,  such  fine  tuning  of  the  signal  measured  in  this  ex'periment  could  not  be 
performed  due  to  the  low  signal-to-noise  difficulty.  In  our  case,  sample  runs  proved  that 
the  inherent  added  complexity  of  quadrature  measurements  were  not  worth  the  extra 
efforts  since  the  signal-to-noise  ratio  was  not  improved  significantly. 

Tlie  RF  amplifier,  with  an  ampUfication  factor  of  approximately  250,  picks  up  the 
imbalance  from  the  fourth  side  of  the  hybrid  tee.  An  RF  diode  after  the  RF  amplifier 
converts  the  AC  signal  into  DC  and  eliminates  the  need  for  a  mixer  and  its  associated 
phase  from  the  experiment.   However,  the  cost  associated  with  this  simplification  is  some 
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loss  of  information.  The  DC  signal  can  not  be  passed  through  a  mixer  to  separate  the 
absorption  and  dispersion  signals.  Thus,  the  tank  cncuit  was  often  tuned  to  insure 
measurement  of  only  the  absorption  signal.  This  magic  tee/RF  diode  method  may  be 
compared  with  the  Q-meter  method  [Fuk  87]  which  measures  only  the  absorption  signal 
but  has  less  sensitivity. 

After  the  RF  diode  and  some  amphfication,  a  band-pass  filter  (10-lOOOHz)  was 
used  to  separate  the  audio  fi'equency  of  the  NMR  signal  from  the  RF  carrier  signal.  A 
second  audio  amplifier  was  then  used  to  raise  the  signal  voltage  to  an  optimal  level  for  the 
lock-in  amplifier.  This  very  special  amplifier  compared  the  amplitude,  frequency,  and 
phase  of  the  incoming  signal  with  that  of  the  audio  modulation  it  generated  and  relayed  to 
the  modulation  magnet. 

One  important  consequence  of  this  method  is  that  the  derivative  of  the  NMR  hne 
shape  was  actually  measured.  Such  derivative  measurement  is  a  consequence  of  using  the 
lock-in  amplifier  as  the  detecting  instrument.  While  the  field  was  ramped  over  three  times 
the  line  width,  the  audio  frequency  modulation  covered  only  a  small  portion  (five  percent 
or  less)  of  the  Une  shape  per  cycle.  Although  narrow  hnes  may  be  measured  without  a 
ramp  or  a  lock- in  amplifier,  the  direct  observation  of  wide  line  shapes  is  more  complex  as 
the  large  inductance  of  the  ramping  (modulation)  magnet  limits  the  sweep  width  or  time. 
One  benefit  of  derivative  measurement  is  that  the  derivative  of  the  Une  shape  is  more 
sensitive  to  small  changes  in  the  hne  shape  than  the  line  shape  itself  All  line  shapes 
presented  here  have  been  numerically  integrated  by  a  computer. 

3.2.2  Tuning  of  the  Tank  Circuit 

To  analyze  the  tank  circuit  in  detail,  we  follow  the  treatment  given  by  Abragam 
[Abr  61].  The  tank  circuit,  with  the  sample  inside  the  cylindrical  coil,  has  an  inductance  L, 


a  volume  Vc  =  Vs  /  q  ,  and  a  current, 

i  =  I cosco  t  ^  Kq{I e' '"' ) 


Equation  3-1! 


which  produces  an  RF  field  Hi(t)  =  2  Hi  coscot  inside  the  coil.    The  flux  across  the  coil, 

O  =  Li,  produces  a  sample  magnetization,  M.   Following  the  standard  procedure,  both  the 

complex  susceptibiUty,  x  ^  ^^'^  magnetization  can  be  broken  into  real,  x  '  ^^'^  ^^^  ^^d 

imaginary  part,  x  "  ^nd  My,  respectively,  where, 

M,  =  Re(2  X  Hi  e^')  Equation  3- 14 

and  the  susceptibUity  is  vanishingly  small  except  near  the  resonance.   The  flux,  O,  througli 
the  coil  is  the  product  of  the  induction,  B  =  H  -4;tM,  and  the  area  of  the  coil,  A. 


Figure  3-2  Model  tank  circuit  used  to  analyze  the  real  tank  circuit  used  in  the  NMR 
experiment. 

A  simplified  tuned  tank  circuit  may  be  modeled  as  shown  in  Figure  3-2.  In  this 
model,  all  capacitance,  including  all  stray  capacitance,  is  assumed,  for  simpUcity,  to  be  in 
parallel  with  the  inductance.  At  high  frequencies  (RF  and  higher)  stray  capacitance  is 
omnipresent  and  comes  from  any  and  all  pairs  of  conductors  (e.g.,  the  windings  of  the  coil 
or  between  the  coil  and  sample  can).  In  practice,  both  parallel  and  series  capacitance  are 
present  and  were  corrected  for  through  the  use  of  a  room  temperature  tuning  box  which 
contained  tunable  parallel   and   series  capacitors.      Tims,   the  impedance   of  the   low 
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temperature  portion  of  the  tank  circuit  was  not  measured  but  rather  tuned  at  room 
temperature  to  50  ohms.  The  tuning  characteristics  of  the  tank  circuit  and  the  role  of  the 
two  capacitors  are  illustrated  in  Figure  3-3. 

Low  temperature  capacitors  were  not  utilized  due  to  the  compUcations  in  tuning 
such  an  arrangement.  However,  the  improvement  in  signal-to-noise  from  using  low 
temperature  capacitors  (estimated  at  a  factor  of  10)  can  make  the  tuning  difficulties  (and 
the  associated  pain)  worthwhile.  Such  an  arrangement  is  currently  being  pursued  as  a 
refinement  of  the  experiment. 

Using  the  above  model,  both  the  shunt  resistance , 

R=I}(o^  Ir  Equation3-15 

and  the  quaUty  factor  of  the  circuit, 

Q  =  (o  LI r  =  Rl CO  L  Equation  3-16 

can  then  be  defined.  The  quaUty  factor,  Q  ,  is  defined  as  the  ratio  of  the  energy  stored  to 
the  energy  dissipated  per  cycle  in  the  LCR  circuit.  The  queue  is  usually  a  number  between 
50  and  1000  for  RF  circuits,  much  less  than  for  microwave  circuits.  The  impedance  of  this 
tank  circuit  model  can  then  be  deduced  from  the  combination  of  the  parallel  and  series 

impedance, 

1  1 

—  = ; — : +  iCco  Equation  3- 1 7 

Z      r+jLco{\  +  Anff)     ■'  ^ 

At  resonance,  LCco^  =  1,  the  impedance  reduces  to 


Z  =  R 


l  +  47rrj-J/Q 


Equation  3-18 


This  expression  can  be  fiirther  reduced  by  realizing  tliat  the  second  term  in  the 

above  numerator  is  much  less  than  unity  so  that  the  quotient  can  be  e.xpanded  by  the 

binomial  theorem.  Tlie  complex  impedance  then  becomes 

Z~R[\-J4;r  qQX]  =  R[[-47r  if  QX"-JK  ff  QX'].         Equation  3-19 
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Figure  3-3   Tuning  characteristics  of  the  CW  NMR  tank  circuit  showing  the  effect  of  both 
the  series  and  parallel  capacitance. 


Great  care  was  taken  to  reduce  the  presence  of  additional  hydrogen  atoms  which 
would  add  to  the  NMR  signal  and  hmit  the  experimental  resolution.  However,  it  was  not 
possible  to  remove  hydrogen  below  the  precision  of  this  experiment;  a  background  signal 
resisted  all  efforts  of  removal.  In  practice,  the  background  signal  was  measured  before  the 
molecular  hydrogen  sample  was  introduced  into  the  sample  chamber.  Furthermore,  during 
the  course  of  an  experimental  run,  the  spin-one  ortho-hydrogen  species  decays  into  spin- 
zero  para-hydrogen.  Thus,  as  the  experiment  progresses,  the  NMR  signal  decays  towards 
the  background  signal.  Wlien  the  signal  from  the  hydrogen  dosed  sample  cell  returned  to 
the  pre-dosed  level,  the  experimental  run  was  ended  and  the  cell  allowed  to  return  to  room 
temperature  to  remove  all  of  the  residual  molecular  hydrogen. 


3.3   Ex^jerimental  Procedure 


Tlie  greatest  difficulty  of  the  experiment  concerned  the  low  signal-to-noise  ratio 
inherent  in  any  two-dimensional  ex-periment.    Even  though  excellent  books  exist  on  the 
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subject  of  low  temperature  experimental  procedures  [e.g.,  White,  PobeU,  and  Lousanama], 
this  section  is  included  to  guide  future  experiments  with  similar  problems.  As  has  been 
often  bemoaned,  "...the  identification  of  the  problem  is  often  one-half  of  the  solution." 

Although  the  previous  experiment  used  the  same  low  temperature  rig,  the  NMR 
set-up  was  not  adequate  for  the  significantly  lower  signal-to-noise  ratio  of  a  two- 
dimensional  experiment.  The  previous  experiments  on  this  rig  studied  three-  and  quasi- 
three  dimensional  hydrogen  systems  and  thus  required  less  sensitivity.  Therefore,  problems 
which  were  present,  but  not  detectable  with  the  lower  gain  required  for  the  earher 
experiments,  were  uncovered.  Furthermore,  the  preceding  experiments  examined  the 
narrower  orientationally-disordered  phase  and  many  of  the  problems  which  were 
uncovered  resulted  from  the  larger  field  sweeps.  Fortimately,  the  ortho-  to  para-hydrogen 
conversion  rate  was  the  lowest  yet  recorded  which  allowed  for  longer  signal  averaging 
times. 

At  temperatures  below  approximately  1  K,  the  width  of  the  Une  shape  increased 
from  approximately  10  to  40  gauss.  An  ideal  magnetic  field  sweep  for  CW  NMR  is  two 
to  four  times  the  width  of  the  line  shape.  Such  a  baseline  to  line  shape  ratio  insures  that 
both  the  signal  and  the  baseline  are  well  determined.  One  difficulty  in  measuring  wide  Une 
shapes  is  background  drift.  During  the  experiment,  it  was  discovered  that  the  RF  detector 
was  sensitive  to  wide  sweeps  of  the  magnetic  field;  that  is,  the  background  would  follow 
the  large  sweeps  in  a  complex  fashion.  The  undesired  linear  (and  largest)  response  was 
minimized  by  averaging  the  signals  from  both  of  the  two  sweep  directions.  Fortunately, 
we  were  able  to  automate  this  procedure  and  thereby  minimize  the  linear  and  other  odd 
terms  of  the  pick-up  response.  However,  the  even  terms  of  the  pick-up  response  were  still 
both  present  and  noticeable.  After  much  investigation  into  electrical  paths  other  than  the 
electromagnetic  response  of  the  sample,  possible  ground  loops  for  example,  we  decided 
that  the  origin  of  the  sweep  pick-up  was  vibrational. 
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There  exists  two  possible  origins  of  the  vibrations:  mechanical  and 
electromagnetic.  Mechanically,  the  resistive  leads  to  the  superconducting  modulation 
magnet  produced  significant  joule  heatmg  for  currents  greater  than  3  amps.  This  heating 
had  two  eflfects:  an  increase  in  the  hquid  helium  bath  boil-oflfrate  and  the  detixniag  of  the 
tank  circuit.  The  rate  of  the  boil-oflFnot  ordy  creates  white  noise  but  also  ahers  the  tuning 
of  the  tank  circuit  since  the  impedance  of  the  cryogenic  coa.xial  cable  is  temperature 
sensitive  even  though  the  cable  resides  inside  an  evacuated  tube  which  is  part  of  the 
vacuum  can.  This  case  was  easily  identified  by  the  signature  dependence  of  joule  heating 
upon  the  currents  in  the  modulation  leads.  Joule  heating  is  independent  of  direction  and 
proportional  to  the  square  of  the  current,  that  is,  the  second  harmonic  of  the  sweep  rate,  a 
non-hnear  response.  This  theory  was  extensively  tested  and,  indeed,  confirmed  by 
simulating  the  effect  with  a  non-inductive  (real  resistance)  heater  and  measuring  both  the 
pressure  of  the  hehum  vapor  above  the  bath  and  the  background  noise  out  of  the  lock-in 
amplifier.  Tlie  old  leads  were  removed  and  new,  less  resistive  leads  were  installed.  Even 
so,  the  noise  was  still  directly  proportional  to  the  level  of  the  hehum  bath  in  the  cryostat 
with  the  lowest  noise  level  occurring  at  low  helium  bath  levels.  Thus,  during  data 
acquisition  the  bath  level  was  rarely  refilled  to  capacity  but  rather  kept  at  a  low  level  near 
the  top  of  the  main  magnet  near  the  quench  danger  zone! 

In  terms  of  electromagnetic  forces,  Ampere's  law  (the  magnetic  force  on  a  wire  is 
proportional  to  the  current  in  the  wire)  must  be  taken  into  account.  For  AC  currents,  the 
forces  on  the  wires  are  vibrational.  Tlie  large  frequency  sweep  required  a  large  DC 
current  through  the  modulation  magnet.  Tims,  the  electromagnetic  forces  were 
proportional  to  the  width  of  the  field  sweep.  Larger  DC  currents  through  the  modulation 
magnet  leads  interact  more  strongly  with  the  constant  AC  magnetic  field  causing  vibration 
of  the  magnet  leads.  Tliis  effect  was  minimized  by  further  twisting  the  pair  of  leads  about 
each  other  and  taping  all  leads  to  rigid  surfaces.  Although  using  twisted  pairs  is  standard 
practice  in  low  temperature  physics,  the  sensitivity  of  the  experiment  required  more  than 
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usual  care.  In  fact,  even  after  the  improvements,  which  reduced  the  noise  floor  by  more 
than  a  factor  often,  the  NMR  rig  made  a  great  accelerometer  and  picked  up  many  noises 
near  the  audio  modulation  frequency  of  111  Hz  (a  frequency  found  to  minimize  this 
effect),  including  voices  and  pumps.  Consequently,  the  quietest  data  was  taken  at  night  in 
the  automated  mode.  A  definite  improvement  in  the  signal-to-noise  ratio  of  the 
experiment  would  be  to  either  replace  the  analog  lock-in  with  one  based  on  digital 
technology  or  to  replace  the  analog  prefilters  with  digital  filters  which  have  a  higher 
dynamic  range.  The  effect  of  these  combined  difficulties  was  to  limit  the  ramp  width 
which  resulted  in  a  poorly  measured  background  for  the  widest  NMR  hne  shapes. 


CHAPTER  4 

NMR  RESULTS  AND  DATA  ANALYSIS 

Hypotheses  nonfmgo 
I  feign  no  hypotheses. 
Sir  Isaac  Newton 

4.1  Overview 

The  NMR  data  was  acquired  using  a  CW  UHF  NMR  bridge  spectrometer 
operating  at  268  MHz  and  interfaced  via  an  IEEE  488  bus  to  an  IBM  compatible  PC.  A 
paper  strip  chart  recorder  was  used  for  back-up  purposes.  A  file  of  the  derivative  line 
shape  was  recorded  on  the  computer  for  each  direction  of  the  individual  magnetic  field 
sweeps.  Included  at  the  end  of  the  file  was  a  list  of  experimental  parameters  including  the 
beginning  and  ending  times  and  temperatures,  ramp  voltages,  RF  power  and  frequency, 
and  lock-in  settings  (modulation  phase;  frequency  and  amplitude;  time  constant; 
sensitivity;  dynamic  range;  filters;  and  offset).  These  values  were  also  obtained  over  the 
IEEE  bus  which  provided  the  connection  for  most  instruments  to  the  PC.  The  PC 
software  changed  these  settings  periodically  in  a  prescribed  manner  for  the  purpose  of 
consistency  checks. 

A  typical  experimental  run  involved  recording  500-3000  files  over  the  course  of  5- 
20  days  while  sweeping  the  temperature  from  0.1-10  K.  Data  was  gathered  for  both 
directions  of  each  CW  NMR  sweep  in  approximately  25  minutes.  The  names  of  these  files 
followed  a  sequential  code  which  was  generated  by  the  PC  software  for  fijture  reference. 
All    software    utilized    for    these    experiments    was    developed    in-house    using    the 
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LabWindows'  programming  enviromnent.  The  CW  ramp  signal  was  also  generated  by  the 
PC  through  the  use  of  a  National  Instruments  AT-MIO-16  digital-to-analogue  (DAC) 
conversion  board.  The  computer  control  of  the  CW  NMR  recording  apparatus  allowed 
for  signal  averaging  and  automated  operation  of  the  experiment  for  as  much  as  23  hours 
per  day.  In  fact,  to  attam  the  highest  signal-to-noise  ratio  obtainable,  data  was  gathered  at 
all  times  during  the  experimental  run  except  during  Uquid  hehum  transfers,  when  the  noise 
of  the  transfer  overwhelmed  the  NMR  signal.  This  experimental  downtime  was  used  to 
back  up  the  data  files. 

As  discussed  above,  the  adsorption  isotherms  were  used  to  measure  the  number  of 
atoms  adsorbed  per  layer  for  a  given  boron  nitride  sample  as  well  as  to  determine  the 
isosteric  heat  of  adsorption.  The  boron  nitride  sample,  already  mside  the  Teflon  holder, 
was  then  moved  to  the  low  temperature  NMR  apparatus.  During  this  process,  a  small 
portion  of  the  fine  boron  nitride  powder  remained  in  the  adsorption  isotherm  apparatus  or 
was  otherwise  lost.  For  this  reason,  additional  adsorption  isotherms  of  Argon  at  77  K 
were  performed  on  the  same  sample  inside  the  NMR  sample  can.  Figure  4-1  provides  a 
schematic  diagram  of  the  low  temperature  sample  can  with  the  NMR  coil  surrounding  the 
Teflon  sample  cell.  Tlie  sample  cell  was  designed  to  meet  two  guidelines:  to  allow  for 
adsorption  of  the  hydrogen  onto  the  boron  nitride  in  equilibrium  and  to  preferentially 
measure  that  part  of  the  hydrogen  sample  adsorbed  onto  the  surface  of  the  boron  nitride. 
Both  goals  were  achieved  by  placing  the  Teflon  sample  cell  inside  the  NMR  coil  which 
was  itself  mside  a  much  larger  brass  sample  can  as  shown  in  Figure  4-1.  Under  certain 
and  very  undesirable  conditions,  the  hydrogen  sample  could  adhere  to  the  brass  walls  of 
the  sample  can  instead  of  onto  the  boron  nhride.  However,  m  this  case  the  NMR  would 
provide  evidence  of  the  wandering  hydrogen  molecules:  the 


'  LabWindows'"'"  is  the  copyright  of  National  Instruments  Corporation,  6504  Bridge  Point  Parkway, 
Austin,  TX  78730-5093. 


91 


To 

Mixing 
Chamber 


Indium 
Seal 


RF  Feed-thru 


Hz/'He  FillLme 


^ 


v^ 


„  A  A  A  „ 


TL 


IT 


H 


H 


H 


U 


wvv 


TL 


RFCoil 


Thermal  Lfaik 
(fine  Cu  wires) 


Sample  Cell 
(teflon) 


Room  Temperature 
Needle  Valve 


Stainless  Steel 
Pin 


Rubber  O-ring 


Vacuum  Port 


Figure  4-1.  The  experimental  cell  arrangement.  Note  the  "cell  within  a  can"  design  where 
the  Teflon  sample  cell  fits  inside  the  brass  NMR  sample  can. 
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amplitude  of  the  hydrogen  NMR  signal  would  be  greatly  diminished^.  Great  care  was 
undertaken  to  maximize  this  design  characteristic  because  it  was  thought  that  NMR  could 
be  used  to  measure  the  two-dimensional  hquid-sohd  transition  thought  to  occur  between  7 
and  14  K.  However,  above  7  K  the  NMR  signal  amphtude  dropped  below  the 
background  signal  causing  this  endeavor  to  be  dropped  in  favor  of  the  lower  temperature 
work  untU  the  problem  with  the  background  could  be  circumrouted.  This  "sample  cell 
within  a  sample  can"  design  is  to  be  contrasted  with  the  more  common  experimental 
arrangement  used  for  the  bulk  studies  where  the  Teflon  (or  Kel-F)  sample  cell  doubles  as 
the  sample  can  and  all  hydrogen  remains  mside  the  NMR  coil. 

Thermal  contact  between  the  cold  finger  of  the  dilution  refi-igerator  and  the 
hydrogen  on  the  surface  of  the  boron  nitride  platelets  was  achieved  by  adding  hehum  to 
the  sample  can.  At  the  higher  sample  cell  loading  temperatures,  14-20  K,  the  hehum 
adsorption  is  minimal  and  the  heUum  provides  gaseous  contact  directly  between  the  walls 
of  the  sample  can  and  the  boron  nitride  platelets.  However,  the  mechanism  for  thermal 
contact  depended  upon  the  sample  temperature.  At  lower  temperatures,  below  1  K, 
where  the  vapor  pressure  above  even  bulk  hehum  is  minimal,  the  hehum  was  assumed  to 
adsorb  onto  all  surfaces  and  thus  to  provide  two-dimensional  gaseous  and  hquid  contact 
between  the  boron  nitride  platelets  and  the  copper  wires  inserted  mto  the  boron  nitride 
powder  for  this  purpose.  Enough  hehum  was  admitted  into  the  cell  to  coat  all  surfaces 
with  at  least  ten  layers.  However,  while  the  hehum  was  fi-ee  to  adsorb  onto  the  surface  of 
the  boron  nitride  before  the  introduction  of  the  hydrogen  sample,  once  gaseous  hydrogen 
was  admitted  to  the  cell  all  hehum  adsorbed  onto  the  boron  nitride  surface  was  assumed  to 
be  replaced  by  hydrogen  due  to  the  much  larger  hydrogen-boron  nitride  binding  energy. 
Heat  capacity  experiments  concerning  the  translational  phases  of  helium  adsorbed  onto  a 


^  Unfortunately,  a  leak  from  the  sample  can  to  the  vacuum  can  has  the  same  NMR  signature  which  lead 
to  ample  confusion  until  the  sample  can  was  made  leak  tight. 


layer  of  hydrogen  which  was  itself  adsorbed  onto  graphite  verifies  this  assumption  were 
later  carried  out  at  the  University  of  Washington  [Ebe  94]. 

Adsorption  isotherms  of  heUum  on  boron  nitride  at  14-20  K  were  performed  in  the 
isotherm  apparatus  to  quaUtatively  measure  the  hehum-boron  nitride  adsorption  energy  at 
these  higher  temperatures.  The  measured  pressure  dependence  upon  the  number  of  atoms 
in  the  system  varied  neghgibly  fi-om  the  ideal  gas  law.  In  other  words,  except  at  the 
lowest  coverages  and  pressures,  the  helium  did  not  adsorb  onto  the  bare  boron  nitride.  In 
contrast,  adsorption  isotherms  at  4  K  did  show  monolayer  adsorption.  At  lower 
temperatures  the  hehum  gas  is  assumed  to  adsorb  onto  both  the  hydrogen  surface  (which 
itself  was  adsorbed  onto  the  boron  nitride)  and  the  brass  walls  of  the  sample  can.  Thus, 
the  hydrogen  was  effectively  sandwiched  between  the  boron  nitride  and  the  adsorbed 
heUum.  The  wetting  properties  of  hcHum  on  weak  substrates  (e.g.  graphite,  cesium,  and 
hydrogen)  is  currently  the  subject  of  many  experimental  endeavors 
[Cro  94,  Pai  94,  Nac  94,  Ebe  94]. 

To  test  the  efiFect  of  the  amount  of  hehum  on  thermal  contact,  differing  amounts  of 
hehum  were  used  from  ex-perimental  run  to  run.  Furthermore,  additional  helium  was 
added  to  the  sample  can  during  the  middle  of  some  nms  without  any  measurable  change  in 
NMR  hne  shape  or  amplitude.  Both  the  shape  and  the  amplitude  of  the  NMR  absorption 
signal  did  depend  upon  temperature  for  all  temperatures  studied  and  no  hysteretic  effects 
were  measured.  Thus,  lacking  any  other  evndence  to  the  contrary,  the  sample  was 
assumed  to  be  in  thermodynamic  equilibrium  with  the  thermometer  momited  on  the  mi.xing 
chamber  of  the  dilution  refrigerator  even  thougli  the  magnetization  followed  much  less 
than  a  Curie-like  dependence  upon  temperature. 

Recently,  Vilches,  et  oi,  have  measured'  the  behavior  of  helium  adsorbed  onto 
hydrogen  adsorbed  onto  graphite  using  heat  capacity  measurements  in  conjunction  with 


^  Private  communications. 
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adsorption  isotherms.  They  observed  no  interlayer  mbdng  for  this  system.  In  other 
words,  the  hydrogen  is  so  tightly  bound  to  the  graphite  surface  relative  to  the  helium  at 
temperatures  below  4  K,  that  the  heUum  can  not  displace  it  and  adsorb  directly  onto  the 
substrate.  However,  they  found  the  opposite  result  when  adsorbing  two  near  pairs  of 
hydrogen  isotopes,  H2/HD,  or  HD/D2,  onto  the  graphite  substrate.  In  this  case,  the 
heavier  isotope  preferentially  adsorbed  onto  the  graphite  surface  at  the  expense  of  the 
Ughter  hydrogen  in  a  manner  that  depended  both  upon  the  temperature  and  the  dosing 
level  [Min  93]. 

Since  heat  leaks  to  the  NMR  cell  needed  to  be  minimized  to  achieve  the  lowest 
temperatures  with  the  dilution  refrigerator,  the  fill  line  for  the  NMR  sample  cell  was  much 
smaller  than  that  for  the  adsorption  isotherm  apparatus.  Thus,  precise  adsorption 
isotherms  measured  on  the  NMR  apparatus  were  impossible  due  to  the  long  pressure 
equihbrium  time,  approximately  6  hours.  Instead,  we  measured  the  number  of  atoms  of 
Argon  required  to  attain  a  final  pressure  of  34  Torr,  conditions  which  result  in  a  kink  in 
the  isotherm.  This  characteristic  kink,  measured  in  our  adsorption  isotherm  apparatus  and 
discussed  by  Migone,  has  been  identified  as  a  two-dimensional  fluid-soHd  transition 
[Mig  93],  occurs  at  1.06  monolayers,  and  was  used  to  determine  the  density  of  atoms  on 
the  boron  nitride  surface.  From  this  measurement  and  knowledge  of  the  associated 
volumes  utihzed,  the  number  of  atoms  per  layer  were  determined  to  within  5  percent. 


4.2  Data  Processing 

Data  processing  involved  averaging  files  recorded  under  the  same  conditions  into 
temperature  and  time  "blocks."  Tlie  data  analysis  programs  allowed  for  a  broad  range  in 
the  size  of  these  data  blocks.  Tlie  finest  grid  of  data  blocks  utilized  in  the  analysis  was 
divided   into    100   milliKelvin   and   six  hour  increments  while  the   coarsest   grid   was 
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approximately  1  K  and  one  day.  An  adequate  signal  to  noise  ratio  for  each  block  was  only 
achieved  when  10-30  or  more  sweeps  (dependmg  upon  temperature  and  concentration) 
were  averaged  together.  The  data  was  simuhaneously  processed  into  plots  showing  the 
magnetization  and  differential  amphtude  versus  temperature  and  time.  However,  these 
graphs  turned  out  to  be  less  insightflil  than  the  averaged  hne  shapes  themselves  due  to  the 
large  scatter  in  the  raw  data  and  the  fact  that  the  NMR  line  shape  is  the  most  flmdamental 
measure  of  ordering. 

The  blocks  of  averaged  CW  NMR  line  shape  derivative  amphtude  could  be 
displayed  either  as  functions  of  time  or  temperature.  Four  sets  of  data  arranged  by 
temperature  range  are  included  as  appendices.  Each  experimental  run  lasted  5-25  days, 
used  four  to  eiglit  layers,  and  ranged  in  temperature  from  0. 1-7K.  The  final  analysis 
involved  visually  looking  for  trends  in  the  CW  NMR  hne  shape  derivative  amplitude  as  a 
flmction  of  temperature,  time,  and  the  amount  of  H2  adsorbed  onto  the  surface. 

Almost  two  dozen  data  collection  runs  were  attempted,  one  dozen  of  which 
yielded  usable  data.  Tlie  remaining  runs  developed  obvious  problems  during  the  course  to 
the  run.  These  problems  were  either  with  the  set-up  (e.g.,  hydrogen  frozen  in  the  capillary 
fill  line  instead  of  on  the  boron  nitride),  electrical  (e.g.,  an  "open"  in  the  NMR  tank 
circuit),  or  mechanical  (e.g.,  a  leak  in  the  dilution  refiigerator).  All  data  obtained  were 
highly  reproducible  under  varying  NMR  conditions,  including  RF  power,  modulation 
frequency,  and  lock-in  settings.  The  results  may  be  summarized  using  just  four  sets  of 
data:  three  for  8  layers  obtained  in  June,  October,  and  November  1993  and  one  for  4 
layers  obtained  in  February  1994.  More  information  is  available  to  the  interested  reader  in 
the  appendices. 
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4.3  H?  CW  NMR  line  shapes 

During  the  year-long  investigation  of  ahnost  one  dozen  samples  of  hydrogen  adsorbed 
onto  BN,  certain  patterns  became  evident.  The  resuUs,  as  described  below,  were  found  to 
be  highly  reproducible  across  samples.  Thus,  in  order  to  simplify  the  matters,  line  shapes 
from  only  four  of  the  dozen  samples  will  be  presented  below.  While  the  NMR  hne  shape 
is  the  fundamental  object  of  interest,  the  resonance  derivative  was  actually  measured.  As 
mentioned  above,  the  resonance  derivative  is  also  more  sensitive  to  small  changes  in  the 
hne  shape  than  the  line  shape  itself  Thus,  aU  figures  present  the  CW  NMR  derivative 
amphtude  as  a  function  of  frequency.  The  integral  of  the  presented  data  is  the  CW  NMR 
line  shape.  No  line  shapes  were  scaled  for  presentation.  The  data  from  four  samples 
measured  over  a  period  of  eight  months  is  presented. 

Figure  4-2  illustrates  the  onset  of  orientational  ordering  using  two  derivative  hue 
shapes  from  an  eight  layer  sample  measured  in  November  1993  as  an  example.  Two 
derivative  signals  for  two  distinct  temperatures,  2.38  and  1.06  K,  are  presented  and  a 
strong  central  peak  is  evident.  The  ordering  is  manifested  by  the  formation  of  two  peaks 
or  "wings,"  spaced  appro.ximately  150  kHz  apart  as  marked  by  the  arrows  in  the  figure. 
Note  that  the  signals  are  not  scaled.  The  magnetization  and,  thus,  the  signal  amplitude 
increase  with  decreasing  temperature  but  at  a  rate  less  than  that  predicted  by  Curie's  Law. 
AU  other  samples  exhibited  similar  behavior. 

Figure  4-3  presents  resonance  derivatives  from  a  sample  with  8  layers  of  hydrogen 
adsorbed  onto  boron  nitride  for  a  temperature  of  approximately  0.23  K.  Tliis  data  was 
recorded  in  June  1993.  Tliree  resonance  derivatives,  all  of  which  have  been  smoothed  for 
clarity,  are  shown.  The  largest  derivative  (long-dash)  was  recorded  2  days  after  the 
introduction  of  the  hydrogen  while  the  innermost  derivative  (short-dash)  was  recorded 
15  days  later.  As  the  ortho-hydrogen  converts  to  para-hydrogen  both  the  magnetization 
and  the  orientational  ordering  decay.  Tlie  remaining  signal  amplitude  represents  the 
background  signal.    Although  this  background  signal  is  from  hydrogen,  it  does  not  arise 
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Figure  4-2.  Temperature  dependence  of  the  CW  NMR  line  shape  derivative  for  sample 
BN1193.  Note  both  the  change  of  signal  ampUtude  as  well  as  the  onset  of  the  formation 
of  two  outer  peaks,  or  "wings,"  near  1  K,  as  marked  by  the  arrows. 
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Figure  4-3.  Time  dependence  of  the  resonance  derivative  for  eight  adsorbed  layers 
(sample  BN693)  at  0.2  K.  The  solid  line  is  the  difference  between  the  early-time  (long 
dash)  and  late  time  (short  dash)  data. 
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Figure  4-4.  Time  dependence  of  the  resonance  derivative  for  sample  BN693  at  0.6  K. 
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from  molecular  hydrogen  since  it  does  not  exhibit  a  time  dependence  on  a  time  scale  of  up 
to  two  weeks.  It  is  surmised  that  atomic  hydrogen  permeated  into  the  copper  coil  or 
substituted  for  fluorine  in  the  tetrafluouroethylene  (Teflon)  sample  holder.  The 
background  signal  includes  an  artifact  from  about  50  to  100  kHz.  This  is  probably  due  to 
another  source  of  hydrogen,  perhaps  in  the  Teflon  (used  as  a  dielectric  spacer)  in  the  end 
of  the  coaxial  cable,  which  is  at  a  shghtly  higher  magnetic  field  and  a  comparable  RF 
intensity.  However,  neither  the  background  nor  the  artifact  manifest  a  strong  temperature 
dependence  and  could  be  regarded  as  a  constant  which  could  be  subtracted  from  all  hne 
shapes  as  has  been  done  in  Figure  4-3.  In  practice,  however,  the  data  was  collected  until 
the  molecular  hydrogen  signal  decayed  below  the  background  signal  thereby  yielding  a 
measurement  of  the  background  signal  for  each  nm.  Fortunately,  the  lower  frequency 
data  shows  no  such  artifact  and,  barring  any  unexpected  asymmetries  in  the  NMR  line 
shape,  data  analysis  can  be  performed  on  tliis  side  of  the  data  alone. 

The  sohd  line  of  Figure  4-3  represents  the  difference  between  the  early  and  late 
time  data.  The  shape  of  the  NMR  response  is  that  of  a  combination  of  two  line  shapes:  a 
wide  Pake  Doublet  (approximately  150  kHz)  and  a  narrower  Gaussian  or  Lorentzian 
central  line  (approximately  50  kHz).  In  the  figure,  arrows  mark  each  set.  Note  that  tlie 
central  to  wing  peak  height  ratio  is  approximately  2: 1  for  this  eight  layer  system  at 
200  mK 

Below  1  K,  the  wing  height  of  the  NMR  derivative  hne  shape  was  found  to 
strongly  depend  on  temperature.  Figure  4-4  presents  resonance  derivatives  from  the  same 
sample  as  above  (BN693)  but  at  a  different  temperature.  Again  four  peaks  are  resolvable 
m  the  data.  However,  for  this  0.6  K  data,  the  wings  have  diminished  greatly  in 
comparison  to  the  0.2  K  data  of  Figure  4-2. 

Figure  4-5  presents  resonance  derivatives  measured  in  February  1994  on  a  four 
layer  sample.  As  above,  the  solid  line  in  the  center  of  the  figure  is  the  difference  between 
the  early  and  late  time  data.  As  with  the  eight  layer  sample,  the  Hne  shape  exhibits  two 
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Figure  4-5.  Time  dependence  of  the  resonance  derivative  signal  for  four  adsorbed  layers 
(sample  BN294B).  The  solid  hne  is  the  difference  between  the  early-time  (short  dash)  and 
late  time  (long  dash)  data.  Note  the  four  peaks  of  nearly  equal  amphtude  in  the  difference 
data. 
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Figure  4-6.  Time  dependence  of  the  resonance  derivative  for  eight  layer  sample  BN1093A 
near  the  transition  temperature  (top).  Note  that  the  wing  amplitude  diminishes  while  the 
central  peak  amplitude  remains  constant  at  lower  temperatures  (bottom). 
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pairs  of  peaks,  inner  (gaussian)  and  outer  (Pake  Doublet).  Note,  however,  that  the  ratio 
of  the  inner  to  outer  peak  heights  has  changed  dramatically,  to  approximately  1:1  for  this 
four  layer  sample. 

Figure  4-6  more  precisely  illustrates  the  role  of  the  ortho-hydrogen  concentration 
upon  the  shape  of  the  resonance  derivative.  The  data  is  from  an  eight  layer  sample 
measured  m  October  1993.  For  simpUcity  of  presentation,  only  the  low  frequency  data  is 
displayed.  In  general,  all  measured  Imes  were  found  to  be  symmetric.  In  the  absence  of 
orientational  ordering,  the  NMR  line  shape  should  decay  uniformly  as  the  number  of 
measured  spins  decreases.  In  sharp  contrast,  the  decay  of  the  measured  amphtude  is 
frequency  dependent  as  illustrated  m  Figure  4-6  (bottom).  The  amphtude  of  the  wings 
decreased  faster  than  that  of  the  central  peak.  The  temperature  of  the  onset  of 
orientational  ordering  can  also  be  seen  to  decrease.  Figure  4-6  (top)  displays  the 
resonance  derivative  amplitude  for  four  days  over  a  nine  day  span.  After  15  hours  the 
wing  height  at  1  K  is  quite  evident.  By  132  hours,  the  wing  has  all  but  disappeared  at  1  K 
(top)  and  even  at  0.6  K  (bottom). 

Figure  4-7  presents  two  resonance  derivatives  from  the  experimental  run  of 
November  1993  (BN1193)  which  flirther  illustrates  the  role  of  ortho-hydrogen 
concentration.  Note  that  over  the  two  day  period  exliibited  in  this  figure  the  central  peak 
actually  grows  at  the  expense  of  the  wing  height  for  a  relatively  fixed  temperature!  Tlius, 
spins  that  contributed  to  the  amplitude  of  the  wings  after  33  hours  contributed  to  the 
central  peak  at  76  hours.  During  this  same  time  period,  the  ortho-hydrogen  concentration 
diminished  from  approximately  65  percent  to  45  percent  and  the  overall  magnetization 
(the  second  integral  of  the  resonance  derivative)  decreased. 

To  further  illustrate  the  evolution  of  the  shape  of  the  resonance.  Figure  4-8 
presents  four  resonance  derivatives  from  a  temperature  and  time  sweep  of  sample 
BN1193.  Tliis  figure  includes  one  of  the  resonance  derivatives  (33  hours,  0.63  K)  from 
Figure  4-7.  The  data  was  obtained  by  allowing  the  ortho-hydrogen  concentration  to  decay 
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Figure  4-7.  Time  dependence  of  the  resonance  derivative  for  eight  layer  sample  BN1193. 
Note  that  the  wing  ampUtude  diminishes  while  the  central  peak  amphtude  mcreases. 


105 


I    I    I    I    I    I    I    i    I    I    I    I    I    I    I    i    i    1    I    I    I    I    I    I    I 

-150  -100  -50 


0 


FREQUENCY  (kHz) 


Figure  4-8.  Time  and  temperaUire  dependence  of  eight-layer  sample  BNl  193. 
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while  decreasing  the  sample  temperature  with  the  dUution  refrigerator.  The  data  is  neither 
scaled  nor  oSset  in  the  figure.  Note  that  the  continuous  evolution  m  the  shape  of  the 
resonance  derivative  is  the  resuh  of  two  effects:  temperature  dependent  magnetization 
(Curie's  Law)  and  ortho-hydrogen  concentration  dependent  orientational  ordering. At  4  K, 
the  magnetization  is  relatively  small  and  there  exists  no  orientational  ordering.  At  2.4  K 
only  the  magnetization  has  increased.  By  1  K,  the  wing  has  popped  up  and  the 
magnetization  continues  to  grow  in  a  Curie-like  manner.  Between  1  K  and  0.6  K,  the 
contribution  from  the  wing  has  iucreased  greatly.  However,  the  magnetization  fails  to 
increase  proportionally.  The  magnetization  was  found  to  be  rather  temperature 
independent  below  IK.   It  is  not  known  where  the  lost  magnetization  went. 

To  explore  the  role  of  ortho-hydrogen  concentration  upon  the  NMR  hne  shape, 
Figure  4-9  presents  derivative  amphtudes  for  sample  BN693  for  a  frxed  temperature 
(Tave  =  0.23  K).  In  this  figure,  the  baseUne  of  each  shape  has  been  offset  by  the  amount  of 
time  since  sample  introduction.  Thus,  zero  time  on  the  right  axis  corresponds  to  the  time 
at  which  hydrogen  adsorbed  onto  the  BN  surface.  No  scaling  has  been  performed  upon 
the  line  shapes.  Two  angled,  broken  Unes  with  a  slope  of  25  hours/kHz  have  been  added 
to  this  figure  and  serve  as  a  guide  to  the  eye.  Note  that  the  wings  spread  from  the  center 
of  the  line  shape  as  the  ortho-concentration  decreases.  Additional  figures  presenting  the 
same  behavior  may  be  found  in  the  appendices. 


4.4  Discussion 

The  observation  of  a  less  than  Curie-like  temperature  dependence  of  the 
magnetization  below  1  K  was  of  great  concern.  Tlie  success  of  the  experimental  design 
rested  in  part  upon  the  increase  in  the  signal  to  noise  ratio  below  1  K.  As  presented  in 
Section  3.1,  the  Unear  approximation  to  the  magnetization  should  break  down  below  1  K 
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and  the  exponential  dependence  of  the  magnetization  become  evident.   In  other  words,  the 
magnetization  should  increase  greatly  below  1  K. 

One  obvious  explanation  involves  a  possible  saturation  of  the  NMR  signal.  Below  the 
saturation  point  of  the  system  the  NMR  amphtude  is  linearly  dependent  upon  the  RF 
amphtude  m  a  positive  fashion.  Upon  saturation,  however,  an  increase  in  the  RF 
amphtude  actually  produces  a  decrease  in  the  NMR  amphtude.  Thus,  the  RF  amphtude 
was  varied  significantly  (orders  of  magnitude)  for  test  purposes  and  no  saturation  was 
Observed.  Furthermore,  the  Une  shape  was  not  found  to  be  RF  level  dependent  for  any 
temperature  above  100  mK.  Below  this  temperature,  both  the  RF  fields  and  the  audio 
frequency  magnetic  modulation  produced  considerable  and  instantaneous  heating. 
Another  possible  explanation  mvolves  the  accuracy  of  the  sample  temperature  readings. 
The  thermal  contact  between  the  thermometer  (mounted  on  the  mLxing  chamber)  and  the 
sample  may  not  have  been  sufiBcient  below  a  certain  temperature.  In  tliis  case,  a  physical 
measurement  on  the  sample  in  addition  to  the  thermometer  provides  the  best  indication  of 
the  actual  sample  temperature.  This  was  recognized  early  during  the  experiments  and  the 
shape  as  well  as  the  amplitude  of  the  NMR  signal  itself  was  used  as  such  a  test.  Wliile 
most  NMR  hne  shapes  were  acquired  continuously  as  the  temperature  was  swept,  some 
were  acquired  after  the  temperature  had  been  adjusted  and  allowed  to  equilibrate  for  some 
time.  No  differences  in  the  magnetization  were  observed.  Furthermore,  temperature 
hysterises  effects  were  not  measured  in  any  sample  or  across  samples.  Tlie  resonance 
shape  for  all  8  layer  samples  at  a  given  temperature  and  time  were  indistinguishable 
regardless  of  the  direction  of  the  temperature  sweep.  Thus,  it  is  surmised  that  the 
temperature  dependence  of  the  magnetization  was  real  for  this  system  and  the  dilution 
refiigerator  thermometer  was  in  fact  in  good  thermal  contact  with  the  sample. 


108 


I  I  I  I  I  I 


^ 
U 
rt 


W 
Q 

h 

Oh 

> 
h 

> 

hH 

W 
Q 


I  I  I  I  I  I  M  M  I  I  I  I  I  I  I  M  I  I  M  I  I  I  I. 

TiVE  —  0.23  K 

H2  BN 

8  Layers 


500 


300^ 


M 


0 
100^ 


100 


-50  50  150 

FREQUENCY  (kHz) 


Figure  4-9.   Time  dependence  of  derivative  amplitude  shape  for  sample  BN693.   Broken 
lines  are  guides  to  the  eye  and  have  a  slope  of  25  hours/kHz. 
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It  is  proposed  that  the  sample  is  composed,  with  respect  to  orientational  ordering, 
of  a  two-component  system.  The  layers  near  the  surface  order  below  1  K  under  the 
additional  mfluence  of  the  magnetic  moments  of  the  surface.  However,  the  higher  layers 
remain  disordered  for  even  the  lowest  measured  temperatures  (100  mK).  This  two- 
component  spin  system  is  borne  out  by  two  results:  the  mner-to-outer  peak  height  ratio 
for  the  four  and  eight  layer  data  as  well  as  the  change  in  the  CW  NMR  hne  shape  as  a 
flmction  of  ortho-hydrogen  concentration.  In  this  case,  the  distribution  of  the 
orientational  order  parameter,  as  illustrated  in  Figure  1-7  part  (h),  has  a  gap  between  the 
highly  ordered  state  and  the  disordered  state.  The  first  two  layers  may  become  ordered 
due  to  a  combmation  of  the  adatom-adatom  (i.e.,  quadrupolar-quadrupolar)  interactions 
and  the  adatom-surface  interaction.  This  ordering  is  observed  for  concentrations  of  at 
least  50  percent  ortho-hydrogen  and  for  temperatures  below  approximately  1  K.  The  data 
for  lower  concentration  samples  is  too  noisy  to  be  analyzed  concretely.  The  dependence 
of  the  ordering  on  the  ortho-hydrogen  concentration  impUes  that  the  surface,  at  most, 
participates  in  the  ordering  but  is  not  the  sole  cause  of  the  ordering.  An  example  of  the 
substrate  induced  ordering  includes  NMR  studies  of  deuterium  on  magnesium  oxide 
[Jeo  93]. 

The  temperature  dependence  for  this  two-state  system  of  the  hydrogen  layers  on 
boron  nitride  is  illustrated  in  Figure  4-10.  The  outer  peaks  which  are  visible  below  1  K 
(top  left)  disappear  upon  warming  (top  right)  for  all  measured  concentrations.  Tlie 
associated  order  parameter,  a,  also  undergoes  a  change.  Figure  4-10  (middle)  illustrates 
the  temperature  dependence  of  the  mean  order  parameter.  Note  that  the  ordered  portion 
of  the  line  (the  wings)  changes  more  abruptly  (around  1  K)  than  the  unordered  portion  of 
the  line  (the  center).  However,  the  data  describing  the  order  of  either  phase  transition  is 
incomplete  and  this  phenomenon  deserves  further  investigation.  In  three  dimensions,  the 
orientational  ordering  transition  is  of  first  order  and  accompanied  by  a  change  in  the 
crystal  structure  as  described  above. 
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The  bottom  of  Figure  4-10  illustrates  the  very  interesting  dependence  of  the  degree 
of  orientational  ordering  upon  the  ortho-hydrogen  concentration.  Note  that  contrar>'  to  aU 
previous  orientational  ordering  resuUs,  the  degree  of  ordering  mcreases  for  the  ordered 
portion  of  the  line  even  as  the  total  ortho-hydrogen  concentration  decreases.  There  exists 
at  least  two  possible  causes  for  this  behavior.  First,  the  ortho-hydrogen  may,  in  tune, 
difilise  through  the  sample  and  form  islands  of  high  ortho-hydrogen  concentration.  Thus, 
while  the  total  ortho-hydrogen  concentration  decreases,  these  islands  have  increasing 
ortho  concentrations  for  a  certain  period  of  tune.  Second,  the  surface  may  preferentially 
adsorb  the  ortho-hydrogen  over  para-hydrogen.  Such  adsorption  behavior  is  well  known 
for  some  substrates  and  has  been  used  to  manufacture  samples  of  molecular  hydrogen  with 
enriched  ortho  concentration.  However,  further  studies,  including  monolayer  experiments, 
are  needed  to  determine  the  actual  cause  of  this  behavior. 

It  is  interesting  to  note  that  the  third  and  higher  layers  showed  no  signs  of  ordering 
even  at  the  lowest  temperatures  e.xamined  (approximately  100  mK),  the  highest  ortho- 
concentration  studied  (about  70  percent),  and  the  highest  number  of  layers  utilized  (eight 
layers).  While  this  work  does  not  rule  out  ordering  for  higher  concentration  samples  at 
even  lower  temperatures  (below  0. 1  K),  the  lack  of  ordering  for  the  higlier  layers  is  in 
contradistinction  to  the  behavior  of  bulk  hydrogen  which  orientationally  orders  below 
approximately  1  K  for  60  percent  ortho-hydrogen  samples.  From  this  work  it  may  be 
inferred  that  "true"  two-dimensional  samples,  that  is,  monolayer  samples,  on  isotropic 
surfaces  may  either  order  only  for  the  highest  ortho-hydrogen  concentrations  and  at  the 
lowest  temperatures  or  may  not  orientationally  order  under  any  conditions  due  the  large 
orientational  zero  point  motion  of  hydrogen.  Either  way,  changing  the  dimensionality  of 
the  system  affects  the  ordering  beiiavior  dramatically.  Tlie  ortho-para-hydrogen 
conversion  rate  for  all  runs  was  similar  and  approximately  1/5  of  the  bulk  hydrogen  rate. 

A  comprehensive  study  of  the  orientational  ordering  of  monolayers  of  hydrogen  on 
graplute  was  conducted  some  years  ago  by  Kubik  and  Hardy  at  the  University  of  British 
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Columbia  at  Vancouver  [Kub  83].  The  surface  of  graphite  shares  the  hexagonal  lattice 
structure  of  boron  nitride  but  is  less  corrugated.  The  graphite  surface  is  also  beheved  to 
be  more  isotropic  than  BN.  Kubik  and  Hardy  utUized  "enriched"  samples  of  hydrogen 
(e.g.,  using  samples  with  greater  than  75  percent  ortho-hydrogen)  and  did  not  find  a 
completely  ordered  state  for  their  system  However,  data  analysis  was  compUcated  by  the 
presence  of  a  rather  large  background  signal  which  varied  in  both  shape  and  amplitude  as 
a  flmction  of  temperature. 


CHAPTER  5 
SUMMARY  AND  IMPLICATIONS 

Boron  nitride  has  proven  to  be  a  useful  substrate  for  the  study  of  two-dimensional 
molecular  systems.  BN  is  homogeneous,  clean,  and  very  stable.  While  BN  possesses  the 
same  hexagonal  planar  molecular  structure  as  graphite,  its  electrical  properties  are  quite 
different.  BN  is  an  insulator  while  graphite  is  a  semi-metal.  This  difference  can  be  very 
useful  for  certain  experiments  as  demonstrated  here.  However,  for  volume  sensitive 
experiments,  such  as  NMR,  the  lower  surface-to-volume  ratio  of  BN  is  a  drawback. 

Although  the  adsorption  isotherms  accomplished  most  of  their  goals,  some 
limitations  were  encountered.  On  the  plus  side,  the  sharpness  of  the  steps  manifest  layer 
by  layer  adsorption.  The  adsorption  potential  and  critical  temperatures  have  also  been 
measured.  However,  the  measurements  were  limited  by  two  elements  in  the  design  of  the 
apparatus:  first,  the  capillary  between  the  loading  volume  and  the  sample  cell  was  too 
narrow;  second,  the  omission  of  a  separate  line  heater  in  addition  to  the  sample  cell  heater 
(which  was  used  for  both  purposes).  Together,  these  two  difficulties  conspired  to  limit 
both  the  accuracy  and  the  precision  of  the  measurements  at  the  lowest  pressures.  Both 
problems  have  been  corrected  on  the  newest  generation  of  the  apparatus.  However,  to 
maintain  continuity  in  the  experiment,  this  new  dipstick  could  not  be  pressed  into  service 
until  the  completion  of  the  current  measurements  on  the  BN  sample. 

BN  also  provides  another  surface  to  compare  and  contrast  the  extensive  work 
performed  over  the  last  few  decades  on  other  substrates.  While  work  on  surfaces  has 
been  wide-ranged  and  includes  aerogels,  vycor,  and  zeolites,  the  most  extensive  work  has 
been  performed  on  the  various  forms  of  graphite.   The  results  obtained  herein  build  upon 
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the  previous  results  and  add  to  the  set  of  known  two-dimensional  and  quasi  three- 
dimensional  systems. 

The  most  interesting  result  of  previous  experiments  may  concern  the  search  for 
superfluidity  in  reduced  dimensions.  The  possibility  of  superfluidity  in  supercooled 
hydrogen  has  been  proposed  for  some  time  and  a  number  of  searches  have  been  made. 
The  adsorption  potential  of  BN  had  been  theorized  to  be  less  than  graphite  and  that  this 
reduction  would  translate  into  a  lowering  of  the  critical  and  triple  point  temperatures  of 
the  adsorbed  hydrogen  films.  However,  such  a  further  suppression  of  these  temperatures 
was  not  found.  The  critical  points,  as  measured  by  the  adsorption  isotherms  were  found 
to  be  somewhat  higher  than  the  analogous  systems  on  graphite.  The  corrugation  of  the 
surface  is  surmised  to  play  some  role  here.  However,  these  experiments  represent  the  first 
effort  to  measure  the  critical  temperatures  and  the  method  used  to  determine  these  points 
tends  to  overestimate  the  critical  temperatures.  Thus,  these  points  should  be  checked  with 
another  experimental  technique  such  as  heat  capacity  measurements. 

No  substeps  were  found  even  though  the  resolution  of  the  adsorption  isotherm  was 
adequate  to  uncover  details  at  a  level  corresponding  to  one  percent  of  a  monolayer.  This 
resolution  was  tested  using  a  small  Argon  on  BN  substep.  Thus,  the  triple  point 
temperatures  for  the  hydrogen  isotopes  adsorbed  onto  BN  could  not  be  measured  using 
the  customary  adsorption  isotherm  technique.  The  lack  of  a  significant  substep  could 
indicate  that  the  phase  diagram  of  the  hydrogens  on  BN  is  significantly  different  than  the 
diagrams  of  the  same  adatoms  adsorbed  onto  graphite.  For  this  and  other  reasons  the 
search  for  superfluidity  of  hydrogen  on  BN  remains  open. 

The  role  of  dimensionality  in  orientational  ordering  has  also  been  explored. 
Evidence  of  orientational  ordering  of  ortho-hydrogen  adsorbed  onto  BN  has  been 
presented  for  temperatures  below  approximately  1  K.  This  two-dimensional  ordering  has 
proven  to  be  both  more  subtle  and  significantly  different  than  its  three-dimensional 
counterpart.  For  example,  the  increase  in  ordering  as  a  fijnction  of  time  for  part  of  the  line 
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is  a  new  phenomena.  Although  the  data  did  not  allow  for  a  determination  of  the  nature  of 
the  phase  transition,  the  ordering  was  shown  to  depend  both  upon  the  number  of  adsorbed 
layers  and  the  concentration  of  ortho-hydrogen.  The  boron  nitride  surface  is  proposed  to 
play  only  a  partial  role  in  the  orientational  ordering  of  the  adsorbed  hydrogen 

Thus,  BN  provides  a  surface  where  the  goal  of  studying  ideal  two-dimensional 
systems  is  partially  attained.  This  is  in  sharp  contrast  to  some  systems  (i.e.,  deuterium  on 
MgO)  where  the  surface  was  shown  to  cause  orientational  ordering  at  relatively  high 
temperatures  [Jeo  93].  Upon  improvements  to  the  sensitivity  of  the  NMR  apparatus,  the 
nature  of  the  orientational  phase  transition  and  monolayer  studies  may  be  investigated. 
Such  improvements,  including  the  addition  of  low-temperature  tuning  capacitors,  are 
currently  underway. 

This  work  can  be  extended  and  the  role  of  the  substrate  further 
investigated.  One  or  more  layers  of  an  inert,  rotationally  symmetric  molecule,  such  as 
krypton  or  argon  (neon  has  been  shown  to  mix  with  adsorbed  H2),  may  be  adsorbed  onto 
the  surface  prior  to  the  adsorption  of  hydrogen.  These  adatoms  would  insulate  the 
hydrogen  from  the  corrugation  and  orientational  anisotropy  of  the  BN  surface.  By  using 
prelayer  adatoms  of  various  sizes  and  by  altering  the  number  of  layers,  the  role  of  the 
surface  can  be  extensively  investigated.  Such  work  is  the  subject  of  current  research. 


LIST  OF  REFERENCES 

[Abr  91]  A.  Abragan,  Principles  of  Nuclear  Magnetism,  Oxford  Science  Publications, 
Oxford,  1961. 

[Ahl  64]  G.  Ahlers  and  W.H.  Ortung,  Phys  Rev.  133  (1964)  A1642. 

[Alk  92]  M.T.  Alkhafaji  and  A.D.  Migone,  Phys.  Rev.  B,  45  (1992)  5729. 

[And  92a]  P.  Anderson,  N.S.  Sullivan,  B.  Andraka,  J.S.  Xia,  and  E.D.  Adams, 
Phys  Lett  A  172  (1992)  85. 

[And  92b]  P.  Anderson,  N.S.  Sullivan,  B.  Andraka,  J.S.  Xia,  and  E.D.  Adams, 
J.  Low  Temp  Phys  89  (1992)  715. 

[Ave  92]  O.  Avenel,  J.  Xu,  J.S.  Xia,  M-F  Xu,  B.  Andraka,  T.  Lang,  P.L.  Moyland, 

W.  Ni,  P.J.  Signore,  CM.  van  Woerkens,  E.D.  Adams,  G.G.  Ihas,  M.W.  Meisel, 
S.E.  Nagler,  N.S.  Sullivan,  Y.  Takano,  D.R.  Talham,  T.  Goto,  N.  Fujiwara,  Phys. 
Rev.  8  46(1992)  1993. 

[Bar  85]  Randall  F.  Barron,  Cryogenic  Systems,  second  ed.  Oxford  University  Press, 
New  York,  1985. 

[Cha  84]  M.H.W.  Chan,  AD.  Migone,  K.D.  Miner,  and  Z.R.  Li,  Phys.  Rev.  B 
30(1984)2681. 

[Cha  85]  F.A.B.  Chaves,  M.E.B.P.  Cortez,  RE.  Rapp,  and  E.  Learner,  Surf.  Sci. 
150(1985)80. 

[Chu  77]  T.T.  Chung  and  J.G.  Dash,  Surf.  Sci.  66  (1977)  559. 

[Cow  84]  B.P.  Cowan  and  A.J.  Kent,  J.  Low  Temp.  Phys.  57  (1984)  589. 

[Cra  94]  T  Crane  and  B.  Cowan,  Physica  B  194-6  (1994)  633. 

[Cro  94]  PA.  Crowell  and  J.D.  Reppy,  Physica  B,  194-196  (1994)  621. 

[CRC  91]  D.R.  Lide,  Ed.,  CRC  Handbook  of  Chemistry  and  Physics,  72nd  edition,  CRC 
Press,  Boca  Raton,  1991. 


116 


117 

[Das  75]  J.G.  Dash,  Films  on  Solid  Surfaces,  Academic  Press,  New  York,  1975 

[Die  82]  R.D.  Diehl,  M.F.  Toney,  ans  S.C.  Fain,  Jr.,  Phys.  Rev.  Lett.  48  (1982)  177. 

[Ebe  94a]  P.S.  Ebey,  Y.M.  Liu,  and  O.  Vilches,  Physica  B,  194-196  (1994)  621. 

[Ebe  94b]  P.S.  Ebey,  Y-M  Liu,  and  O.E.  Vilches,  Physica  B  ,  194-196  (1994)  635. 

[Ebe  94c]  P.S.  Ebey,  private  communications. 

[Eck  86]  R.E.  Ecke,  J.  Ma,  D.  Migone,  and  T.S.  Sullivan,  Phys.  Rev.  B  33  (1986)  1746. 

[Pre  85]  H.  Freimuth  and  H.  Wiechert,  Surf.  Sci.  162  (1985)  432. 

[Fre  90]  H.  Freimuth,  H.  Wiechert,  H.P.  Schildberg,  and  H.J.  Lauter,  Phys.  Rev.  B  42 
(1990)  587. 

[Fuk  81]  E.  Fukushima  and  S.  Roeder,  Experimental  Pulse  NMR:  a  Nuts  and  Bolts 
Approach,  Addison-Wessley,  Reading  MA,  1981. 

[Gin  72]  V.L.  Ginzburg  and  A.  A.  Sobyanin,  JETP  Lett.  15  (1972)  242. 

[Goo  93]  D.  Goodstein,  P.  Day,  M.  LaMadrid,  and  M.  Lysek,  Mat.  Res.  Soc.  Symp. 
Proc.  280(1993)319. 

[Got  90]  J.M.  Gottlieb  and  L.W.  Baich,  Phys.  Rev.  B  41  (1990)  7195. 

[Ham  83]  J.J.  Hamilton  and  D.L.  Goodstein,  Phys.  Rev.  B  28  (1983)  3838. 

[Har  79]  A.B.  Harris  and  A.J.  Berlinsky,  Can.  J.  Phys.  48  (1979)  1852. 

[Jeo  93]  E.K.  Jeong,  B.  Ouyang,  RE.  Norberg,  PA.  Fedders,  and  M.S.  Conradi, 
Phys.  Rev.  Lett.  69  (1993)  2983. 

[Kri  84]  J.  Krim,  Ph.D.  dissertation.  University  of  Washington,  Seattle  (1984) 

[Kub  78]  P.R.  Kubik  and  W.N.  Hardy,  Phys.  Rev.  Lett.  41  (1978)  257. 

[Kub  85]  P.R.  Kubik,  W.N.  Hardy,  and  H.  Glattli,  Can.  J.  Phys.  63  (1985)  614. 

[Lah  82]  J. P.  Laheurte,  J.C.  Noiray,  M.  Obadia,  and  J. P.  Romagnan, 
Surface  Science  122  (1982)  330. 

[Lan  69]  L.D.  Landau  and  EM.  Lifshitz,  Statisical  Physics  (Addison-Wesley, 
Reading,  MA,  1969),  Section  59. 


118 


[Lar  79]  Y.  Larher,  Mol.  Phys.  38  (1979)  789. 

[Laz  37]  B.  Lazarew  and  L.  Schubnikow,  Phys.  Z.  Sowjet.  11  (1937)  445. 

[Liu  93]  Y.M.  Liu,  Ph.D.  dissertation,  University  of  Washington,  Seattle  (1993). 

[Liu  92a]  F.C.  Liu,  Ph.D.  dissertation.  University  of  Washington,  Seattle  (1992). 

[Liu  92b]  F.C.  Liu,  Y.M.  Liu,  and  O.E.  Vilches,  J.  Low  Temp.  Phys.  89  (1992)  649. 

[Lys  92]  M.J.  Lysek,  M.  LaMadrid,  P.  Day,  and  D.  Goodstein,  Langmuir  8  (1992)  898. 

[Ma  89]  Jian  Ma,  Ph.D.  dissertation.  University  of  Washington,  Seattle  (1989). 

[Mar  83]  H.J.  Maris,  G.M.  Seidel,  and  T.E.  Huber,  J.  Low  Temp.  Phys.  51  (1983)  477. 

[Mar  93]  D.  Marx,  S.  Sengupta,  and  P.  Nielaba,  J.  Chem.  Phys.  99  (1993)  6031 
and  references  therein. 

[Men  59]  K.  Mendelssohn,  Progress  in  Cryogenics,  Vol.  1,  Academic  Press, 
New  York,  1959. 

[Mig  93]  A.D.  Migone,  M.T.  AJkhafaji,  G.  Vidali,  and  M.  Karimi, 
Phys.  Rev.  B  47  (1993)  6685. 

[Mig  94a]  A.  Migone,  P.  Shrestha,  M  AJkhafaji,  M.  Lukowitz,  and  G.  Yang,  preprint. 

[Mig  94b]  A.  Migone,  private  communications. 

[Mil  66]  R.L.  Mills,  A.F.shuch,  and  D.A.  Depatie,  Phys.  Rev.  Lett.  17  (1966)  1131. 

[Mis  92]  G.  Mistura,  H.C.  Lee,a  nd  M.H.W.  Chan,  J.  Low  Temp.  Phys. 
89,  633  (1992) 

[Mot  85]  F.C.  Motteler  and  J.G.  Dash,  Phys,  Rev.  B  31  (1985)  346. 

[Muc  66]  K.F.  Mucker,  S.  Talhouk,  P.M.  Harris,  and  R.A.  Erikson.,  Phys.  Rev.  Lett. 
16(1966)  176. 

[Mur  93]  M.  Maruyama,  M.  Bienfait,  F.C.  Liu,  Y.M.  Liu,  O.E.  Vilches,  F.  Rieutord, 
Surf.  Sci.  283(1993)333. 

[Nac  94]  P.J.  Nacher,  B.  Demolder,  and  J.  Dupont-Roc,  Physica  B  194-196  (1994)  975. 


119 

[Pai  94]  C.G.  Paine  and  G.M.  Seidel,  Physica  B  194-196  (1994)  969. 

[Phi  93]  J.M.  Phillips,  Q.M.  Zhang,  and  J.Z.  Larese,  Phys  Rev  Let,  71  (1993)  2971. 

[Pip  83]  J.  Piper,  J.A.  Morrison,  C.  Peters,  and  Y.  Ozaki,  J.  Chem.  Soc.  Far. 
Trans.  79  (1983)  2863. 

[Pob  92]  Frank  Pobell,  Matter  and  Methods  at  Low  Temperatures,  Springer- Verlag, 
Beriin,  1992. 

[Pur  46]  E.MN.  Purcell,  H.C.  Torrey,  and  R.V.  Pound,  Phys.  Rev.  69  (1946)  37. 

[Reg  79]  J.  Regnier,  A.  Thorny,  and  X.  Duval,  J.  of  Coll.  Int.  Sci.  70  (1979)  105. 

[Ric  88]  Robert  C.  Richardson  and  Eric  N.  Smith,  Experimental  Techniques  in 

Condensed  Matter  Physics  at  Low  Temperatures,  Addison- Wesley  Publishing 
Co.,  Redwood  City,  CA,  1988. 

[Rol  72]  R.J.  Rollefson,  Phys.  Rev.  Lett.  29  (1972)  410. 

[Ros  83]  T.F.  Rosenbaum,  S.E.  Nagler,  and  P.M.  Horn,  Phys.  Rev.  Lett. 
50(1983)  1791. 

[Sch  93]  M.  Schick  and  O.E.  Vilches,  Phys.  Rev.  B  48  (1993)  9910. 

[Shr  94]  P.  Shrestha,  M.T.  Alkhafaji,  M.M.  Lukowitz,  G.  Yang,  and  A.D.  Migone, 
Langmuir  10(1994)3244. 

[Sil  80]  I.F.  Silvera,  Rev.  Mod.  Phys.  52  (1980)  393. 

[Sli  90]  C.P.  Slichter,  Principles  of  Magnetic  Resonance,  third  ed.,  Sringer- Verlag, 
Beriin,  1990 

[Spr  94]  D.T.  Sprague,  N.  Alikacem,  P.A.  Sheldon,  and  R.B.  Hallock,  Physica 

B,  194-6(1994)629. 

[Sul  67]  N.S.  Sullivan,  Ph.D.  dissertation,  Harvard  University,  Cambridge,  1967. 

[Sul  72]  N.S.  Sullivan  and  R.V.Pound,  Phys.  Rev.  A,  6  (1972)  1 102. 

[Sul  83]  N.S.  Sullivan  and  J.M.  Vaissiere,  Physical  Review  Letters  51  (1983)  658. 

[Suz  77]  J.  Suzanne  and  M.  Bienfait,  J.  Phys,  (Paris)  38  (1977)  C4-93. 


120 


[Suz  91]  J.  Suzane,  J.M.  Gay,  M.  Muruyama,  P.  Zeppenfeld,  D.  Degenhardt,  HJ.  Lauter, 

F.  Rieutord,  G.  Goddins,  Excitations  in  Two-Dimensional  and  Three-Dimensional 
Quantum  Fluids,  A.G.F.  Wyatt  and  H.J.  Lauter  (eds.).  Plenum,  New  York,  1991. 

[Tab  92]  P.  Taboreck  and  J.E.  Rutledge,  Phys  Rev.  Lett.  68,  201  (1992) 

[Tel  83]  J.L.  Tell  and  H.J.  Maris,  Phys.Rev.  B  28  (1983)  5122. 

[Tho  70]  A  Thorny  and  X.  Duval,  J.Chim.  Phys.  Physicchim.  Biol.  67  (1970)  1101. 

[Vid  91]  G.  Vidali,  G  Ihm,  H-Y  Kim,  and  M.W.  Cole,  Surf.  Sci.  Rep.  12,  (1991)  133. 

[Vil  91]  O.E.  Vilches,  F.C.  Liu,  D.L.  Kingsbury,  J.  Ma,  M.  Bienfait,  J.  Suzanne,  J.M. 

Gray,  M.  Maruyama,  P.  Zeppenfeld,  D.  Degenhardt,  H.J.  Lauter,  F.  Rieutord,  and 

G.  Coddens,  Excitations  in  Two-Dimensional  and  Three-Dimensional  Quantum 
Fluids,  A.G.F.  Wyatt  and  H.J.  Lauter  (eds),  Plenum,  New  York,  1991. 

[Vil  92]  O.E.  Vilches,  J.  Low  Temp.  Phys.  89  (1992)  267 

[Vil  94]  O.E.  Vilches,  Y-M.  Liu,  P.S.  Ebey,  and  F-C  Liu,  Physica  B 
194-196(1994)665. 

[Vol  89]  U.G  Volkmann  and  K.  Norr,  Surf.  Sci.  221  (1989)  379. 

[Whi  79]  Guy  K.  White,  Experimental  Techniques  in  Low-Temperature  Physics,  third 
ed.  Clarendon  Press,  Oxford,  1979. 

[Wie  91a]  H.  Wiechert,  Physica  B  169  (1991)  144. 

[Wie  91b]  H.  Wiechert,  m  Excitations  in  Two-Dimensional  and  Three-Dimensional 
Quantum  Fluids,  A.G.F.  Wyatt  and  H.J.  Lauter  ed..  Plenum  Press, 
New  York,  1991. 

[Wliu  93]  W.  Liu  and  S.C.  Fain,  Jr.,  Phys.  Rev.B  47  (1993)  15965. 

[You  90]  H.S.  Youn  and  G.B.  Hess,  Phys.  Rev.  Lett.  64,  918  (1990) 

[Zim  92]  G.Zimmerli,  G.  Mistura,  and  M.H.W.  Chan,  Phys.  Rev.  Lett.  68,  60  (1992). 


APPENDICES 

These  appendices  are  included  to  document  the  CW  NMR  derivative  line  shapes  as 
a  function  of  both  time  and  temperature.  Although  at  first  glance,  presenting  the  data  in 
two  separate  formats  may  seem  redundant,  it  is,  in  fact,  necessary  in  order  to  deconvolute 
the  temperature  dependence  of  the  line  shape  from  the  time  dependence  (caused  by  the 
conversion  of  ortho-hydrogen  to  para-hydrogen).  The  raw  data  consisted  of  many 
(thousands)  of  low  signal  to  noise  line  shapes  including  some  line  shapes  with  obvious 
problems.  In  the  process  of  reducing  the  raw  data  into  blocks  of  average  temperatures 
and  times,  both  the  noisiest  and  most  problematic  raw  data  (approximately  ten  percent) 
was  discarded.  Even  after  this  data  reduction,  it  is  practical  to  present  only  a  summary  of 
the  more  definitive  data.  Thus,  this  appendix  includes  only  data  near  the  orientational 
ordering  transition. 

There  exits  three  purposes  for  this  appendix:  first,  to  detail  the  conditions  under 
which  the  data  were  obtained;  second,  to  provide  the  interested  reader  with  more  explicit 
data  and,  third,  to  preserve  the  data  in  a  compact  and  lasting  format.  As  mentioned  in  the 
text,  almost  two  dozen  experimental  runs  were  attempted  over  a  two  year  period  which 
yielded  about  one  dozen  sets  of  NMR  measurements  of  hydrogen  on  boron  nitride. 
However,  most  of  these  runs  had  experimental  difficulties  which  limited  the  usefulness  of 
the  data.  Since  complete  knowledge  of  the  line  shape  is  only  possible  for  a  broad 
distribution  of  sample  temperatures  and  ortho-hydrogen  concentrations,  multiple 
overlapping  data  sets  must  be  analyzed.  Data  from  four  samples  taken  over  the  most 
productive  eight  month  period  is  presented  in  the  following  sections. 
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The  following  graphs  display  the  most  interesting  line  shapes  for  each  sample  in 
two  formats:  first,  as  a  function  of  time,  that  is,  the  baseline  for  each  line  shape  is  shifted 
by  the  average  amount  of  time  expired  since  sample  introduction;  second,  as  a  function  of 
temperature,  that  is,  the  baseline  for  each  line  shape  is  shifted  by  the  average  sample 
temperature.  The  tables  list  the  time  in  hours  from  the  introduction  of  the  sample,  the 
average  temperature  of  the  sample,  and  the  number  of  CW  NMR  line  shapes  averaged  in 
each  time  and  temperature  block. 


APPENDIX  A 


DETAILS  FROM  SAMPLE  BN693 


This  appendix  presents  selected  data  from  eight  layer  sample  BN693.  This  sample 
was  adsorbed  onto  the  BN  at  7AM,  10  June  1993.  The  discontinuities  near  the  ends  of 
the  line  shape  are  a  result  of  the  drifting  main  magnetic  field.  Spurious  signal  pick-up,  as 
discussed  in  the  text,  limited  the  magnetic  field  sweep  to  less  than  about  300  kHz.  Upon 
aligning  the  CW  NMR  line  shapes  with  respect  to  each  other,  the  lack  of  data  for  the 
wings  of  one  side  of  some  of  the  line  shapes  becomes  evident.  After  subtracting  the  last 
measured  line  shape  from  the  previous  ones,  the  discontinuities  appears  on  both  sides  of 
the  line  shape  for  some  of  the  data  sets. 
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TIME( 

Tic 

3urs) 

TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

42.2 

56.5 

115.3 

130.0 

150.1 

215.7 

223.6 

1 

40.9 

57.3 

91.5 

103.5 

138.9 

152.6 

208.9 

231.3 

2 

68.7 

82.1 

132.4 

156.7 

205.4 

234.2 

3 

73.0 

133.1 

165.5 

202.0 

236.3 

4 

132.3 

159.8 

180.9 

198.5 

239.2 

5 

160.7 

178.3 

197.0 

6 

7 

196.7 

8 

194.8 

196.1 

9 

192.6 

10 

191.3 

11 

190.8 

12 

190.3 

13 

188.3 

AVE 

41.5 

63.8 

86.80 

109.4 

133.3 

157.6 

188.4 

202.6 

232.9 

TEMP 

TIM 

E  BLOCK 

BLOCK 

9 

10 

11 

12 

13 

14 

15 

16 

17 

0 

382.5 

392.4 

424.2 

1 

371.5 

396.5 

422.3 

2 

370.8 

424.3 

3 

367.6 

419.7 

4 

249.4 

277.7 

326.2 

347.0 

364.7 

424.7 

5 

249.8 

284.0 

350.7 

364.4 

432.6 

436.0 

6 

287.5 

323.6 

341.1 

364.2 

7 

322.9 

361.1 

8 

360.1 

9 

305.4 

357.5 

10 

287.8 

303.0 

357.0 

11 

301.2 

356.3 

12 

298,5 

335.1 

355.9 

13 

289.7 

305.4 

323.2 

350.7 

AVE 

249.6 

285.4 

302.7 

326.3 

353.7 

369.4 

394.4 

424.6 

436.0 

Table  A-1.  Average  time,  in  hours,  from  sample  introduction 
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TEMPERATURE  (K) 


TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

0.25 

0.23 

0.22 

0.24 

0.24 

0.22 

0.22 

1 

0.26 

0.28 

0.41 

0.32 

0.29 

0.34 

0.31 

0.35 

2 

0.64 

0.57 

0.61 

0.63 

0.64 

0.59 

3 

0.84 

0.84 

0.89 

0.85 

0.92 

4 

1.06 

1.13 

1.16 

1.05 

1.16 

1.15 

5 

1.26 

1.31 

1.45 

1.29 

6 

7 

1.85 

8 

2.21 

2.21 

9 

2.30 

10 

2.64 

11 

2.87 

12 

3.14 

13 

7.13 

TEMP 

TIME  BLOCK 

BLOCK 

10 

11 

12 

13 

14 

15 

16 

17 

AV 
E 

0 

0.22 

0.23 

0.22 

0.23 

1 

0.39 

0.29 

0.34 

0.33 

2 

0.67 

0.62 

0.62 

3 

0.84 

0.88 

0.86 

4 

1.20 

1.20 

1.18 

1.25 

1.13 

1.15 

5 

1.26 

1.26 

1.38 

1.39 

1.49 

1.34 

6 

1.62 

1.61 

1.62 

1.60 

1.62 

7 

1.89 

1.96 

1.90 

8 

2.08 

2.17 

9 

2.40 

2.35 

2.35 

10 

2.60 

2.60 

2.73 

2.64 

11 

2.89 

2.90 

2.89 

12 

3.07 

3.20 

3.12 

3.14 

13 

5.22 

3.73 

4.54 

4.07 

4.94 

Table  A-2.  Temperature,  in  Kelvin,  of  the  sample. 


126 


NUMBER 


TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

11 

5 

66 

6 

9 

31 

16 

1 

55 

13 

63 

18 

51 

18 

20 

14 

2 

15 

8 

11 

8 

7 

8 

3 

22 

16 

45 

15 

4 

4 

5 

6 

29 

11 

22 

10 

5 

1 

8 

1 

2 

6 

7 

1 

8 

9 

1 

9 

8 

10 

2 

11 

2 

12 

2 

13 

11 

TEMP 

TIMEB] 

LOCK 

BLOCK 

10 

11 

12 

13 

14 

15 

16 

17 

0 

42 

8 

23 

1 

3 

18 

16 

2 

2 

7 

3 

13 

18 

4 

53 

22 

44 

1 

15 

5 

4 

1 

1 

11 

1 

6 

1 

1 

1 

1 

7 

1 

8 

8 

15 

9 

2 

3 

10 

1 

13 

1 

11 

10 

3 

12 

2 

1 

3 

13 

3 

24 

7 

8 

Table  A-3.  Number  of  averaged  line  shapes. 
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Figure  A-1.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  0.23  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  top  figure  shows  the  measured  averaged  line  shapes 
while  the  bottom  figure  displays  the  same  with  the  background  subtracted.  The  dashed 
line     is    a    guide    to     the     eye     and     has     a     slope     of    25     hours     per    kHz. 
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Figure  A-2.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  0.33  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  top  figure  shows  the  measured  averaged  line  shapes 
while  the  bottom  figure  displays  the  same  with  the  background  subtracted.  The  dashed 
line     is     a    guide    to    the     eye     and     has     a     slope     of    25     hours     per    kHz. 
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Figure  A-3.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  0.62  Kelvin.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time 
elapsed  since  formation  of  the  sample.  The  top  figure  shows  the  measured  averaged  line 
shapes  while  the  bottom  figure  displays  the  same  with  the  background  subtracted.  The 
dashed  line  is  a  guide  to  the  eye  and  has  a  slope  of  25  hours  per  kilohertz. 
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Figure  A-4.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  0.86  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  top  figure  shows  the  measured  averaged  line  shapes 
while  the  bottom  figure  displays  the  same  with  the  background  subtracted.  The  dashed 
line  is  a  guide  to  the  eye  and  has  a  slope  of  25  hours  per  kHz. 
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Figure  A-5.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  1.15  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  top  figure  shows  the  measured  averaged  line  shapes 
while  the  bottom  figure  displays  the  same  with  the  background  subtracted. 
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Figure  A-6.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  1.34  and  1.90  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time 
elapsed  since  formation  of  the  sample. 
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Figure  A-7  Dependence  of  CW  hMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  2.35  and  2.64  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time 
elapsed  since  formation  of  the  sample. 
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Figure  A-8.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion 
time  at  2.89  and  3.14  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time 
elapsed  since  formation  of  the  sample. 
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Figure  A-9.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  cx)nversion 
time  at  4.94  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample. 
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Figure  A- 10.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  low 
temperature  data  taken  over  300  hours  apart..  The  baseline  of  each  line  shape  is  offset  by 
the  average  sample  temperature.  The  latter  time  signal  is  essentially  background.  Note 
the  full  disappearance  of  the  peak  at  -75  kHz  and  the  diminished  peak  at  +75  kHz. 
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Figure  A-11.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  high 
temperature  data  taken  almost  200  hours  apart.  The  baseline  of  each  line  shape  is  offset 
by  the  average  sample  temperature.  Note  the  lack  of  outer  peaks  in  both  figures. 


APPENDIX  B 


DETAILS  FROM  SAMPLE  BN1093A 


This  appendix  presents  selected  data  from  sample  BN1093A.  The  sample  was 
adsorbed  onto  the  BN  at  noon,  19  October  1994.  This  experimental  run  attempted  to 
overcome  the  effects  of  an  inadequately  measured  base  line.  Ideally,  one  would  measure 
the  magnetic  response  of  the  sample  for  a  width  at  least  three  times  the  resonance. 
However,  this  was  impossible  for  this  experiment  due,  in  part,  to  the  limited  width  of  the 
magnetic  field  sweep  and,  in  part,  to  the  wide  line  shapes  being  measured.  As  an 
alternative  method  of  measuring  the  baseline,  the  center  of  the  field  sweep  was  moved 
away  from  the  center  of  the  resonance  and  toward  the  baseline.  The  lower  frequency 
baseline  was  chosen  to  avoid  the  more  complicated  background  signal  at  slightly  higher 
frequencies.  Theoretically,  since  the  molecular  hydrogen  line  was  expected  to  be 
completely  symmetrical  with  respect  the  center  of  the  line,  no  information  is  lost.  In 
practice,  however,  only  one  phase  of  the  signal  was  measured.  Much  effort  was  made  to 
keep  the  spectrometer  tuned  to  measure  the  absorption  signal  while  ignoring  the 
dispersion  signal.  Measuring  a  mixture  of  absorption  and  dispersion  creates  a  skewed  line 
shape  and  some  of  the  spin  weight  may  have  been  either  missed  or  erroneously  included. 
Thus,  while  measuring  one  side  of  the  resonance  yields  improved  line  shapes  for  the 
purpose  of  measuring  Pake  doublet  information,  it  is  at  the  expense  of  adequately 
measuring  the  total  magnetization  which  yields  the  conversion  rate.  For  this  reason,  all 
other  samples  were  measured  symmetrically. 
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TIME  (LI 

ours) 

TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

0 

74.7 

91.2 

151.8 

1 

30.5 

37.3 

67.4 

102.2 

144.0 

156.2 

2 

25.0 

41.5 

61.4 

105.2 

106.5 

136.9 

159.1 

3 

20.4 

45.4 

59.5 

108.1 

133.9 

168.7 

4 

14.7 

53.7 

119.3 

131.6 

168.7 

5 

49.2 

117.0 

6 

48.0 

119.1 

7 

10.5 

8 

11.4 

48.2 

126.1 

9 

10.9 

50.2 

125.5 

10 

124.6 

11 

49.4 

124.0 

12 

50.0 

13 

AVE 

17.63 

47.28 

65.74 

99.54 

118.91 

139.63 

163.15 

TEMP 

TIME  BLOCK                                    | 

BLOCK 

7 

8 

9 

10 

11 

12 

0 

1 

2 

3 

215.6 

4 

181.1 

216.0 

5 

185.2 

6 

7 

214.7 

8 

211.1 

311.3 

9 

207.1 

10 

11 

189.6 

205.3 

12 

13 

194.5 

202.7 

AVE 

187.60 

210.36 

311.32 

Table  B-1.  Average  time,  in  hours,  from  sample  introduction. 
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TEMPERATURE  (K) 

TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

0 

0.23 

0.23 

0.23 

1 

0.44 

0.43 

0.38 

0.31 

0.35 

0.33 

2 

0.62 

0.63 

0.72 

0.56 

0.70 

0,63 

0.63 

3 

0.94 

0.92 

0.86 

0.87 

0.90 

0.90 

4 

1.07 

1.11 

1.11 

1.04 

5 

1.36 

1.39 

6 

1.56 

1.55 

7 

1.97 

8 

2.20 

2.13 

2.15 

9 

2.29 

2.36 

2.46 

10 

2.66 

11 

2.79 

2.94 

2.78 

12 

3.13 

13 

3.90 

TEMP 

TIMEB 

LOC 

K 

BLOCK 

7 

8 

9 

10 

11 

12 

AVE 

0 

0.230 

1 

0.370 

2 

0.640 

3 

0.96 

0.906 

4 

1.23 

1.08 

1.099 

5 

1.30 

1.350 

6 

1.560 

7 

1.99 

1.99 

1.980 

8 

2.04 

2.04 

2.00 

2.106 

9 

2.34 

2.34 

2.361 

10 

2.660 

11 

2.78 

2.86 

2.840 

12 

3.132 

13 

3.90 

5.85 

4.874 

Table  B-2.  Sample  temperature,  in  Kelvin,  for  each  time  and  temperature  block. 
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NUMBER 


TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

0 

4 

39 

12 

1 

21 

22 

24 

13 

33 

13 

2 

9 

7 

4 

6 

3 

9 

6 

3 

13 

9 

3 

8 

10 

45 

4 

24 

12 

26 

7 

13 

5 

10 

20 

6 

1 

4 

7 

1 

8 

1 

1 

1 

9 

2 

4 

3 

10 

1 

11 

1 

1 

12 

3 

13 

TEMP 

TIMEBl 

.OCK 

BLOCK 

7 

8 

9 

10 

11 

12 

AVE 

0 

0.230 

1 

0.370 

2 

0.640 

3 

1 

0.906 

4 

12 

4 

1.099 

5 

23 

1.350 

6 

1.560 

7 

1 

1.980 

8 

22 

2.00 

2.106 

9 

3 

2.361 

10 

2.660 

11 

1 

1 

2.840 

12 

3.132 

13 

20 

2 

4.874 

Table  B-3.  Number  of  averaged  line  shapes  for  each  temperature  and  time. 
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Figure  B-1.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
0.23  and  0.37  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and  has  a  slope  of  22 
hours  per  kHz. 
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Figure  B-2.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
0.64  and  0.91  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and  has  a  slope  of  22 
hours  per  kHz. 
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Figure  B-3.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
1.10  and  1.35  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and  has  a  slope  of  22 
hours  per  kHz. 
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Figure  B-4.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
1.98  and  2.1 1  K.  The  baseline  of  each  line  shape  is  ofiFset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and  has  a  slope  of  22 
hours  per  kHz. 
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Figure  B-5.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
2.36  and  2.84  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and  has  a  slope  of  22 
hours  per  kHz. 
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Figure  B-6.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  resonances 
measured  at  two  average  times,  17.6  and  47.3  hours,  after  sample  formation.  The  baseline 
of  each  line  shape  is  offset  by  the  average  sample  temperature. 
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Figure  B-7.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  resonances 
measured  at  two  average  times,  65.7  and  99.5  hours,  after  sample  formation.  The  baseline 
of  each  line  shape  is  offset  by  the  average  sample  temperature. 
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Figure  B-8.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  resonances 
measured  at  two  average  times,  118.9  and  139.6  hours,  after  sample  formation.  The 
baseline  of  each  line  shape  is  offset  by  the  average  sample  temperature. 
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Figure  B-9.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  resonances 
measured  at  two  average  times,  163.1  and  187.6  hours,  after  sample  formation.  The 
baseline  of  each  line  shape  is  offset  by  the  average  sample  temperature. 
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Figure  B-10.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  resonances 
measured  210.4  hours  after  sample  formation.  The  baseline  of  each  line  shape  is  offset  by 
the  average  sample  temperature. 


APPENDIX  C 
DETAILS  FROM  SAMPLE  BNl  193 

This  appendix  presents  selected  data  from  sample  BNl  193.  The  hydrogen  sample 
was  adsorbed  onto  the  BN  at  7AM,  10  November  1993.  The  data  for  this  sample  is 
included  because  it  conveniently  illustrates  some  of  the  interesting  early  time  features  also 
exhibited  by,  but  less  obvious,  in  other  data.  However,  it  is  noted  that  the  background 
signal  was  not  adequately  measured  because  the  experimental  run  ended  prematurely. 


Table  C-1.  Average  time,  in  hours,  from  sample  introduction. 

TIME  (Hours) 


TEMP 

1 

riME  BLOCK                                 1 

BLOCK 

0 

1 

2 

3 

4 

0 

48.1 

1 

41.7 

85.8 

2 

32.7 

75.7 

89.5 

3 

21.5 

57.0 

67.2 

94.9 

4 

21.5 

56.8 

62.3 

101.5 

107.3 

5 

56.0 

102.7 

111.0 

6 

111.9 

7 

54.8 

8 

9 

13.7 

10 

13.1 

11 

12.5 

12 

13 

11.1 

AVE: 

18.01 

52.42 

68.39 

94.41 

110.07 
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153 


Table  C-2.  Temperature,  in  Kelvin,  for  each  time  and  temperature  block. 

TEMPERATURE  (K) 


TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

AVE 

0 

0.208 

0.208 

1 

0.359 

0.490 

0.424 

2 

0.63 

0.575 

0.596 

0.600 

3 

0.883 

0.991 

0.829 

0.897 

0.900 

4 

1.057 

1.093 

1.037 

1.124 

1.197 

1.102 

5 

1.285 

1.269 

1.381 

1.312 

6 

1.286 

1.286 

7 

1.940 

1.940 

8 

9 

2.378 

2.378 

10 

2.663 

2.663 

11 

2.861 

2.861 

12 

13 

4,076 

4.076 

Table  C-3.  Number  of  sweeps  averaged  for  each  block. 

NUMBER 


TEMP 

TIV 

E  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

0 

14 

1 

21 

23 

2 

3 

41 

8 

3 

7 

1 

13 

20 

4 

35 

3 

13 

29 

6 

5 

1 

1 

20 

6 

2 

7 

3 

8 

9 

2 

10 

2 

11 

2 

12 

13 

6 
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Figure  C-1.  E>epcndence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  for 
temperature  blocks  one  and  two.  The  baseline  of  each  line  shape  is  offset  by  the  amount 
of  time  elapsed  since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and 
has  a  slope  of  20  hours  per  kHz. 
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Figure  C-2.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  for 
temperature  blocks  three  and  four.  The  baseline  of  each  line  shape  is  offset  by  the  amount 
of  time  elapsed  since  formation  of  the  sample.  The  dashed  line  is  a  guide  to  the  eye  and 
has  a  slope  of  20  hours  per  kHz. 
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Figure  C-3.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  blocks 
zero  and  one.  The  baseline  of  each  line  shape  is  offset  by  the  average  sample  temperature. 
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Figure  C-4.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  blocks 
two  and  three.  The  baseline  of  each  line  shape  is  offset  by  the  average  sample 
temperature. 


158 


O 
H 

CL. 
IS 

> 

W 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1  1 1 1 1 1 1 1 1  r  1 1 1 1 1 1 1 1 1 1 1 1  1 1 1 1 1 1 1 

Tate  =    110.1   Hours 


-2.4 


-0.4^ 


^ 


I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    r    I    I    I — r~>    r— i—T— r-T— n— r 

150  -50  50  150 

FREQUENCY   (kHz) 


-0.5 


Figure  C-5.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  block 
four.  The  baseline  of  each  line  shape  is  offset  by  the  average  sample  temperature. 


APPENDIX  D 
DETAILS  FROM  SAMPLE  BN294B 

The  appendix  presents  selected  raw  data  from  sample  BN294B.  The  hydrogen 
sample  was  adsorbed  onto  the  BN  at  2:45PM,  22  February  1994.  The  dosing  for  this  data 
set  was  the  most  accurate  of  all  of  the  experimental  runs.  Immediately  before  the  run,  an 
Argon  isotherm  at  77  K  was  performed  in  place  on  the  NMR  sample,  and  the  number  of 
atoms  adsorbed  for  a  final  pressure  of  34  Torr  was  measured.  As  noted  in  the  text, 
previous  isotherms  performed  in  the  isotherm  apparatus  noted  the  liquid  to  solid  transition 
marked  by  a  substep  at  this  point.  For  the  purposes  of  this  thesis,  monolayer  coverage  is 
measured  with  respect  to  this  condition.  Four  layers  of  hydrogen  were  dosed  into  the 
sample  cell  for  sample  BN294b.  It  is  believed  from  the  ease  with  which  the  hydrogen 
sample  passed  through  the  fill  line  capillary  that  virtually  all  of  the  hydrogen  went  into  the 
NMR  cell.  Earlier  hydrogen  samples  sometimes  froze  in  the  fill  line  and  had  to  be  melted 
and  re-admitted.  In  all  cases,  the  hydrogen  dosing  was  followed  by  the  admittance  of 
gaseous  helium  to  determine  whether  the  hydrogen  had  passed  the  fill  line.  In  some  cases, 
a  part  of  the  hydrogen  sample  remained  in  the  fill  line  and  had  the  effect  of  a  plug.  In 
these  cases,  the  NMR  magnetization  was  used  as  a  measure  of  the  amount  of  hydrogen 
present  in  the  NMR  cell.  A  fijrther  purpose  of  admitting  helium  into  the  NMR  cell  was  to 
provide  a  thermal  link  between  the  cold  finger  of  the  dilution  refrigerator  and  the 
individual  boron  nitride  platelets.  Thus,  helium  was  also  admitted  to  the  cell  before  the 
hydrogen  sample  not  only  to  provide  a  thermal  link  but  also  to  check  for  fill  line  flow. 
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Table  D-1.  Time  in  hours  since  formation  of  sample. 


TIME  (Hours) 

TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

51.6 

53.7 

2 

47.4 

55.6 

188.2 

3 

15.9 

45.1 

58.1 

147.9 

163.7 

188.0 

4 

24.3 

38.6 

61.9 

147.2 

156.4 

175.4 

5 

25.6 

66.5 

6 

26.9 

33.3 

70.5 

7 

27.8 

73.7 

78.2 

8 

11.6 

29.1 

83.7 

145.3 

9 

5.9 

90.8 

143.5 

10 

5.5 

94.6 

142.9 

11 

5.1 

96.5 

142.3 

12 

98.9 

141.7 

13 

100.5 

102.2 

134.4 

AVE 

16.51 

40.84 

62.86 

91.89 

102.25 

143.16 

160.05 

183.83 

TEMP 

TIME  BLOCK 

BLOCK 

8 

9 

10 

11 

0 

1 

275.2 

2 

279.2 

3 

205.6 

264.9 

281.6 

4 

204.1 

261.1 

5 

216.1 

264.2 

6 

219.1 

7 

222.9 

8 

228.1 

9 

233.3 

10 

236.7 

11 

237.7 

12 

238.1 

13 

AVE 

211.22 

232.79 

263.38 

278.66 
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Table  D-2.  Sample  temperature,  in  Kelvin,  for  each  time  and  temperature  block. 


TEMPERATURE  (K) 

TEMP 

TIME  BLOCK 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

0 

1 

0.47 

0.47 

2 

0.61 

0.60 

3 

0.91 

0.77 

0.88 

0.95 

0.94 

4 

1.14 

1.09 

1.11 

1.12 

1.08- 

5 

1.33 

1.36 

6 

1.78 

1.59 

1.61 

7 

1.89 

1.87 

1.95 

8 

2.11 

2.12 

2.11 

2.05 

9 

2.34 

2.34 

2.34 

10 

2.63 

2.63 

2.60 

11 

2.83 

2.93 

2.84 

12 

3.13 

3.16 

13 

3.30 

3.87 

4.34 

TEMP 

TIME  BLOCK 

BLOCK 

7 

8 

9 

10 

11 

AVE 

0 

1 

0,39 

0.445 

2 

0.75 

0.59 

0.634 

3 

0.84 

0.95 

0.87 

0.90 

0.891 

4 

1.06 

1.06 

1.08 

1.092 

5 

1.41 

1.32 

1.355 

6 

1.63 

1.65 

7 

1.86 

1.893 

8 

2.15 

2.107 

9 

2.35 

2.343 

10 

2.62 

2.621 

11 

2.86 

2.868 

12 

3.10 

3.129 

13 

3.835 
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Table  D-3.  Number  of  averaged  line  shapes  for  each  time  and  temperature  block. 


NUMBER 


TEMP 

TIME  BLOCKS  (DAYS) 

BLOCK 

0 

1 

2 

3 

4 

5 

6 

0 

1 

9 

5 

2 

7 

6 

3 

35 

5 

10 

1 

56 

4 

7 

35 

15 

2 

12 

5 

4 

15 

6 

3 

1 

9 

7 

3 

10 

9 

8 

10 

2 

29 

8 

9 

2 

20 

2 

10 

1 

4 

2 

11 

2 

6 

2 

12 

9 

2 

13 

2 

10 

30 

TEMP 

TW 

E  BLOCK                1 

BLOCK 

7 

8 

9 

10 

11 

0 

1 

16 

2 

1 

9 

3 

11 

42 

19 

7 

4 

11 

12 

3 

5 

10 

2 

6 

12 

7 

13 

8 

14 

9 

19 

10 

2 

11 

1 

12 

2 

13 
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Figure  D-1.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
0.44  and  0.63  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample. 
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Figure  D-2.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
0.89  and  1.09  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample. 
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Figure  D-3.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
1.36  and  1.89  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample. 
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Figure  D-4.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
2.11  and  2.34  K.  The  baseline  of  each  line  shape  is  oflfset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample. 
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Figure  D-5.  Dependence  of  CW  NMR  derivative  line  shapes  on  ortho-para  conversion  at 
2.62  and  2.87  K.  The  baseline  of  each  line  shape  is  offset  by  the  amount  of  time  elapsed 
since  formation  of  the  sample. 
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Figure  D-6.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  blocks 
zero  and  one.  The  baseline  of  each  line  shape  is  offset  by  the  sample  temperature. 
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Figure  D-7.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  blocks 
two  and  three.  The  baseline  of  each  line  shape  is  offset  by  the  sample  temperature. 
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Figure  D-8.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  blocks 
five  and  six.  The  baseline  of  each  line  shape  is  offset  by  the  sample  temperature. 
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Figure  D-9.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  blocks 
seven  and  eight.  The  baseline  of  each  line  shape  is  offset  by  the  sample  temperature. 
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Figure  D-10.     Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time 
blocks  nine  and  ten.  The  baseline  of  each  line  shape  is  offset  by  the  sample  temperature. 
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Figure  D-1 1.  Temperature  dependence  of  CW  NMR  derivative  line  shapes  for  time  block 
eleven.  The  baseline  of  each  line  shape  is  offset  by  the  sample  temperature. 
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